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FOR  AN  UNMANNED  RESEARCH  VEHICLE 
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Preface 


The  purpose  of  this  study  wejs  to  build  a  robust  flight  control  system  for 
Lambda,  an  unmanned  research  vehicle.  My  interest  in  unmanned  aerial  vehicles 
led  me  to  choose  this  topic  and  I  am  thankful  to  have  had  the  opportunity.  Al¬ 
though  many  controllers  have  been  synthesized  for  Lambda,  this  is  the  first  design 
applying  Quantitative  Feedback  Theory  to  a  mathematical  model  developed  from 
Lambda’s  fight  test  data  and  based  on  Lambda’s  entire  flight  envelope.  Simula¬ 
tions  show  encouraging  performance  and  I  hope  to  see  my  controller  implemented 
on  Lambda  in  the  near  future. 

My  sincere  thanks  go  to  the  many  people  who  helped  me  during  the  research. 
First,  I’d  like  to  thank  my  fellow  classmates  in  the  Guidance  and  Control  Section 
for  their  help,  advice  and  support.  I’ve  learned  and  benefited  from  every  one  of 
them  and  wish  them  all  the  best  of  luck  in  their  Air  Force  careers.  I’d  also  like  to 
thauk  those  students  who  worked  on  the  project  previously  and  had  a  great  impact 
on  this  effort.  In  particular.  I’d  like  to  thank  Capt  David  Wheaton  for  proving  the 
concepts  and  sharing  his  results  with  me,  and  First  Lieuteneint  Gerald  Swift  for 
developing  the  mathematical  model  from  flight  test  data.  My  thanks  also  to  the 
Control  Systems  Development  Branch  of  the  Wright  Laboratory  and  especially  to 
Capt  Stuart  Sheldon  for  providing  the  thesis  topic  and  giving  welcome  advice  and 
direction.  It  has  been  a  pleasure  working  with  Professor  Constantine  Houpis  and  I’ve 
enjoyed  the  unique  opportunity  to  work  with  Professor  Isaac  Horowitz  employing  his 
Quantitative  Feedback  Theory.  Finally,  I’d  like  to  thank  my  wife,  Nancy  and  my 
son  for  their  understanding  and  patience.  I  couldn’t  have  done  it  without  them. 

Donald  Jerry  Lacey  Jr. 
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Abstract 

This  thesis  describes  the  application  of  the  digital  multiple-input  multiple- 
output  (MIMO)  Quantitative  Feedback  Theory  (QFT)  technique  to  the  design  of 
a  three  axis  rate  controller  for  the  Lambda  Unmanned  Research  Vehicle.  The  re¬ 
sulting  robust  controller  performs  well  throughout  the  flight  envelope  without  gain 
scheduling.  It  results  from  research  done  at  the  Air  Force  Institute  of  Technology  . 
The  Lambda  URV  was  designed  by  the  Control  Systems  Development  Branch  of  the 
Wright  Laboratory  at  Wright- Patterson  AFB,  OH  for  flight  testing  aircraft  control 
hardware  and  software.  The  flight  control  system  is  built  using  a  small  perturbation 
linear  model  developed  from  flight  test  data.  The  actuators,  also  modelled  from 
aircraft  test  data,  are  second  order  in  roll  and  pitch  and  first  order  in  yaw.  Nineteen 
separate  plants  are  used  to  represent  the  flight  envelope  of  the  aircraft  resulting  from 
variations  in  speed,  altitude,  center  of  gravity  location,  and  weight.  The  sample  rate 
is  50  Hertz.  The  pitch  channel  is  decoupled  from  the  lateral-directional  channel  re¬ 
sulting  in  a  single-input  single-output  (SISO)  system  for  the  pitch  channel  and  a 
two-by-two  MIMO  system  for  the  lateral-directional  channel.  The  design  employs 
the  Nichols  Chart  and  is  accomplished  in  the  w'-domain. 
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A  ROBUST  DIGITAL  FLIGHT  CONTROL  SYSTEM 
FOR  AN  UNMANNED  RESEARCH  VEHICLE 
USING  DISCRETE  QUANTITATIVE 
FEEDBACK  THEORY 


I.  Introduction 


1,1  Background 

The  Control  Systems  Development  Branch  of  the  Wright  Laboratory  at  Wright- 
Patterson  AFB,  OH  (WL/FIGL)  developed  the  Lambda  unmanned  research  vehicle 
(URV)  for  in-flight  testing  of  aircraft  control  hardware  and  software.  The  Wright 
Laboratory,  at  the  request  of  Congress,  is  using  Lambda  to  develop  and  test  the 
common  core  avionics  package  destined  to  be  used  on  all  unmanned  aerial  vehicles 
(UAVs)  in  the  Department  of  Defense  [25].  In  order  to  provide  the  necessary  test 
data.  Lambda  requires  a  stable,  robust  flight  control  system  which  performs  well 
throughout  the  entire  flight  envelope.  Capt  David  Wheaton,  in  his  masters  thesis, 
designed  a  flight  control  system  for  Lambda  using  digital  Quantitative  Feedback 
Theory  [26].  The  Wright  Laboratory  is  pleased  with  the  design;  however,  it  is  based 
on  an  earlier  mathematical  model  of  Lambda  derived  solely  tVom  the  physical  di¬ 
mensions  of  the  aircraft  design.  Flight  test  data  indicates  that  the  aircraft  behaves 
much  differently  than  the  physically  derived  predictions.  A  new,  more  representative 
model,  derived  from  night  test  data,  was  created  by  Lt  Gerald  Swift  [24]  and  is  used 
for  this  design.  A  new  flight  control  system  providing  the  robustness  of  Wheaton’s 
design  but  using  an  updated  model  is  needed.  The  Wright  Laboratory  is  sponsoring 
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this  thesis  because  they  require  a  robust  controller  for  Lambda.  Using  Quantitative 
Feedback  Theory  in  the  design  will  provide  the  robustness  required. 

Quantitative  Feedback  Theory  (QFT)  was  developed  by  Dr  Horowitz  in  1979. 
It  is  an  extension  of  his  work,  which  began  in  1958  [10],  concerning  the  control  of 
plants  with  uncertain  paramet  ers.  QFT  is  a  frequency  design  technique  incorporating 
plant  uncertainty  early  in  the  design  process  resulting  in  a  robust  controller  [8]. 
Robust  automatic  flight  control  systems  in  a  URV,  and  UAVs  in  general,  reduce 
pilot  work  load,  increaise  safety,  and  aid  in  the  recovery  of  the  vehicle.  The  automatic 
flight  control  system  gives  the  ground-based  pilot  the  ability  to  hold  a  desired  flight 
condition.  Additionally,  functions  normally  provided  by  an  on-board  pilot,  such  as 
turn  coordination,  can  be  performed  by  the  automatic  flight  control  system. 

1.2  Problem 

The  challenge  is  to  produce  a  robust  digital  flight  control  system  for  Lambda 
using  QFT  and  deliver  a  set  of  controller  difference  equations  to  the  sponsor.  The 
sponsor,  the  Wright  Laboratory,  is  responsible  for  installing  and  testing  the  system 
on  Lambda.  The  system  must  uncouple  the  roll,  pitch  and  yaw  responses.  For 
example,  a  command  to  roll  the  airplane  should  not  produce  t  yaw  or  pitch  response. 
The  system  must  meet  all  performance  specifications  throughout  the  flight  envelope. 
The  flight  envelope  includes  speeds  from  45  to  110  knots,  center  of  gravity  locations 
between  21.8  and  32.4%  Mean  Aerodynamic  Cord  (MAC),  vehicle  weights  from  181 
to  215  pounds,  and  altitudes  up  to  5000  feet.  Speed  will  vary  between  76  and  185 
feet-per-second.  The  system  must  perform  several  autopilot  functions  including: 
attitude  hold,  altitude  hold,  heading  hold,  and  out-of-sight  maneuvering. 

1.3  Summary  of  Current  Knowledge 

Quantitative  Feedback  Theory  is  a  design  technique  aimed  at  meeting  both  sta¬ 
bility  and  performance  specifications  [20]  and  guarantees  system  performance  over  a 
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wide  range  of  uncertainties  in  the  plant  and  environment  [26].  In  1963,  Dr  Horowitz 
first  pioposed  a  Nichols  chart  design  technique  that  accounted  for  uncertainties  “up 
front,”  early  in  the  design.  In  1979,  after  years  of  work,  he  refined  the  technique  and 
coined  the  name  Quantitative  Feedback  Theory  or  QFT  [26],  which  later  became 
known  as  Method  One.  From  the  beginning,  QFT  was  criticized  for  its  tendency  for 
overdesign  and  less  than  optimum  realization.  However.  D  '  Horowitz  jicver  claimed 
his  technique  produced  optimum  designs.  He  simply  claimed,  “The  technique  guar¬ 
antees  a  satisfactory  design  ...  by  a  comparatively  straightforward  systematic  pro¬ 
cedure  [8].”  In  fact,  QFT  is  one  of  the  very  few  methods  claiming  a  straightforward 
approach  that  addresses  robust  performance  [20].  In  1982,  Dr  Horowitz  developed 
the  QFT  Improved  Method,  called  Method  Two,  which  reduced  the  overdesign  [9]. 
According  to  Horowitz  [10],  plant  and  disturbance  uncertainties  are  the  reasons  for 
feedback  in  the  system.  This  philosophy  is  fundamental  to  QFT  [8].  Quantitative 
Feedback  Theory,  a  development  of  conventional  control  theoiy,  is  growing  in  ac¬ 
ceptance  in  spite  of  criticism  from  mod-»rn  control  theorists.  Modern  control  theory, 
or  optimal  control  theory,  completely  ignored  plant  uncertainties  until  recently  [27]. 
Additionally,  history  has  shown  QFT  designs  to  succeed  where  other  techniques  have 
failed.  Specifically,  QFT  designers  working  on  the  F-4B  were  able  to  produce  a  flight 
control  system  with  excellent  stability  margins  using  a  single  fixed  gain  compensator. 
The  original  compensator  used  a  gain  scheduling  technique  to  handle  variations  in 
the  plant  and  en»^ironment  but  did  not  work  well  due  to  its  additional  complexity 
and  inherent  computational  lag  [27].  Successful  QFT  controllers  have  been  synthe¬ 
sized  for  the  YF16CCV,  Cl  35,  AFTi-  FfG.  STOL,  .^nd  F4B  aircraft.  “There  is  a 
great  lack  of  similar  design  examples  by  other  modern  robustness  method[s].  ,  .But 
where  quantitative  comparison  was  possible,  QFT  design  was  much  more  economical 
in  loop  compensation  gam  and  bandwidth  [27].”  Air  Force  Institute  of  Technology 
(  AITT)  students  have  been  designing  flight  control  systems  using  QFT  since  1982. 
The  following  is  a  partial  list  of  theses  using  QFT: 
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•  J.  M.  Adams  —  Digital  QFT  Design  for  the  AFTI/F- 16  [1]  (1988) 

•  P.  B.  Arnold  —  FCS  Reconfiguration  Using  QFT  [2]  (1894) 

•  B.  T.  Clough  —  Reconfiguration  for  a  STOL  Aircraft  Using  QFT  [4]  (1985) 

•  J.  S.  Coucouies  —  Effects  '  f  )i.\.rHizing  QFT  Designs  [5]  (1985) 

•  S.  W.  Hamilton  —  Dl  v  QF'. '' .gn  for  an  Unmanned  Research  Vehicle  [o] 
(1987) 

•  B.  S.  Migyanko  —  Integrate*',.  F’-jit/Propnlsion  Control  for  a  STOL  Aircraft 
Using  QFT  [15]  (1986) 

•  K.  N.  Neumann  —  Digital  Ql  T  Design  jar  the  Control  Reconfigurable  Combat 
Aircraft  [19]  (1988) 

•  P.  T.  Ott  —  URV  Reconfiguration  using  Digital  QFT  [21]  (1988) 

•  H.,  H.  Russell  —  Analog  QFT  Design  for  the  KC-135  [22]  (1984) 

•  D.  L.  Schneider  —  AFTI/F-16  FCS  Design  Using  Digital  QFT  [23]  (1986) 

•  D.  G.  Wheaton  —  AFCS  Design  for  a  URV  Using  Discrete  QFT  [26]  (1990) 

1.4  Assumptions 

This  thesis  uses  a  linear  time-invariant  (LTl)  model  of  Lambda  which  is  known 
to  be  a  nonlinear  system.  Further,  the  equations  of  motion  and  the  aircraft  dynamics 
are  Laplace  transformable.  In  addition,  the  following  conditions  are  assumed: 

•  Small-angle  perturbation  models  are  valid. 

•  Aircraft  mciss  is  constant. 

•  The  aircraft  is  rigid:  that  is,  no  bending  or  ilutter  modes. 

•  The  earth  is  an  inei  Lial  reference  frame. 

•  The  atmosphere  is  fixed  in.  relation  to  the  earth. 

•  The  commanded  inputs  and  the  commanded  outputs  are  measurable. 
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•  Three-axis  rate  signals — pitch  rate  (9),  roll  rate  (p)  and  yaw  rate(r) — and  the 
Euler  angles —  pitch  angle  (0),  roll  angle  (</>),  and  yaw  angle  (^) — are  available 
on  the  Lambda  DRV. 

•  A  digital  sampling  rate  of  50  Hz  is  I'sed. 

1.5  Scope 

The  design  is  limited  to  the  development  of  the  set  of  c^"'.troller  difference 
equations  needed  to  control  Lambda  throughout  the  expected  flight  envelope.  T  he 
Wright  Laboratory  will  implement  the  design  on  Lambd  ..  Although  the  success  of 
the  design  will  be  proven  when  the  automatic  flight  cou  vol  system  is  installed  and 
flowii  on  Lambda,  this  design  ends  with  a  successful  computer  sunulation. 

1.6  Standards 

The  controller  is  required  to  have  a  45°  phase  margin  throughout  the  flight 
envelope  in  all  three  axes — pitch,  roll,  and  yaw.  In  addition,  the  aircraft  is  required 
to  meet  figures-of-merit  including  specific  rise  time  and  overshoot  requirements.  See 
Table  (1.1).  The  perfo:  ance  requirements  are  the  same  as  in  the  previous  work 
done  by  Wheaton  [26].  Frequency  response  tracking  bounds  for  pit^h,  roll,  yaw  are 
included  in  Figure  (1.1).  Time  response  tracking  bounds  are  included  in  Figure  (1.2) 


1.7  Approach/Methodology 

This  design  begins  with  a  mathematical  description  of  the  airplane,  followed 
by  definition  of  the  flight  envelope  and  flight  conditions.  Individual  flight  conditions 
include  variations  in  flight  speed,  center  of  gravity  location,  vehicle  weight,  and  alti¬ 
tude.  Analysis  of  the  model  reveals  the  lateral-directional  dynamics  can  be  separated 
from  the  longitudinal  dynamics.  Since  separation  is  possible,  the  design  reduces  to  a 
single-input  single-output  (SISO)  QFT  system  design  and  a.  two-by-two  multiple- 
input  multiple-output  (MIMO)  QFT  system  design.  Analysis  of  the  system  indicates 
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Caiii(ilB}  CiinfdB)  Caln(ilB) 


Figure  1.2.  Pitch,  Roll,  and  Yaw  Time  Response  Tracking  Bounds 
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Table  1.1.  Figures-of-Merit  for  the  QFT  Response  Models 


Model 

Ts  (sec) 

7cmd;.n(^) 

0.84 

1.56 

1.0 

0.15 

2.05 

1.25 

PLT)[s) 

Pctndfai^} 

0.87 

1.56 

1.0 

pubU) 
PcmdfT  ai^) 

0.44 

0.78 

1.0 

r(*) 

3.48 

6.22 

1.0 

ruBKo) 
^cmdri  n(^) 

1-76 

3.13 

1.0 

that  unstable  and  nonminimum  phase  conditions  exist.  In  contrast  with  Wheaton’s 
work,  a  weighting  matrix  to  blend  the  multiple  inputs  into  three  effective  inputs  [26] 
is  already  included  in  the  mathematical  model.  Since  inversion  of  the  plant  matrix  is 
necessary  [8],  the  plants  must  be  square,  (the  same  number  of  outputs  as  inputs)  in 
this  case,  three-by-three.  Flight  data,  used  for  the  mathematical  model,  produced 
square  plants. 

A  three-by-three  MIMO  system  QFT  design  requires  the  synthesis  of  three 
separate  transmission  loops;  a  two-by-two  plant  requires  two  loops.  Design  starts 
with  the  loop  with  the  smallest  bandwidth.  Once  obtained,  the  first  loop  is  used 
in  developing  the  second  and  third  loops.  The  technique  of  designing  loops  using 
elements  of  previously  designed  loops  is  part  of  Dr  Horowitz’s  improved  method 
(Method  Two).  It  reduces  the  inherent  overdesign,  a  characteristic  of  the  QFT 
method  [9].  Once  the  three  loops  are  shaped,  the  design  is  simulated  on  the  computer 
to  validate  the  performance  expectations.  In  designing  the  controller  for  the  SISO 
system,  the  steps  u.sed  for  application  of  the  discrete  QFT  technique  are;, 

•  Choose  the  Flight  Conditions 

•  Determine  the  Plant  Transfer  Functions,  P,(s) 

•  Transform  P,(s)  to  the  w'  Domain,  Pi(w') 


•  Determine  the  Frequency  Response  Data 

•  Extract  the  Template  Data  from  the  Frequency  Response  Data 

•  Plot  the  Templates 

•  Choose  a  Nominal  Plant,  P^ 

•  Use  the  Templates  to  form  the  Stability  and  Performance  Bounds 

•  Shape  the  Nominal  Loop,  Lo(w') 

•  Extraci  the  Compensator,  G(w'),  from  Lo(w') 

•  Synthesize  the  Prefilter,  F(w') 

•  Verify  Loop  Shaping  in  the  w'-Plane 

•  Transform  F(w')  to  F(z)  and  G(w')  to  G(z) 

•  Simulate  using  G(z),  F(z),  and  P,-(s) 

A  slightly  modified  procedure  must  be  followed  when  designing  the  controller  for  the 
MIMO  sysi 'm  using  the  discrete  QFT  technique.  The  steps  are: 


•  Choose  the  Flight  Conditions 


•  Determine  the  Plant  Transfer  Function  Matrix,  P(s) 


•  Invert  P(s) 


P-^(S) 


Pn('®)  PnC'S) 
Pni^)  P2i(«) 


•  Calculate  the  Q(s)  Matrix 


Q(s)  = 


1 

pJiW 

1 

pLI*) 

<7ii 

912 

p«(*)  - 

<721 

922 

(1.1) 


(1.2) 
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•  Transform  Q(s)  to  the  w'  Domain,  Q(w')  which  includes  the  necessary  zero 
order  hold  (ZOH). 

•  Determine  the  Frequency  Response  Data 

•  Extract  the  Template  Data  from  the  Frequency  Response  Data 

•  Plot  the  Templates  for  One  Input-Output  Pair 

•  Choose  a  Nominal  Plant 

•  Use  the  Templates  to  Form  the  Stability  and  Performance  Bounds 

•  Shape  the  Nominal  Loop,  Lo(w') 

•  Extract  the  Compensator,  G(w'),  from  Lo(w') 

•  Synthesize  the  Prefilter,  F(w') 

•  Shape  the  Remaining  Loops 

•  Form  the  Remaining  G(w')s  and  F(w')s 

•  Verify  Loop  Shaping  in  the  w'-Plane 

•  Transform  the  F(w')3  to  F(z)s  and  the  G(w')s  to  G(z)s 

•  Simulate  using  the  G(z)s,  F(z)s,  and  P(s)s 

The  next  step  adds  the  automatic  flight  control  modes.  The  first  mode  is  turn 
coordination.  There  are  at  least  four  methods  of  turn  coordination  [3].  One  is  chosen 
based  on  the  sensors  available  and  the  desires  of  the  Wright  Laboratory.  Wheaton 
added  and  tested  a  wings  leveler  to  his  design  [26].  This  effort  extends  the  design 
by  allowing  the  pilot  to  set  a  bank  angle.  Auto  speed  and  altitude  hold  modes  are 
the  final  modes  added  to  Lambda.  Finally,  the  system  is  simulated  on  M.ATRIXx- 

1.8  Materials  and  Equipment 

The  design  and  analysis  is  performed  using  M  ATRIXx  on  the  SP.ARC  station  2 
computers  available  at  the  Air  Force  Institute  of  Technology  and  at  the  Wright 
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Laboratory.  MATRIX^  is  a  matrix  calculator  with  graphics  and  a  system  build 
capability  which  can  be  used  to  simulate  the  complete  system  [14], 


1.9  Other  Support 

The  Wright  Laboratory,  owners  of  Lambda  and  sponsor  of  the  thesis,  provided 
performance  specifications  and  aircraft  dynamics  models.  The  thesis  committee,  con¬ 
sisting  of  Dr  Houpis,  Dr  Horowitz,  Dr  E.ving,  and  Capt  Sheldon,  provided  direction, 
advice,  and  constructive  feedback  during  and  after  the  design  effort. 


1-11 


11.  The  Aircraft  Model 


Figure  2.1.  Lambda  URV 


This  chapter  discusses  the  Lambda  aircraft  and  the  model  used  in  the  design 
of  the  QFT  rate  controller. 

2.1  The  Aircraft 

Lambda,  shown  in  Figure  2.1,  is  a  .small  remotely-piloted  airplane  with  a  14 
foot  wing  span  and  weight  of  approximately  200  pounds.  See  Table  2.1  for  Lambda’s 
characteristics.  It  uses  a  pusher  propeller  behind  the  fuselage  and  in  front  of  a 
conventional  aft  tail.  The  horizontal  tail  consists  of  a  horizontal  stabilizer  and  a 
split  elevator.  A  vertical  tail  is  located  on  either  end  of  the  horizontal  tail,  and 
consists  of  a  vertical  stabilizer  and  rudder.  The  wings  are  slightly  tapered;  each 
has  three  trailing,  moveable  control  surfaces.  Lambda  closely  resembles  the  Pioneer 
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Figure  2.2.-  Three  Views  of  the  Lambda  URV 


UAV  successfully  employed  by  the  United  States  Navy  in  the  Persian  Gulf  War.  The 
Pioneers  are  slightly  larger  and  twice  as  heavy  but  share  the  configuration  of  an  aft 
tail,  center  mounted  engine,  and  pusher  propeller.  Lambda  has  ten  separate  control 
surfaces.  They  are: 

•  Left  Elevator 

•  Right  Elevator 


•  Left  Rudder 


Table  2.1.  Lambda  Descriptive  Data 


Dimensions  Wing  Span .  14  ft 

Wing  Area . .  19  sq  ft 

Length .  9.6  ft 

Height  With  Landing  Gear .  3  ft 

_ _ Propeller  Diameter .  2.3  ft _ 

Weights  Maximum  Fuel . . .  14  lb 

Maximum  Payload  .  15  lb 

_ Maximum  Flight  Weight  .  200  lb _ 

Performance  Maximum  Level  Speed  at  Sea  Level  115  mph  (100  knots) 

Stall  Speed  .  63  mph  (55  knots) 

_ Stall  Speed  With  Flaps .  52  mph  (45  knots) 

Engine  Power  .  18  hp 

Type  .  2  cycle,  2  cylinder 

Control  Limits  Elevator  Deflection  Limit .  ±15° 

Rudder  Deflection  Limit .  ±25° 

Flap  Deflection  Limit  .  20°  down 

Aileron  Deflection  Limit .  u:  15° 


•  Right  Rudder 

•  Left  Aileron 

•  Right  Aileron 

•  Left  Outer  Flap 

•  Right  Outer  Flap 

•  Left  Inner  Flap 

•  Right  Inner  Flap 


Even  though  each  control  surface  can  be  operated  independently,  the  flight  data,  used 
to  create  the  mathematical  model  upon  which  this  design  is  based,  were  obtained 
using  the  control  surfaces  together.  All  four  of  the  flap  control  surfaces  were  operated 
together  as  were  the  two  rudders.  The  split  elevators  were  deflected  together  while 
the  ailerons  were  operated  differentially.  Much  of  the  flexibility  built  into  Lambda 
was  discarded  by  employing  the  surfaces  in  this  manner.  Future  projects  should 
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consider  implementing  control  schemes  that  rely  on  independent  control  surface  de¬ 
flections.  Since  the  mathematical  model  does  not  include  flap  actuator  dynamics,  the 
flaps  are  not  used  as  a  control  surface.  Again,  a  future  design  might  incorporate  the 
flaps  as  an  active  part  of  the  flight  control  system  once  the  flap  actuator  dynamics 
are  known  and  included  in  the  model. 

This  control  system  is  based  on  the  design  specifications  of  Lambda.  It  is  ex¬ 
pected  that  Lambda  will  gain  weight.  Indeed,  today’s  model  of  Lambda  is  heavier 
than  it  was  when  first  built,  which  was  heavier  than  specified.  The  weight  increase  re¬ 
sults  from  adding  equipment,  such  as  video  cameras  and  additional  sensors.  For  this 
reason,  this  design  includes  weights  up  to  215  pounds.  Additionally,  since  Lambda 
is  designed  to  fly  at  speeds  up  to  100  knots,  speeds  up  to  110  knots  are  included 
to  encompass  all  of  the  expected  plant  variation.  The  heavier  weights  are  very  de¬ 
manding  of  the  technique  and  contribute  greater  plant  variation  compared  with  the 
higher  speeds  which  contribute  little  variation. 

2.2  The  Model 

Lt  Gerald  Swift  developed  the  small  perturbation  model  based  on  a  series  of 
test  flights  conducted  by  the  Wright  Laboratory  [24].  His  model  incorporates  the 
traditional  stability  derivatives  into  a  state  space  representation  of  the  aircraft.  Lt 
Swift  provided  a  MATLAB  macro  file  which  creates  the  state  space  representation 
along  with  the  stability  derivatives  when  given  a  set  of  flight  conditions  (see  Appendix 
A).  The  aircraft,  excluding  actuator  and  sensor  dynamics,  can  be  described  using 
nine  states.  They  are:  q,  p,  r,  9,  <{)  yij),  a,  u,  /3.  Of  these,  q,  9,  a,  and  u  are  used 
to  describe  the  longitudinal  channel,  while  the  states  p,  r,  <^,  0,  and  ^  describe  the 
lateral-directional  dynamics.  The  aircraft  model  state  space  representation  is: 

x(<)  =  A(t)x(<)  -f  B(t)u(<)  (2.1) 
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y(0  =  Cx(i)  (2.2) 

Assuming  zero  initial  conditions  (xq  =  0)  and  a  trim  condition  (xq  =  0),  the  Laplace 
transform  of  the  plant  is  the  perturbation  model: 


X(s)  =  A{s)X{s)  +  B{s)U(s) 


(2.3) 


where: 


Y(5)  =  C{s)X{s) 


A(s)  = 


A(s)/on5  0 

0  A(s),a< 


B(s)  = 


B(s)/onff 

B(3)/ot 


C{s) 


C(s)/on5 


(2.4) 


(2.5) 

(2.6) 


The  Lambda  longitudinal  and  lateral-directional  dynamics  are  modelled  as  decou¬ 
pled.  In  other  words,  pitch  motion  doesn’t  induce  a  roll  or  yaw  response  and  a  yaw 
or  roll  does  not  induce  a  pitch  response.  While  the  model  has  decoupled  dynamics, 
actually  ,  there  is  some  unmodelled  coupling.  However,  this  design  uses  the  best 
model  available  for  Lambda. 


S.3  Actuators  and  Sensors 

The  actuators  used  in  the  mathematical  model  of  Lambda  are  second-order  in 
the  pitch  and  roll  channels  and  first-order  in  the  yaw  channel.  The  model  used  for 
the  pitch  channel  actuator  is: 


87 

52  -f  18s  +  120 

The  model  used  for  the  roll  channel  is: 

84 

s2  -I-  18s  -f- 120 


(2.7) 


(2.8) 
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The  model  used  for  the  yaw  chaunel  is: 

5.58 
s  +  5.2 


(2,9) 


The  sensors  used  to  measure  the  rates  are  identical  gyroscopes  with  transfer  functions 
modelled  as: 


50 

5  +  50 


(2.10) 


2.4  Center  of  Gravity 

Lt  Swift  also  provided  the  formula  for  calculating  the  center  of  gravity  (eg) 
location: 

(2-11) 

where  Xcg  is  the  location  of  the  center  of  gravity  in  inches  from  the  nose  of  the 
aircraft. 

2.5  Flight  Envelope 

The  Wright  Laboratory  built  Lambda  to  fly  within  the  flight  envelope: 

•  Speeds  from  45  to  100  Knots 

•  Maximum  Flight  Weight  200  Pounds 

•  Altitude  up  to  5000  Feet 
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III.  Longitudinal  SISO  Design 


This  chapter  covers  the  design  of  the  longitudinal  channel  of  Lambda  using 
QFT.  The  design  steps  are  outlined  in  Section  1.7  of  Chapter  1, 


3.1  Flight  Conditions 

The  flight  conditions  are  selected  from  the  design  flight  envelope  of  Lambda 
with  considerations  given  to  the  growth  of  the  program  and  expected  changes.  There 
are  four  variables  considered  in  the  selection  of  flight  conditions.  They  are: 

•  Center  of  Gravity  Location 

•  Flight  Speed 

•  Altitude 

•  Aircraft  Weight 

The  minimum  and  maximum  values  of  each  variable  create  a  set  of  sixteen  flight 
conditions.  The  sixteen  flight  conditions  produce  two  widely  spaeed  plant  groupings. 
The  difference  in  speed  contributes  the  most  to  the  wide  plant  variation,  with  the 
other  variables  contributing  smaller  variations.  Three  additional  speeds  are  added 
to  fill  in  the  gap  between  the  lower  and  higher  speeds.  A  weight  of  200  pounds,  a 
center  of  gravity  location  in  the  center  of  the  range,  along  with  the  lowest  altitude, 
are  chosen  for  the  additional  conditions.  A  total  of  nineteen  separate  plants  are  used. 
See  Table  3.1  for  the  conditions  used  in  the  design. 
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Table  3.1.  The  Aircraft  Flight  Conditions 


Plant  # 

eg 

(%  MAC) 

Qoo 

(lbs/ft2) 

Weight 

(lbs) 

1 

21.8 

no 

35.05 

181 

2 

21.8 

no 

35.05 

215 

3 

21.8 

no 

40.08 

4 

21.8 

no 

40.08 

imi 

5 

32.4 

no 

35.05 

181 

6 

32.4 

no 

35.05 

215 

7 

32.4 

no 

40.08 

181 

8 

32.4 

no 

40.08 

215 

9 

27.1 

61 

12.61 

200 

10 

27.1 

77 

20.08 

200 

11 

27.1 

93 

29.90 

200 

12 

32.4 

45 

5.915 

181 

13 

32.4 

45 

5.915 

215 

14 

32.4 

45 

6.765 

181 

15 

32.4 

45 

6.765 

215 

16 

21.8 

45 

5.915 

181 

17 

21.8 

45 

5.915 

215 

18 

21.8 

45 

6.765 

181 

19 

21.8 

45 

6.765 

215 

3.2  s-Plane  Transfer  Fu.  ctions 

The  longitudinal  plant  transfer  functions  are  developed  from  MATRIX^  us¬ 
ing  the  MATRIX^  System  Build  capability.  Each  condition  is  used  as  an  input 
into  the  MATLAB  macro  file  provided  to  the  Wright  Laboratory  by  Lt  Swift  (see 
Appendix  A).  The  outputs  include  A  and  B  matrices  for  the  longitudinal  and 
lateral-directional  dynamics  of  Lambda.  The  A  and  B  matrices  for  the  longitudinil 
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channel  are  of  the  form: 


X„ 

Zu 

u-Z* 

Mu  Zu 
u-Za 

0 


X. 

z^ 

U-Za 

M»fMAZ» 

U-Za 

0 


-32.174  cos(0) 
-32.174  sin  0 
-32.174M<i  sin0 
0 


(3.1) 


The  state  space  representation  is  used  in  a  MATRIX;^  System  Build  simulation  to 
produce  a  linearized  state  space  model  which  includes  the  actuators  and  sensors.  The 
resulting  state  space  representation  is  then  transformed  to  a  transfer  function  which 
relates  the  output  to  the  input.  The  transfer  functions  are  used  in  the  design  of  the 
QFT  controller  which  includes  a  compensator  as  well  as  a  prefilter  (see  Figure  3.1). 
The  following  transfer  function  represents  the  longitudinal  channel  at  the  nominal 
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Figure  3.1.  QFT  Controller  Block  Diagram 


conditions: 

qjs)  _ 1.5774  X  10S(5)(5  -f  0.054837)(s  +  3.7340) _ 

«eJev(s)  («  +  0.0017819 ±  i0.22985)(s  +  3.2050  ±i8.0129)(s  +  9  ±  j6.2450)(«  +  50.000) 

A  complete  set  of  transfer  functions  used  in  the  design  is  included  in  Appendix  B. 
3.3  V  Plant  Transfer  Rinctions 

Lambda  uses  digital  flight  control  hardware  requiring  a  digital  (discrete)  flight 
control  system  design.  Although  an  additional  transform  is  necessary,  the  w'-plane 
is  used  in  this  design  because  it  allows  use  of  s-plane  design  tools,  such  as  the  Nichols 
Chart  in  the  design  of  a  discrete  controller.  Other  options  include  designing  in  the 
z-domain  or  designing  in  the  s-domain  and  then  discretizing  the  results.  Designing 
directly  in  the  z-domain  is  tedious  and  very  sensitive  to  small  numerical  errors  be¬ 
cause  the  entire  left  half  of  the  s-plane  is  compressed  into  a  unit  circle.  Designing 
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either  in  the  w'-plane  or  the  s-plane,  the  psuedo-continuous  approach  (PCX),  in¬ 
volves  the  property  of  warping  when  transforming  the  s  or  w'-plane  controller  into 
the  z-domain.  The  degree  of  “warping”  decreases  as  the  sampling  time  decreases 
as  shown  in  Chapter  VI  of  reference  [13].  A  good  degree  of  correlation  between 
the  w'-  and  z-domain  is  maintained  for  a  longer  sample  time  than  for  the  s-  to 
z-domain  correlation.  Thus,  for  these  reasons,  the  w'-plane  is  an  excellent  compro¬ 
mise.  The  w'-plane  representation  is  normally  obtained  by  transforming  the  s-plane 
plant,  with  a  zero  order  hold  (ZOH),  into  the  z-plane  then  performing  the  bilinear 
transformation: 


/  2  12-1] 

T  L2-f  ij 

(3.6) 

2  +  w'T 

2-w'T 

(3.7) 

where  T  is  equal  to  the  sampling  interval.  The  sampling  interval 

is  obtained  from: 

T-  1 

Sampling  Frequency 

(3.8) 

Lambda  has  a  sampling  frequency  of  50  Hz  which  is  equivalent  to  < 

a.  sampling  interval 

of  0.02  seconds. 

As  a  review,  the  following  relationships  between  s,  z,  and  the  w'-plane  are 

presented  [13]: 

5  =  O's^plane  "b  j^s-plane  —  n  -b  jiO 

(3.9) 

>  ,  •  •  2  r^-i] 

W  —  (T^i  plane  "b  plane  —  U  JV  — 

i  IZ  -f  L. 

(3.10) 

■  Tw'h-2  1 
_Tw'-|-2j 

(3.11) 

(3.12) 

2  =  =1  0  1  LuT 

(3.13) 
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plcne^l 


(3.14) 


y.-p/anJ^  <0,297 


are  satisfied  then: 


where: 


5  «  w 


5  —  O' s— plane  i  jit^a-plane 


(3.15) 


(3.16) 


(3.17) 


and  (Ta-piane  and  tjJa-piane  are  the  components  of  the  fastest  pole  or  zero  in  the  plant. 
The  fastest  pole  in  the  longitudinal  plants  is  the  pole  associated  with  the  sensor  lo¬ 
cated  at  —50.  The  highest  frequency  in  Equation  (3.5)  is  associated  with  the  complex 
pole  pair  at  (s  -f-  3.205  ±j8.013).  Using  cta-pUne  —  ~50  and  T=0.02  Equation  (3.14) 


becomes: 


^<^s-planeT  ^  _  (— 50)(0.02)  _  Q  25  <  2 

L  2  .  2 


(3.18) 


Using  Wa-p/ane  =  8.01  and  T=0.02  Equation  (3.15)  becomes: 


UJa-planeT  (8.01)(0.02) 


=  0.0801  <  0.297 


(3.19) 


Therefore  the  condition  of  Equation  (3.15)  is  definitely  satisfied.  In  Equation  (3.14), 
the  “much  less”  requirement  may  not  be  met.  The  ratio  of  2  to  I  = 

|lz52M2:22iJ^  is  8:1  which  is  sufficient  providing  the  next  highest  real  pole  is  suffi¬ 
ciently  small.  In  this  case,  the  next  highest  real  pole  is  located  at  s=  —9.  Using 
Equation  (3.14): 


^a-plane'J' 

2 


r-[ 


(-9)(0.02) 


=  0.0081  <  2 


(3.20) 
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which  meets  the  requirement  by  a  wide  margin.  Therefore,  even  though  some  warp¬ 
ing  occurs  due  to  the  pole  at  —50  there  is  little  warping  due  to  the  remaining  poles 
and  s  «  w'  is  assumed.  This  indicates  that  the  warping  in  going  from  the  w'-plane 
to  the  s-plane  should  be  minimal. 

3.3.1  The  Hofmann  Algorithm  In  this  design,  the  w -plane  representation 
of  the  plants  are  found  using  the  Hofmann  algorithm  [7].  The  Hofmann  algorithm 
allows  conversion  of  the  s-domain  representation  through  the  z-plane  and  into  the 
w -plane  without  the  normal  numerical  difficulties  [26].  An  important  feature  of 
Hofmann’s  algorithm  is  the  automatic  inclusion  of  the  zero  order  hold.  A  ZOH  is 
assumed  to  exist  between  the  controller  and  the  aircraft.  A  block  diagram  of  the 
sampled-data  system  is  included  (see  Figure  3.2).  The  Hofmann  algorithm  allows 


Figure  3.2.  Sampled  Data  System 
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conversion  of: 


and: 


to: 


and: 


sX(s)  =  AX(s)  +  BU(s) 

Y(s)  =  CX(s)  +  DU(s) 

wT’ 

w'X(w')  =  A*X(w')  +  1- V  B*U(w') 
Y(w')  =  CX(wO  +  DU(wO 


(3.21) 

(3.22) 

(3  23) 

(3.24) 


A  MATRIXjif  executable  file  is  used  to  perform  the  calculation  and  is  included 
below:^ 


//Hofmann. cmd 

t=l/60;  //sample  rate  of  Lambda 
Ca,b,c,d]=split(8,n3); 
ao=[eye-l/3*(a*t/2)**2t'2/16*(a*t/2)**4. . . 
-17/316*(an/2)>».*6... 

+62/2835* (a*t/2)**8. . . 

-1382/166926* (a*t/2) **10] ; 

aw=a*ae; 

be=ae*b; 

bH=be-t/2*aw*be ; 


'This  MATRIXx  executable  file  differs  from  the  one  included  in  Capt  Wheaton’s  Thesis  [26].. 
The  executable  file  included  in  Wheaton’s  thesis  contains  typographical  errors. 
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dw=d-t/2*c*be; 

SH= [aw.bHjcH.dw] ; 

Since  there  are  nineteen  plants  included  in  this  design,  the  Hofmann  algo¬ 
rithm  is  performed  nineteen  times.  This  is  acc.>mplished  by  using  the  programming 
capability  of  MATRIXjf  and  placing  the  above  code  in  an  iterative  loop.  The  follow¬ 
ing  is  a  w'-plane  representation  of  the  transfer  function  at  the  nominal  conditions: 

g(w')  Kv/'{v/'  +  0.005484)(w'  +  3.732)(w'  -  100)(w'  +  116.2)(w'  -  130.3)(w'  +  718.7) 

SeUv  ~  (w'  -f  0.001782  ±  i0.2299)(w'  +  3.225  ±  i8.022)(w'  -I-  9.010  ±  jf6.203)(w'  +  46.21) 

(3.28) 

{K  =  1.3414  X  10"^) 

The  above  equation  contains  non- minimum  phase  (NMP)  zeros  at  100  and  130.  They 
result  from  the  w'  transformation,  which  always  produces  a  transfer  function  with 
the  same  number  of  zeros  as  poles.  The  NMP  zero  at  100  is  due  to  the  sampling  rate. 
The  NMP  zero  at  130  is  due  to  the  excess  number  of  poles  over  zeros  in  the  s  domain 
transfer  function.  According  to  Horowitz  [10],  if  the  number  of  poles  exceeding  the 
number  of  zeros  is  two,  there  will  be  two  additional  zeros  created  and  one  will  be 
NMP.  A  complete  set  of  w -plane  transfer  functions  is  included  in  Appendix  B. 

3.4  Frequency  Response  Data 

The  QFT  technique  uses  templates  that  represent  the  plant  parameter  uncer¬ 
tainty,  at  the  specified  design  frequencies.  The  plant  uncertainty  for  this  thesis  is 
represented  by  the  nineteen  LTI  plants.  For  this  reason,  the  frequency  response  in¬ 
formation  for  all  of  the  plants  at  the  design  frequencies  is  necessary.  While  obtaining 
the  frequency  response  data,  a  problem  was  discovered  in  the  way  MATRIX;f  handles 
phase  information.  See  Appendix  E  for  a  discussion  of  the  problem.  The  frequency 
response  data  of  the  nineteen  longitudinal  plants  are  summarized  in  Figure  3.3. 
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Figure  3.3.  Frequency  Response  of  the  Nineteen  Longitudinal  Plants 
3.5  Plant  Templates 

Plant  templates  are  developed  from  the  frequency  response  data.  First  the 
template  frequencies  are  chosen  on  the  basis  of  plant  variation.  Where  the  plants 
vary  significantly,  the  template  frequencies  must  be  closely  spaced.  Where  the  plants 
are  closely  spaced,  the  frequencies  can  be  more  widely  spaced. 

For  this  design  the  frequencies  u=  0.001,  0.005,  0.01,  0.05,  0.1,  0.2,  0.3,  0.4, 
0.5,  0.6,  0.7,  1,  2,  3,  4,  5,  6,  8,  10,  20,  50,  100  (rad/sec)  are  used;  where  v  is  the 
w'  analog  of  u;  in  the  s-plane,  and  is  given  by  Equation  3.12. 
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Template  for  Frequency  (v)  =  5  [red/eee] 

Figure  3.4.  Example  of  Template  Used  in  Longitudinal  SISO  Design 

Next,  the  template  data  are  “stripped  out”  of  the  frequency  response  data. 
Twenty  two  sets  of  template  data  are  formed,  a  set  for  each  template  frequency. 
Each  set  contains  magnitude  and  phase  information  for  the  nineteen  plants  at  the 
template  frequency.  Only  the  magnitude  and  phase  differences  are  important  to  the 
QFT  technique.  In  plotting  the  templates,  the  data  is  normalized  by  subtracting  the 
smallest  phase  from  the  phase  data  and  the  smallest  magnitude  from  the  magnitude 
data.  The  templates  are  then  plotted  for  use  in  determining  the  stability  bounds. 
An  example  template  is  shown  in  Figure  3.4  for  u  =  5.  The  template  is  defined  by  a 
set  of  points,  each  representing  a  plant.  The  spread  of  the  points  is  due  to  the  plant 
variation  at  the  template  frequency.  In  the  example  template,  the  lowest  point  is  not 
located  at  the  origin  because  a  small  magnitude  and  phase  is  added  to  the  template 
data  to  distinguish  all  points  from  the  plot  axes.  The  circle  represents  the  nominal 
plant.  Choosing  the  nominal  plant  is  discussed  in  the  next  section.  The  data  used 
to  plot  the  twenty-two  templates  are  included  in  Appendix  D. 
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3.6  Nominal  Plant 


Any  one  of  the  plants  can  be  chosen  as  the  nominal  plant;  however,  choosing  the 
plant  with  the  least  number  of  right-half-plane  (RHP)  poles  and  zeros  may  make  the 
loop  shaping  easier  [10].  The  location  of  the  nominal  plant  on  each  of  the  templates 
is  needed  before  the  bounds — tracking,  disturbance,  and  stability — can  be  plotted  on 
the  Nichols  Chart  (NC).  Plotting  the  bounds  on  the  NC  is  key  to  the  QFT  technique. 
Because  the  plant  uncertainty  is  captured  in  the  templates  and  the  templates  are 
used  to  form  the  bounds,  the  single  transmission  loop  shape  meeting  or  exceeding 
the  plotted  NC  bounds  guarantees  each  plant  considered  will  meet  or  exceed  the 
tracking,  disturbance,  and  stability  performance  specifications.  The  nominal  plant 
chosen  for  the  longitudinal  SISO  design  is  flight  condition  one.  Additionally,  the  same 
flight  condition  is  used  in  the  lateral-directional  MIMO  system  design  discussed  in 
Chapter  IV. 

3.7  Stability  Bounds 

A  benefit  to  the  designer  using  a  predetermined  (fixed)  sampling  rate  is  the 
need  to  consider  only  the  stability  bounds  during  the  transmission  loop  shaping 
[10].  The  maximum  gain  in  the  loop  is  primarily  constrained  by  the  NMP  zeros. 
Therefore,  the  objective  is  to  “pack  in”  as  much  transmission  loop  gain  as  possible. 
If  the  tracking  and  performance  specifications  are  not  met,  a  change  in  the  sampling 
rate  must  be  made  to  allow  higher  loop  transmission  gain.  The  military  standard 
requirement,  from  MIL-STD-1797A  [16],  specifies  a  phase  margin  7  =  45°  to  be 
used  to  form  the  NC  stability  contour  Ml-  The  7  =  45°  requirement  corresponds 
to  a  3  dB  Ml  contour  on  the  NC.  Figure  3.5  is  an  example  of  a  NC  with  the  3  dB 
contour  highlighted. 

3.7.1  The  All-Pass-Filter  Technique.  The  all-pass-filter  (apf)  technique  is 
used  throughout  this  design  [12]  to  remove  the  NMP  elements  from  the  design  while 
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Phkte  [deg] 


Figure  3.5.  Nichols  Chart  With  3  dB  Contour  Highlighted 


still  accounting  for  their  elFects.  The  apf  technique  complicates  the  design;  however, 
it  allows  the  designer  to  rapidly  determine  the  maximum  possible  loop  transmission 
which  is  worth  the  additional  complication.  The  phase  of  the  apf  is  included 
when  plotting  the  stability  bounds.  The  apf  used  in  the  longitudinal  SISO  design  is: 

(w' -  100)(w' -  130) 

“  (w'  +  100)(w'  + 130) 

The  amount  of  phase  associated  with  the  template  frequencies  is  summarized  in 
Table  3.2.  At  low  frequencies,  below  v  =  0.5  {rad/sec)(u;  =  0.4999958[rad/sec]), 
the  shift  is  negligible.  At  higher  frequencies,  above  v  =  0.5  (rad/sec),  the  shift  is 
incorporated  in  the  stability  bounds.  The  moving  Nichols  chart  technique,  which 
is  slightly  different  from  the  method  usually  employed,  is  used  to  plot  the  stability 
bounds  on  the  NC. 
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Table  3.2.  All-Pass-Filter  Used  in  the  Longitudinal  SISO  System  Design 


Template  Frequency  v  (rad/sec) 

Phase  <^A  (deg) 

0.001 

0.002 

0.005 

0.01 

0.05 

0.10 

0.20 

0.40 

0.30 

0.61 

0.40 

0.81 

0.50 

1.01 

0.60 

1.21 

0.70 

1.41 

1.00 

2.02 

2.00 

4.04 

3.00 

6.06 

4.00 

8.08 

5.00 

10.1 

6.00 

12.1 

8.00 

16.1 

o 

o 

20.2 

20.0 

40.0 

50.0 

94.9 

100 

165 

3.1.2  The  Moving  Nichols  Chart  Technique.  Two  transparent  NCs  are  cre¬ 
ated.  The  3  dB  Ml,  contour  is  then  plotted  on  each  NC.  Instead  of  cutting  out 
the  plant  templates,  a  template  which  is  plotted  on  a  NC  with  the  same  scale  as 
the  transparent  NCs,  is  placed  beneath  one  of  the  transparent  NCs,  refered  to  as 
the  top  NC  (TNC),  (see  Figure  3.6).  In  this  way,  the  template  orientations  can  be 
maintained  by  aligning  the  grid  lines  of  the  TNC  with  the  grid  lines  on  the  template 
NC.  For  a  particular  location  of  the  template,  with  the  template  contour  touching 
the  stability  bound,  the  location  of  the  nominal  plant  point  is  marked  on  the  TNC 
(see  Figure  3.6).  The  template  is  then  moved  to  a  new  location  with  the  template 
contour  just  touching  the  Ml  contour  at  another  location.  Again  the  location  of 
the  nominal  plant  point  is  marked  on  the  TNC.  This  procedure  continues  until  a 
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With  3  dB  Ml  Contour 
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Figure  3.6.  Plotting  Stability  Bound  Using  Moving  Nichols  Chart  Technique 


stability  boundary  is  drawn  through  these  points  on  the  TNG.  Once  completed,  the 
bound  is  marked  with  the  appropriate  template  frequency.  This  procedure  is  fol¬ 
lowed  for  each  template  up  to  u  =  0.5.  For  frequencies  above  v  =  0.5,  the  second 
transparent  NC  is  used  and  is  refered  to  as  the  reference  Nichols  Chart,  RNC.  It 
is  used  to  indicate  the  phase  lead  necessary  to  account  for  the  lead  contributed 
by  the  NMP  elements  removed  from  the  design  by  the  apf.  The  phase  contribution 
must  be  incorporated  in  plotting  the  stability  bounds.  Before  the  stability  bounds 
arc  plotted,  the  Ml  contour  must  be  shifted  to  the  right  on  the  NC.  This  is  readily 
accomplished  by  placing  the  RNC,  with  the  Ml  contour  marked,  under  the  TNC. 
The  RNC  is  then  shifted  to  the  right  by  the  number  of  degrees  contributed  by  4>k 
associated  with  the  template  frequency  v,-.  Using  the  “shifted”  Ml  contour  on  the 
RNC,  stability  bound  plotting  can  continue.  Next,  the  template  is  placed  under  the 
RNC.  The  stability  bounds  are  obtained  by  moving  the  template  boundaries  about 
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Figure  3.7.  Plotting  Stability  Bound — Including  Phase  Contribution  of  apf 


the  shifted  contour  and  marking  the  nominal  plant  point  on  the  TNC  as  before  (see 
Figure  3.7).  The  resulting  stability  bound  is  plotted  through  these  points  on  the 
TNC  which  reflects  this  shift  to  the  right  of  the  stability  bound  by  (f>A,  the  number 
of  degrees  needed  to  account  for  the  NMP  character  of  the  plant  at  the  template 
frequency.  As  the  process  continues  for  each  template  frequency,  the  RNO  is  shifted 
further  to  the  right  by  the  appropriate  ^a-  This  continues  until  the  180°  line,  of  the 
RNC  is  on  the  0°  line  on  the  TNC  (see  Figure  3.8).  The  bottom  of  the  last  stability 
bound  dictates  the  amount  of  loop  transmission  gain  available  to  the  designer.  A 
NC  showing  shifted  stability  bounds  is  included  in  Figure  3.9. 

3.8  Loop  Shaping 

The  stability  bounds  are  used  in  shaping  the  loop  transmission.  The  design 
can  proceed  using  a  variety  of  techniques.  The  two  techniques  outlined  below  are 
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Figure  3.8.  Plotting  Stability  Bound — Plotting  the  Last  Bound 

the  “backward  loop  shaping  technique”  and  “minimum  order  compensator  design 
technique” .  A  high  order  compensator  is  not  desirable  for  Lambda  and  the  “back¬ 
ward  loop  shaping  technique”  may  produce  a  high  order  compensator,  however,  both 
loop  shaping  design  techniques  are  used  in  this  design.  Although  the  final  design  is 
based  on  the  “minimum  order  compensator  technique”,  the  “backward  loop  shaping 
technique”  is  used  to  determine  the  theoretical  maximum  loop  transmission.  In  the 
lateral-directional  MIMO  design,  covered  in  Chapter  IV,  the  “backward  loop  shaping 
technique”  is  used  as  a  performance  baseline, 

3.8.1  Backward  Loop  Shaping  Technique.  This  technique  takes  full  advantage 
of  the  minimum  phase  (MP)  character  of  the  plant  with  all  of  the  NMP  elements 
removed  by  the  apf.  In  MP  systems  the  phase  lag  is  correlated  with  the  slope  of  the 
gain  (dB/octave). 
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Figure  3.9.  Nichols  Chart  With  Stability  Bounds  Plotted 


This  characteristic  allows  the  designer  to  quickly  determine  the  maximum  loop 
transmission  gain  possible.  The  maximum  loop  transmission  gain  is  determined  by 
starting  at  the  frequency  corresponding  to  the  highest  stability  bound  and  working 
down  (backward)  in  frequency.  The  first  step  is  to  select  a  point  close  to  the  zero 
degree  line  and  slightly  below  the  highest  frequency  boundary  on  the  NC  (see  Fig¬ 
ure  3.9).  The  point  marked  “Ideal”  is  the  beginning  point  of  the  design.  Using  the 
fact  that  the  dB/octave  slope  is  determined  by  the  amount  of  phase  lag,  the  dB 
increase  to  the  next  lower  bound  is  determined  by  calculating  the  average  phase  lag 
between  the  bounds.  The  average  phase  lag  between  the  two  points  is  determined 
by  averaging  the  phase  lag  of  each  of  the  bounds: 


Phase  Lagav  = 


Phase  L^gBound  i  +  Phase  Lagg^^^nj  2 
2 


(3.30) 
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The  MP  character  of  the  design  allows  the  amount  of  dB  increase  to  be  determined 
based  on  the  average  phase  lag.  The  relationship  between  phase  and  the  gain  slope 
is: 

90°phaseLag  ~  6  dB/Octave  (3.31) 


This  relationship  is  used  to  determine  the  slope  at  the  average  phase  lag  of  the  two 
bounds.  For  example,  if  the  average  phase  lag  between  the  two  bounds  is  10°  the 
slope  is: 


(lOphaseLag) 


6  dB/Octave 


90phase  Lag 


=  0.667  dB/Octave 


(3.32) 


The  number  of  octaves  between  the  bounds  is  determined  from: 


Frequency  (Bound  l)  2(^Iuniber  of  Octaves) 

Frequency(Bound2) 


Or  more  directly  from: 


In 


FVequencyfB.und  n 
Frequency(Bound3) 


In  2 


=  Number  of  Octaves 


For  example,  the  number  of  octaves  between  v  =  50  and  u  =  40  is: 


ln(52) 

-r-^  =  0.332  Octaves 
In  2 


(3.33) 


(3.34) 


(3.35) 


The  maximum  dB  value  at  the  next  bound  is  the  product  of  the  number  of  octaves 
and  the  dB/octave  slope.  For  the  examples  above,  the  maximum  rise  in  dB  from  the 
bound  at  u  =  50,  to  the  bound  at  u  =  40,  is  approximately  0.22  dB.  This  procedure 
is  repeated  until  a  point  above  the  lowest  frequency  bound  is  reached  or  when  the 
designer  determines  the  design  is  not  possible.  If  the  design  is  not  possible,  the 
designer  chooses  a  lower  frequency  bound,  that  is,  for  this  example,  a  value  of  f  <  50 
(rad/sec)as  a  new  starting  point.  When  a  satisfactory  loop  shape  is  determined  by 
this  method,  it  results  in  the  highest  loop  transmission  gain  and  the  best  performance 


3-19 


possible.  The  maximum  loop  transmission  provides  the  most  benefit  from  feedback 
[10].  However,  the  loop  transmission  formed  using  this  procedure  may  produce  a 
higher  order  compensator  than  desired.  If  the  compensator  is  limited  to  a  lower 
order,  it  can  be  reduced  by  absorbing  some  of  the  nominal  plant  poles  and  zeros 
into  the  loop  transmission  transfer  function.  If  the  order  is  still  too  high,  the  entire 
plant  can  be  included  for  the  final  transmission  loop  design  as  described  below.  By 
using  the  above  procedure,  the  loop  shaping  information  obtained  by  excluding  the 
nominal  plant  poles  and  zeros  is  useful  in  determining  the  theoretical  maximum  loop 
transmission  gain  that  is  achievable.  In  addition,  the  higher  order  design  may  be 
used  as  a  performance  baseline  for  comparison  with  lower  order  designs. 

3.8.2  Minimum  Order  Compensator  Design  Technique.  The  previous  tech¬ 
nique  is  most  useful  if  the  loop  is  formed  from  scratch.  That  is,  the  nominal  plant  is 
not  included  in  the  design  of  the  transmission  loop.  The  drawback  is  the  possibility 
of  an  unacceptably  high  order  compensator.  The  order  can  be  reduced  somewhat 
by  absorbing  as  many  poles  and  zeros  of  the  nominal  plant  as  possible  into  the 
loop  transmission  equation.  Therefore,  a  satisfactory  lower-order  compensator  is 
obtained  by  including  the  nominal  plant  in  the  design.  The  nominal  plant  is  first 
multiplied  by  the  inverse  of  the  all-pass-filter,  (apf)“^ ,  which  produces  a  MP  transfer 
function.  Thus,  the  resulting  MP  transfer  function,  designated  by  P„,o,  is  plotted  on 
the  NC  or  on  a  set  of  Bode  plots.  Loop  shaping  begins  at  the  lowest  frequency  and 
progresses  until  a  point  below  the  highest  frequency  bound  is  reached.  The  result  of 
the  backward  loop  shaping  technique  is  used  as  the  design  goal.  Loop  shaping  by 
this  technique  produces  the  following  MP  transmission  loop  shape  and  is  used  for 
the  remainder  of  the  longitudinal  SISO  design: 

L  =  0-0183(w^  +  2)(w^  +  5)(w^  -H  10)(w^  -h  50)P,„o 

(w02(w'  +  160±;T20)(w'-i-300)  ^  ’  '' 
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Figure  3.10  is  the  MP  loop  transmission  Lmo  plotted  on  the  NC.  A  NC  showing  a 
plot  of  all  the  plants  including  the  NMP  elements  is  shown  in  Figure  3.11. 


PhMe  [4eg] 

(■4’l00X«'»^130Mlw2)rmo  Niebol*  Chart 

Figure  3.10.  NC  Showing  MP  Transmission  Loop  After  Loop  Shaping 


8.9  w'  Compensator 

The  w'  compensator,  G(w'),  is  found  by  dividing  the  loop  transmission  by  the 
nominal  plant.  If  the  plant  is  included  in  the  loop  transmission,  G(w')  is  simply 
the  poles,  zeros,  and  gain  added  during  loop  shaping.  If  the  nominal  plant  is  not 
included  in  the  transmission  loop,  the  zero  and  poles  not  absorbed  in  L„,o  become 
poles  and  zeros,  respectively,  of  G(w').  The  longitudinal  SISO  G(w')  for  this  design 
is: 


,  ^  5.5(w72  +  l)(w75  -H  l)(w78  -f  I)(w7l0  +  l)(w750  +  1) 

^  ’  (w')2(w7200H1.6w7200  +  I)(w7300  +  1)  ^ 
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Figure  3.11.  NC  Including  All  Nineteen  Plants  and  NMP  Elements 


or  equivalently; 


G(w') 


0.0183(w'  +  2)(w'  +  5)(w'  +  8)(w'  +  10)(w'  +  50) 
(w')2(w'  +  160  ± il20)(w'  +  300) 


(3.38) 


3.10  Prefilter  Design 

The  prefilter,  F(w'),  is  obtained  by  plotting  the  frequency  response  data  of  the 
closed  loop  system:’ 

L(w')  \  /  G(w')P(w')  \ 

1  +  L(w')y  “  Vl  +  G(wOP(w')  j 

The  frequency  response  bounds,  developed  from  the  specifications,  are  plotted  over 

the  closed  loop  response.  A  family  of  curves  is  produced.  See  the  first  plot  in 

Figure  3.12.  Poles  and  zeros  are  placed  in  front  of  the  closed  loop  system  until  all 

of  the  responses  are  contained  within  the  performance  bounds,  as  in  the  second  plot 

in  Figure  3.12.  The  prefilter  for  the  longitudinal  SISO  QFT  design  is: 
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Figure  3.12.  Closed  Loop  Response  Without  and  With  Prefilter,  F(w') 


,  _  0.001875(w^  +  200)^ 
^  ’  (w'  +  5)(w'  +  15) 


(3.40) 


where: 


lim  =  1 

w'-*0 


(3.41) 


S.ll  w'  Simulation 

The  w'-plane  design  is  simulated  to  verify  the  performance  before  proceeding 
with  the  design.  The  w'  simulations  are  performed  using  the  System  Build  capability 
of  MATRIXx  (see  Figure  3.13).  Time  domain  step  responses  are  plotted  along  with 
the  performance  bounds  (see  Figure  3.14)  which  indicate  that  a  robust  design ’s 
achieved,  that  is,  all  responses  lie  between  the  upper  and  lower  specified  i  racking 
bounds. 
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Figure  3.13.  System  Build  w'  Simulation 


Figure  3.14.  System  Build  w'  Simulation  Frequency  Response 
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3.12  z-Plane  Transformations 

F(w')  and  G(w')  are  transformed  to  the  z-plane  for  digital  implementation 
on  Lambda.  The  z-plane  representations  are  obtained  by  performing  the  bilinear 
transformation;  see  Equation  (3.6).  The  MATRIX^  user  function  below  readily 
performs  this  transformation: 


//  [numz ,  denz]  =r2z  (tium ,  den ,  dt ) 

//Function  to  take  a  -C\rm  «}’  domain  numerator  and  denominator  into 
//the  z-domain. 
sz=size(num) ; 
n=sz(2)-l; 

//multiply  the  (2/dt)  factors  into  the  num  and  den 
for  i=l:n; . . . 

num(i)=num(i)*(2/dt)**((N+l)-i) ; . . . 
den(i)=den(i)*(2/dt)**((H+l)-i) ; . . . 
end; 

//build  the  transformation  matrix 
tm=ones ( (n+1) , (n+1) ) ; 
facta=Cl  -1] ;factb=Cl] : 
for  j=l:(n+l);... 
roB=l; . . . 

for  i=j :n;row=conv(roR,f acta) ; end; . . . 
for  i=l: (j-1) ;row=conv(row,factb) ;end; . . . 
for  i=l: (n+1) ;tm(j ,i)=roB(i) ;end; . . . 
end; 
tm, 

//multiply  the  num  and  den  vectors 

//with  the  matrix 

numz=numi‘tm; 

denz=den*tm; 

//normalize  the  transfer  function 
numz=numz/denz ( 1 ) ; denz=denz/ denz ( 1 ) ; 


In  addition  to  performing  the  transformation,  the  above  function  divides  the 
numerator  and  denominator  by  the  denominator  leading  coefficient  resulting  in  a 
normalized  z-plane  representation: 


F(z) 


0.0139752^  +  0.00931682  +  0.0015528 
22  -  1.64392  -f  0.66874 


(3.42) 


G(z)  =  -  374.312^^  +  566.962^  -  417.222^  +  147.362  -  19.538 

2^  -  0.768292“  -  0.8780523  +  0.1707322  -  0.365852  +  0.10976  ’ 
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3.12.1  Comparison  of  vf'  and  z-Plane  F’s  and  G’s  .  The  required  differ¬ 
ence  equations  are  obtained  from  the  z-plane  representations  of  the  design  and  are 
included  in  Appendix  F,  The  w'  and  z-plane  compensators  and  filters  must  be 
comparable  up  to  frequencies  greater  than  (rad/sec)  [11].  The  sampling  fre¬ 
quency  on  Lambda  is  50  Hertz,  which  means  =  314  (rad/sec).  Therefore,  the 
w'  and  Z-plane  representations  must  be  fairly  close  up  to  a  w  =  105  (rad/sec). 
The  frequency  response  comparison  of  G(w')  and  G(z)  are  plotted  in  Figure  3.16 
for  comparison.  In  addition,  the  frequency  response  of  F(w')  and  F(z)  are  plotted 
in  Figure  3.15.  The  plots  show  a  good  comparison  within  the  required  frequency 
bandwidth  0  <  w  <  105  (rad/sec). 


Figure  3.15.  Comparison  of  F(w')  and  F(z)  Magnitude  Frequency  Response 


3.13  Simulations 

Simulation  is  the  final  step  in  the  design.  The  design  is  simulated  using  the 
best  model  available.  In  this  case,  the  best  model  available  is  the  same  one  used  in 
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the  design  with  the  additional  non-linear  elements  added.  The  non-linear  elements 
are  due  to  the  software  and  hardware  limiters  present  on  Lambda  (see  Figure  3.17). 
Figure  3.18  is  the  MATRIX^  System  Build  representation  of  Lambda  without  the 
non-linear  elements  while,  Figure  3.20  includes  the  non-linear  elements.  MATRIX^ 
has  the  ability  to  simulate  a  system  with  both  continuous  and  discrete  elements. 
This  “hybrid”  simulation  capability  is  used  to  verify  the  performance  of  the  design. 
The  design  is  simulated  without  the  non-linear  elements.  The  results  are  plotted 
in  Figure  3.19.  Adding  the  limiters  produces  the  simulation  results  shown  in  Fig¬ 
ure  3.21.  The  simulations  predict  the  performance  expected  from  Lambda  in  actual 
flight.  They  show  good  performance  at  all  flight  conditions  for  a  step  command  of 
10°  per  second  for  at  least  two  seconds.  After  two  seconds,  the  deflection  limits  on 
the  actuators  begin  to  affect  the  response. 
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Figure  3.17.  System  Build  Model  of  Lambda  Including  Non-Linear  Elements 
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Figure  3.18.  Model  Used  for  Hybrid  Simulation  of  Longitudinal  SISO  QFT  Design 
3.14  Summary 

A  single  QFT  controller  was  synthesized  that  produces  satisfactory  perfor¬ 
mance  throughout  the  flight  envelope  of  Lambda.  It  is  important  to  note  that  no 
gain  scheduling  is  necessary.  Response  of  the  system  remains  stable  in  spite  of  the 
inherent  non-linearities  present  due  to  hardware  and  software  control  limiters.  The 
hybrid  simulation  with  non-linear  elements  (see  Figure  3.21)  is  the  final  test  of  the 
design  as  it  includes  all  that  is  known  about  the  plant.  Some  of  the  performance 
limitations  can  be  overcome  if  Lambda’s  speed  is  held  constant.  In  the  simulations, 
Lambda  can  hold  a  pitch  rate  constant  until  the  speed  drops.  As  the  speed  drops, 
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Figure  3.19.  Hybrid  Simulation  of  the  Longitudinal  SISO  QFT  Design 


more  elevator  is  needed  than  is  available  to  the  controller.  Since  the  simulation  re¬ 
sults  are  good,  within  the  performance  limitations  of  Lambda,  it  is  assumed  that  the 
aircraft  will  perform  well  with  this  controller  design. 
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IV.  Lateral-Directional  MEMO  Design 

This  chapter  covers  the  synthesis  of  the  QFT  compensators  and  prefilters  for 
the  two-by-two  lateral-directional  MIMO  plants,  and  follows  the  steps  outlined  in 
Chapter  1  Section  1.7.  The  MIMO  QFT  design  technique  differs  only  slightly  from 
the  SISO  QFT  technique,  and  only  the  differences  are  covered  in  detail.  Method 
Two  requires  sequential  design  of  the  yaw  and  roll  loops;  however,  they  are  presented 
simultaneously. 

4-1  Flight  Conditions 

The  same  flight  conditions  used  in  the  longitudinal  SISO  design  are  used  in 
this  lateral-directional  MIMO  QFT  design. 

4.2  Plant  Transfer  Function  Matrix 

The  plant  transfer  function  matrix,  P(w'),  is  obtained  from  a  MATRIX^  sim¬ 
ulation  of  the  aircraft  at  the  selected  flight  conditions.  The  A  and  B  matrices  are 
of  the  form: 
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and: 


B/at  = 


ibiibJ 


ibJibJ 


L^r+fb]N^.  N,,„,+[b]L,, 


ibJibJ 


ibJibJ 


The  transfer  function  matrix  resulting  from  the  MATRIX^  simulation,  which 
includes  the  actuator  and  sensor  dynamics,  is  a  two-by-two  matrix: 


P_  a  ^  M^) 


where  paen  is  the  roll  rate  sensed  by  the  roll  gyro  and  Taen  is  the  yaw  rate  sensed  by 
the  yaw  gyro.  and  Srudcmd  commanded  aileron  and  rudder  deflections, 

respectively.  At  the  nominal  conditions,  the  elements  of  P(s)  are: 

_ 3.5051  X  10^(s)(3  +  0.9516  ±j3.8448) _  . 

{s  +  0.0141)(s  +  0.9734  ±  j3.9l81)(s  +  6.7492)(5  +  9.0000  ±  ;6.2450)(5  +  50.0000)  ^  ^ 

= 

_ -72.7664(g)(s  H-  3.9245)(s  -  83.4616) _  , 

{$  +  0.0141)(5  +  0.9734  ±  i3.9181)(s  +  6.1996)(s  +  6.7492)(5  +  50.0000)  ^  '  ' 

_ -9.5335  X  103(g-  1.8300)(s  +  2.3301)(s  + 22.9473) _ 

(s  +  0.0141)(5  +  0.9734  ±  i3.9181)(s  +  6.7492)(s  +  9.0000  ±  ;6.2450)(s  +  50.0000)  ^  ' 


-  SrIZjg)  - 


_ 5.4760  X  10^(8  +  0.0658  db  j0.3201)(s  +  6.7753) _ 

(8  +  0.0141)(s  +  0.9734  ±  j3.9181)(8  +  6.1996)(s  +  6.7492)(s  +  50.0000) 
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A  complete  set  of  transfer  functions  is  located  in  Appendix  B. 


4.3  w'  Plant  Transfer  Functions 

The  s-plane  transfer  functions  of  Equation  4.3  are  transformed  to  the  w'  do¬ 
main  as  described  in  Section  3.3  using  the  Hofmann  routine  which  takes  into  account 
the  ZOH  yield  Pe(w'). 

4.4  MIMO  Plant  Inversion,  P"^(w') 

The  MIMO  QFT  technique  requires  an  mxm  MIMO  plant  be  reduced  to 
equivalent  multiple-input  single-output  (MISO)  loops  in  order  to  perform  a 
QFT  design  [12].  The  lateral-directional  dynamics  of  Lambda  are  represented  by 
a  2x2  MIMO  plant  which  is  converted  to  four  equivalent  loops.  Since,  a  diagonal 
prefilter  F  is  specified,  the  loops  on  the  diagonal  have  two  inputs — a  desired  input 
and  a  disturbance  input  and  are  MISO  loops.  The  off-diagonal  elements  have  only 
disturbance  inputs  and  are  represented  as  SISO  loops  (see  Figure  4.1).  The  desired 
input  is  the  control  input  from  the  compensator  g,,-,  while  the  disturbance  input  is 
due  to  the  cross  coupling  in  the  MIMO  plant.  This  design  uses  a  compensator  matrix 
G  and  the  prefilter  matrix  F,  which  are  diagonal,  that  is: 


gn 

0 

0 

G(u;')  = 

0 

g22 

0 

(4,8) 

0 

0 

g33 

and: 

fii 

0 

0 

F(«;0  = 

0 

f22 

0 

(4.9) 

0 

0 

fss 

Therefore,  only  the  main  diagonal  MISO  loops,  elements  922  and  733,  have  desired 
control  inputs.  Disturbance  inputs  exist  in  each  equivalent  loop. 
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Figure  4.1.  Equivalent  Loops 


The  transfer  functions  used  in  the  equivalent  loops  are  obtained  by  first  invert¬ 
ing  the  plant  matrix  Pe(w')  and  then  inverting  each  element  of  Pe~’(w').  Symboli¬ 
cally: 


PliK) 

pW®')  Pk(®') 


(4.10) 


and: 


Q(®')  = 


1 

1 

9n 

9l2 

phK) 

PlaK) 

= 

1 

Pjjfw') 

1 

Paafw')  J 

921 

922 

(4.11) 


The  elements  of  Q(w’)  are  the  transfer  functions  of  the  equivalent  loops  (see  Fig¬ 
ure  4.1).  Once  the  q,-,s  are  obtained,  the  design  procedure  is  the  same  as  the  SISO 
design  procedure  with  the  plant  transfer  function  replaced  with  q,-,-  (see  Figure  4.2). 


4.4-1  The  Improved  QFT  Method.  The  Improved  Method  [9]  is  used  to  re¬ 
duce  the  uncertainty  of  the  MISO  disturbance  inputs.  The  design  starts  with  the 
loop  with  the  least  uncertainty  due  to  the  cross  terms.  The  resulting  loop  transmis¬ 
sion  is  used  in  the  subsequent  loops.  In  this  case,  the  second  loop  is  designed  with 
a  modified  transfer  function,  q22e- 

Because  the  yaw  channel  has  the  lowest  bandwidth  ,  the  design  of  the  MISO 
loop  cissociated  with  element  (3,3)  of  Q  is  completed  first,  followed  by  the  MISO  loop 
associated  with  element  (2,2).  The  following  equation  is  used  to  find  the  effective 
transfer  function,  q22e: 


922e  =  922 


(l  +  Ts) 

1+  L3  —  732 


(4.12) 


where: 

933922 

732  =  - 

932923 


(4.13) 


Figure  4.2.  w'-domain  QFT  Design  for  Fxjuivalent  MISO  Loops 


4-5 


4.5  Frequency  Response  Data 


The  frequency  response  data  needed  to  obtain  the  plant  templates,  is  obtained 
for  all  nineteen  plants.  Figures  4.3  and  4.4  are  the  Bode  plots  of  q33(w')  and  q22e(w') 
respectively,  from  which  the  frequency  response  data  are  obtained. 


Figure  4.3.  Frequency  Response  of  the  Nineteen  Equivalent  MISO  Plants,  q33, 


Figure  4.4.  Frequency  Response  of  the  Nineteen  Equivalent  MISO  Plants,  q22e, 
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4-6  Plant  Templates 

The  plant  template  data  for  each  template  are  taken  from  the  frequency  re¬ 
sponse  by  stripping  out  the  magnitude  and  phase  of  each  plant  at  the  template 
frequency.  The  plant  templates  are  then  plotted  in  obtain  the  required  stability 
bounds.  A  sample  of  a  plant  template  is  included  in  the  previous  chapter. 

4.7  Nominal  Plants 

Since  all  q,i  are  all  open  loop  stable,  the  same  conditions  used  to  create  the 
nominal  plant  for  the  longitudinal  SISO  design  are  used  to  generate  the  nominal 
plants  used  in  the  lateral-directional  MIMO  system  design.  The  nominal  conditions 
are  the  first  conditions  listed  in  Table  3.1.  The  nominal  equivalent  MISO  loop 
transfer  functions  used  in  the  lateral-directional  MIMO  design  are: 
q33{w')= 

3.1345  X  lO-'^Cw'  +  0.005)(w'  +  0.1523)(w'  -  100)(w'  -t-  115.9)(ii;  -  130.9)(w'  -t-  678.4)  . 

(w' -f  0.06488  ±i0.3224)(w'-i-6.7661)(w' -I- 9.0107  ±i6.2026)(w'-f  46.37)  ^  ' 

and: 

q22e(w')  = 

1.6304  X  10-3(w'  +  0.1519)(w'  -  100)(w'  +  149.05)(w'  -  208.14) 

(w'-f  0.9748  ±j3.9179)(w'-l-6.7389)(w' -1-46.212)  ^  ’  '' 

A  complete  set  of  equivalent  MISO  loop  transfer  functions  used  in  the  design  are 
included  in  Appendix  C. 

The  NMP  elements  of  the  equivalent  plants  are  removed  to  facilitate  the  ap¬ 
plication  of  the  MP  design  technique,  as  in  the  longitudinal  SISO  design.  The  apf’s 
used  in  the  lateral-directional  MIMO  design  are: 
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(w' +  100)(w' +  208) 
33  -  (y  _  ioo)(w'  -  208) 

and: 

(w'  +  100)(w'  +  131) 
22  -  (y  _  ioO)(w'  -  131) 


(4.16) 


(4.17) 


4-8  stability  Bounds 

Stability  bounds  are  obtained  using  the  moving  Nichols  chart  technique  dis¬ 
cussed  in  Section  3.7.2. 

4.9  Loop  Shaping 

Loop  shaping  is  performed  using  both  techniques  discussed  in  Section  3.8.  The 
backward  loop  shaping  technique  generates  a  maximum  loop  gain  within  the  limits 
caused  by  the  NMP  character  of  the  sampled  data  plant.  Once  the  loop  transmission 
is  obtained,  an  attempt  to  reduce  the  order  of  the  compensator  is  made  by  absorbing 
poles  and  zeros  of  the  nominal  plant  into  the  loop.  Failing  to  reduce  the  compensator 
sufficiently,  the  loop  shape  is  reformed  using  the  second  loop  shaping  technique:  the 
minimum  order  compensator  design  technique.  Since  the  owners  of  Lambda,  and 
sponsors  of  this  thesis,  desire  low-order  compensators,  design  by  both  methods  is 
necessary.  Only  the  designs  resulting  from  the  minimum  order  compensator  design 
technique  are  presented. 

Figures  4.5  and  4.6  are  the  MP  transmission  loops  for  qaa  and  q22e,  respectively, 
using  these  minimum  order  compensators.  The  NMP  transmission  loops  for  all 
nineteen  flight  conditions  are  included  in  F'igures  4.7  and  4.8.  From  the  figures, 
it  is  apparent  that  the  fifth  order  compensators  are  sufficient  to  maintain  stability. 
However,  the  higher  order  compensators  allow  more  loop  transmission  gain  and  result 
in  better  tracking  and  disturbance  rejection. 
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30 


Phase  [deg] 

Figure  4.6.  NC  of  MP  Transmission  Loop  Shape  for  q22e 
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4-10  w'  Compensator 

The  compensators  are  extracted  from  the  loop  transmissions  by  dividing  out 
the  nominal  plants.  The  resulting  w'  compensators  are: 


g(w')33  = 


188.89(w^  +  1.5)(w^  +  5)(w'  +  15)(w^  +  20)(w^  +  33) 
(w')(w'  +  i.7)(w'  +  60)(w'  +  100)(w'  +  208) 


(4.18) 


and: 


g(w')22e  = 


86.097(w'  +  0.5)(w'  +  1.4  ±  il.4283)(w'  +  7  ±  il2.124) 
(w')3(w'  +  75±jl29.9) 


(4.19) 


4,11  Prefilter  Design 

Synthesis  of  the  prefilters  follows  the  procedure  discussed  in  Section  3.10.  The 
closed-loop  magnitude  frequency  responses  of  L33  and  L22e  are  plotted  along  with 
the  frequency  response  tracking  bounds.  Poles  and  zeros  are  chosen  for  fj,-(w')  so 
that  all  of  the  responses  of: 


f(w0 


g(w0g.f(w0 

1  +  giyf')qii{^') 


are  contained  within  the  bounds  (see  Figures  4.9  and  4.10). 


(4.20) 
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4- IS  z-Plane  Transformations 

The  next  step  is  to  transform  the  compensators  and  prefilters  into  the  z-plane. 

Again,  a  MATRIXjf  user  function  is  used  for  the  transformation  (see  page  3- 
25).  The  resulting  z-plane  representations  are: 


0.014908^2  -  0.014884^ 
(2)33  -  ^2  _  5  9861^  ^  0.98610 


(4.21) 


and: 


g(^)33 


37.484^®  -  142.552'*  +  214.132®  -  158.6622  +  57.9432  -  8.3350 
2®  -  1.88322'*  +  0.7193523  +  0.2444322  -  0.0805472 


(4.22) 


and: 


1(2)22 


0.0849672  -  0.045752 
2  -  0.96078 


(4.23) 


and  finally: 


21.7242®  -  100.472'*  +  I86.6O2®  -  +  81.6822  -  15.426 

2®  -  2.47372^  +  1.789523  -  0.5263222  -  0.578952  -  0.36842  ^  ^ 


4.IS.I  Comparison  of  W  and  z-Plane  f’s  and  g’s  .  The  w'  and  z-plane 
representations  of  the  lateral-directional  compensators  and  prefilters  are  compared. 
As  in  the  longitudinal  case,  they  must  be  comparable  up  to  —  105  (rad/sec)(see 
Figures  4,11,  4.12,  4.13,  and  4.14),  As  in  the  design  for  the  longitudinal  channel, 
there  are  differences  in  the  w'  and  z-plane  representations.  However,  they  are  in 
good  agreement  up  to  the  required  frequency. 


4-13 


Figure  4.13.  Comparison  of  f22(w')  and  £22(2)  Magnitude  Frequency  Response 


Prequancy  [rad/acc] 

Figure  4.14.  Comparison  of  g22(w')  and  g22(z)  Magnitude  Frequency  Response 
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4- IS  Simulations 

The  system  was  simulated  on  MATRIXjt  and  the  performance  was  verified 
by  comparing  the  step  response  to  the  performance  bounds.  Once  the  performance 
was  verified,  the  nonlinear  limiter  elements  were  added  to  the  simulation.  Again,  the 
performance  was  verified.  The  results  of  the  simulations  are  contained  in  Figures  4.15 
and  4.16,  and  the  results  of  the  simulations  containing  the  nonlinear  elements  are 
contained  in  Figures  4.17  and  4.18.  In  Figures  4.15  and  4.16,  the  units  are  radians  per 
second  for  the  roll  and  yaw  rates  and  radians  for  the  aileron  and  rudder  deflections. 
A  unit  step  input  was  used  for  the  simulations.  In  Figures  4.17  and  4.18,  the  units 
for  the  roll  and  yaw  rates  are  degrees  per  second,  while  the  units  of  the  aileron 
and  rudder  deflection  are  degrees.  A  step  of  10°  was  used  in  the  simulations  which 
included  the  nonlinear  elements.  In  the  linear  system,  the  size  of  the  input  does 
not  effect  the  response  of  the  system.  For  convenience,  a  unit  step  is  used.  After 
the  nonlinear  elements  are  added  to  the  simulations,  however,  the  level  of  the  input 
directly  affects  the  system  response.  For  example.  Lambda  commands  are  limited 
by  hardware  and  software  limiters  to  a  15°  maximum  aileron  deflection.  Thus,  an 
aileron  command  of  20°  will  have  the  same  effect  as  a  15°  command. 

4.14  Summary 

With  the  exception  of  generating  the  MISO  equivalent  transfer  functions,  the 
QFT  design  for  the  lateral-directional  MIMO  portion  of  Lambda  follows  the  same 
procedure  as  the  longitudinal  SISO  channel.  The  design  is  effective  but  somewhat 
limited  by  the  low-order  compensator  requirement.  Again,  adding  the  nonlinear 
elements  in  the  model  limited  the  response  but  did  not  produce  instabilities.  Better 
performance  could  be  obtained  through  faster  sampling,  higher  bandwidth  actuators, 
and  higher-order  compensation,  all  of  which  would  allow  higher  loop  transmission 
gain.  In  spite  of  the  physical  and  imposed  limitations,  the  QFT  controller  performs 
well. 
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Figure  4.17.  Simulation  of  Lateral-Directional  MIMO  System  With  Step  Roll  Rate 
Command  Including  Nonlinear  Elements 
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Figure  4.18.  Simulation  of  Lateral-Directional  MIMO  System  With  Step  Yaw  Rate 
Command  Including  Nonlinear  Elements 


V.  Autopilot  Design 


This  chapter  presents  an  autopilot  design  for  Lambda  based  on  the  QFT  rate 
controller  discussed  in  earlier  chapters.  Included  are  the  design  requirements,  syn¬ 
thesis,  and  simulations. 

5.1  Requirements 

The  Control  Systems  Development  Branch  of  the  Wright  Laboratory  requires 
a  flight  control  system  to  aid  the  Lambda  test  pilot.  The  QFT  rate  controller  meets 
part  of  the  requirement  by  stabilizing  the  aircr^dt  and  guaranteeing  performance 
throughout  the  flight  envelope  by  reducing  the  effective  plant  variation.  Thus  tests 
can  be  conducted  without  redesigning  or,  at  least,  installing  a  new  controller  for  each 
test  flight.  But,  autopilot  functions,  such  as  turn  coordination  and  attitude  control, 
are  also  desired. 

This  design  extends  Capt  Wheaton's  [26]  wing  leveler  to  include  turn  coordi¬ 
nation  and  roll  attitude  control.  It  is  an  outer  loop  for  the  QFT  rate  controller  using 
specifications  from  MIL-C-18244A  [17].  The  standard  requires  the  roll  attitude 
remain  within  0.5®  of  the  commanded  roll.  Smooth  and  rapid  return  to  commanded 
attitude  after  a  disturbance  is  also  specified.  A  disturbance  of  20®  will  return  to  a 
commanded  attitude  with  one  overshoot  or  less  with  a  peak  overshoot  less  than  4®  . 

5.2  Synthesis 

The  autopilot  is  built  as  an  outer  loop  aiound  the  QFT  rate  controller,  (see 
Figure  5.1).  It  is  designed  directly  on  MATRIX^  using  System  Build.  For  the 
purpose  of  simulation,  the  '-rale  space  model  of  Lambda  is  modified  to  generate 
airspeed  variations  u  as  an  output.  The  simulation  requires  total  airspeed  which  is 
made  up  of  the  sum  of  the  nominal  airspeed  U  with  the  airspeed  variations  u.  Ciong 
is  changed  and  the  state  space  is  reconstructed.  On  the  aircraft,  the  velocity  sensor 
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Figure  5.1.  Block  diagram  of  Autopilot  Design  Using  QFT  Controlled  Lambda 


provides  total  velocity  17+  u.  Roll  attitude  ^  is  obtained  by  integrating  the  roll 
rate  output  p.  Additional  integrators  generate  simulation  outputs  for  pitch  and  yaw 
attitude  9  and  respectively.  On  the  aircraft  the  attitude  is  sensed  with  positional 
gyros.  A  disturbance  input  is  added  along  with  variable  control  inputs.  Since  the 
design  must  work  throughout  the  flight  envelope,  the  nominal  speed  of  the  aircraft 
is  applied  as  an  input  to  the  system.  The  Lambda  test  pilot  desires  the  ability  to 
command  bank  angle  [18].  Based  on  previous  test  flights,  a  wing  leveler  on  the 
autopilot  is  also  desirable.  Rate  control  of  pitch  and  yaw  is  desirable,  particularly 
during  landing.  Landing  is  normally  performed  using  pitch  rate,  making  it  the  more 
natural  choice.  Yaw  can  be  controlled  with  rate  or  attitude.  Since  there  is  no 
advantage  to  either  methods  and  rate  control  of  yaw  will  be  simpler  to  implement 
because  the  turn  coordination  scheme  in  Reference  [3]  provides  a  yaw  rate  in  response 
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Figure  5.2.  Lambda  Autopilot  Using  QFT  Compensated  Aircraft 

to  a  commanded  roll,  the  autopilot  design  maintains  rate  control  of  pitch  and  yaw 
but  allows  the  pilot  to  directly  command  <j>. 

The  QFT  compensated  Lambda  does  not  require  any  additional  compensation 
to  enable  the  aircraft  to  respond  to  the  pilot’s  command  of  (f>.  A  simple  unity- 
feedback  loop  is  added  connecting  the  ^  output  to  the  Additional  circuitry 

is  needed  for  pitch  and  yaw  turn  coordination  (see  Figure  5.2).  The  pitch  and  yaw 
coordination  equations  are  obtained  from  Reference  [3]  and  are  repeated  below. 

The  pitch  coordination  equations  are: 


Qcoord 


— tan  ^  sin  ^ 
(/  +  u 


(5.1) 


and: 

record  =  frr~  sin  <f>  (5.2) 

where:  q  is  the  pitch  rate  (rad/sec),  g  is  the  acceleration  due  to  gravity  (32.1  ft/sec^), 
U  is  the  nominal  speed  (ft/sec),  and  r  is  the  yaw  rate  (rad/sec).  <f)cmd  enters  the 
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controller  in  degrees  and  is  converted  to  radians.  The  equations  are  simplified  for 
implementation  on  Lambda  using  a  small  angle  approximation: 


and: 


^coord  — 


C^  +  u 


<i>‘‘ 


(5.3) 


^coord  — 


U  +  n 


4' 


(5.4) 


The  (j)  input  to  the  coordination  circuits  is  taken  from  the  pilot’s  command  to 
prevent  the  coordination  circuits  from  responding  to  bank  angle  disturbances  caused 
by  outside  influences  such  as  gusts.  The  nominal  speed  is  entered  from  outside  the 
controller,  while  the  changes  in  speed  u  are  fed  back  (for  simulation).  The  model 
used  in  the  autopilot  simulations  is  the  QFT-compensated  Lambda  that  includes 
the  nonlinear  limiters  (see  Figure  3.17). 


5.3  Simulations 

Simulation  was  performed  on  MATRIX^  •  Each  plant  condition  was  included 
in  the  simulation  through  program  iteration.  The  first  simulation  represented  in 
Figure  5.3  is  the  response  of  the  autopilot  to  a  step  disturbance  of  20°  .  The  au¬ 
topilot  returns  the  aircraft  to  the  commanded  roll  angle  of  0°  within  approximately 
2  seconds,  with  little  or  no  overshoot.  A  small  error  in  yaw  angle  is  created  by  the 
maneuver,  which  must  be  taken  out  by  the  rudder  and  ailerons.  The  error  remains 
constant  until  the  rudder  actuator  is  saturated  after  about  four  seconds.  Interven¬ 
tion  by  the  pilot  is  required.  This  characteristic  is  due  to  the  limit  on  the  control 
surface  and  the  lower-order  compensator  design. 

The  second  simulation  represented  in  Figure  5.4  shows  the  response  of  the 
autopilot  to  a  -10°  roll  command  with  a  -1-10°  roll  disturbance.  Again,  the  aircraft 
model  responds  within  specifications.  In  this  simulation,  it  is  evident  that  the  system 
is  more  than  a  wings  lever  as  it  returns  the  aircraft  to  the  commanded  .oil  angle  of 
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-10°  .  As  before,  the  pilot  must  intervene  to  keep  from  saturating  the  controls  after 
about  four  seconds.  In  this  case,  the  airspeed  decreases  as  expected,  because  the  turn 
coordination  computer  is  commanding  additional  pitch  rate.  The  additional  pitch  is 
required  to  keep  the  flight  path  horizontal  throughout  the  turn.  As  Lambda  slows, 
more  elevator  is  needed  to  generate  the  required  lift.  Thus,  the  airspeed  decreases 
steadily.  An  automatic  throttle  control  may  prevent  this  situation  by  holding  the 
airspeed  constant. 

The  third  simulation  (see  Figure  5.5)  shows  the  autopilot  response  to  a  con¬ 
stant  30°  roll  command.  In  this  case,  the  autopilot  very  accurately  responds  to 
the  command;  the  roll  angle  is  established  in  about  2.5  seconds  and  held  constant. 
Again,  control  saturation  of  the  elevator  and  rudder  actuators  are  reached.  In  this 
simulation,  the  rudder  saturates  for  some  conditions  as  early  as  4.6  seconds,  while 
for  the  high  speed  conditions  it  does  not  saturate  within  10  seconds.  The  elevator 
actuator  saturates  at  6.5  seconds  for  a  few  of  the  conditions.  Again,  an  automatic 
throttle  control  may  extend  the  capability  of  holding  a  desired  roll  angle. 

Figures  5.6  and  Figures  5.7  show  the  autopilot  response  to  a  constant  10°  pneh 
rate  command  and  a  2  second  10°  pulse  pitch  command,  respectively.  The  response 
is  the  same  as  the  QFT  rate  controller  without  the  autopilot  outer  loop,  which  is  to 
be  expected.  The  2  second  10°  pulse  pitch  command  is  more  realistic  as  it  results  in 
a  20°  pitch  angle. 

The  10°  constant  yaw  rate  command  simulation  in  Figure  5.8  and  the  2  second 
10°  pulse  yaw  rate  simulation  in  Figure  5.9  have  the  same  response  as  the  rate 
controller  without  the  autopilot. 

The  next  simulation.  Figure  5.10,  shows  the  autopilot  response  to  a  command 
consisting  of  a  4  second  !-10°  roll  angle  command  followed  bv  a  4  second  -10°  roll 
angle  command.  The  model  follows  the  commands  well  and  doesn’t  saturate  the 
control  actuators. 
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Figure  5.3.  Autopilot  Response  to  a  20®  Roll  Disturbance  Input 
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Figure  5.4.  Autopilot  Response  to  a  -10®  Roll  Command  With  a  +10®  Roll  Distur¬ 
bance 
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The  last  three  simulations  involve  a  ramp  roll  position  input. 

Figure  5.11  is  the  autopilot  response  to  a  very  slow  (F/sec)  ramp.  The  rudder 
actuator  begins  to  saturate  for  some  of  the  plants  after  18  seconds.  Only  four  of  the 
plants  experience  the  saturation  within  the  20  second  simulation. 

Figure  5.12  response  is  due  to  a  ramp  with  a  peak  of  20/dg  and  a  slope  of 
5°/sec.  The  aircraft  tracks  the  ramp  well  with  a  small  (0.5  sec)  delay.  At  5  seconds, 
all  of  the  nineteen  plants  achieve  the  desired  roll  angle  and  hold  it  throughout  the 
simulation  well  within  ±  0.5°  specification  of  MIL-STD-18244A.  Four  of  the  plants 
saturate  the  elevator  after  15  seconds.  Eight  plants  exhibit  rudder  saturation  starting 
after  nine  seconds. 

The  final  simulation,  Figure  5.13,  demonstrates  the  autopilot  response  to  a 
triangular  shaped  roll  input.  The  peak  input  of  20  °  is  realized  at  four  seconds  after 
the  start  of  the  simulation.  At  the  peak,  the  slope  is  reversed  until  the  input  is  zero 
at  eight  seconds.  In  this  case,  no  actuators  saturate. 
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Figure  5.5.  Autopilot  Response  to  a  Constant  30®  Roll  Command 
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Figure  5.7.  Autopilot  Response  to  a  2  Second  10®  Pulse  Pitch  Command 


Figure  5.8.  Autopilot  Response  to  a  Constant  10®  Yaw  Rate  Command 
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Figure  5.11.  Autopilot  Response  to  Ramp  Roll  Input  With  a  20“  Peak  and  Slope 
of  l“/sec 
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Figure  5.12.  Autopilot  Response  to  a  Ramp  Roll  Input  With  a  20“  Peak  and  a 
Slope  of  5“/sec 
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5.4  Summary 

A  satisfactory  autopilot  design  was  easily  synthesized  around  the  QFT  con¬ 
troller  due  to  the  reduced  uncertainty,  stability,  and  predictable  performance  pro¬ 
vided  by  the  controller.  The  autopilot  remained  stable  and  responsive  to  commands 
until  control  saturation  was  reached.  After  one  or  more  of  the  controls  became 
saturated,  the  model  did  not  go  unstable.  An  autopilot  can  be  a  great  aid  to  the 
pilot  but,  in  this  design,  a  pilot  is  still  necessary.  The  better  performance  could 
be  expected  from  a  higher-order  controller  with  a  higher  loop  transmission  gain. 
In  addition,  this  autopilot  uses  one  of  four  turn  coordination  schemes  presented  in 
Reference  [3].  Another  scheme  may  perform  better  or,  possibly,  the  sine  and  tangent 
functions  could  be  used  instead  of  a  small  angle  approximation.  Of  course,  a  higher 
loop  transmission  gain  would  be  obtained  with  an  increased  sample  rate  or  increased 
order  compensator.  The  main  benefit  of  the  higher  gain  would  be  reduced  variation 
in  the  responses  particularly  in  the  yaw  channel. 
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VL  Summary,  Conclusions,  and  Recommendations 

This  chapter  presents  a  summary  of  the  thesis,  draws  conclusions,  and  recom¬ 
mends  future  actions. 

6.1  Summary 

The  QFT  design  technique  was  applied  to  the  Lambda  URV  and  autopilot  func¬ 
tions  were  added  as  an  outer  loop  to  the  QFT  controller.  In  contrast  with  previous 
efforts  that  used  models  developed  from  computer  predictions  from  design  specifi¬ 
cations,  this  thesis  employs  a  model  developed  from  flight  test  data.  In  the  model, 
aircraft  dynamics  of  the  Lambda  model  separate  into  two  independent  channels — the 
longitudinal  channel  and  the  lateral-directional  channel.  The  longitudinal  channel  is 
a  SISO  system  while  the  lateral-directional  channel  is  a  two-by-two  MIMO  system. 

The  longitudinal  SISO  QFT  design  follows  the  steps  outlined  in  Chapter  I 
Section  1.7  and  proceeds  from  the  choice  of  flight  conditions  through  simulations 
using  the  SISO  QFT  technique  in  the  w'domain.  Flight  conditions  were  chosen  with 
great  care  to  insure  a  realistic  usable  design  with  minimal  over-design.  Stability 
bounds  were  plotted  using  graphical  techniques,  while  performance  bounds  are  not 
used  due  to  the  reasons  explained  in  Chapter  III  Section  3.7.  The  loop  shaping 
is  done  with  the  entire  nominal  plant  included  in  the  loop  as  part  of  a  minimum 
realization  technique.  The  compensator  is  extracted  by  dividing  the  plant  out  of  the 
resulting  transmission  loop.  After  the  prefilter  is  formed,  the  design  is  transformed 
to  the  z-plane  and  simulated  using  the  hybrid  simulation  capability  of  MATRIX^  • 

The  lateral-directional  MIMO  QFT  design  proceeds  in  a  similar  manner.  Only 
the  differences  in  the  design  techniques  are  discussed  in  detail.  The  flight  conditions 
used  in  the  longitudinal  SISO  design  are  the  same  ones  used  in  the  longitudinal 
SISO  QFT  design,  including  the  choice  of  the  nominal  plant.  An  equivalent  transfer 
function  mai  rix  is  developed  and  inverted  to  create  the  Q  matrix  with  elements  qu 
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which  are  the  equivalent  MISO  loop  transfer  functions.  Once  the  qu  are  obtained, 
the  design  follows  the  SISO  QFT  design  steps.  The  transmission  loop  is  formed 
using  two  separate  techniques — a  backward  loop  shaping  technique,  Section  (3.8.1), 
and  a  minimum  order  compensator  design  technique.  The  backward  loop  shaping 
technique  is  applied  first  to  establish  a  theoretical  maximum  performance  baseline 
for  the  design.  Then,  application  of  the  minimum  order  compensator  design  tech¬ 
nique  produces  a  fifth  order  compensator  required  for  implementation  on  Lambda’s 
digital  flight  control  computer.  After  the  first  loop  of  the  MIMO  system  is  designed, 
the  subsequent  loops  are  designed  using  modified  qa's.  Comparing  the  results  of 
simulations  with  specifications  completes  the  design. 

An  autopilot  is  added  cis  an  outer  loop  built  around  the  QFT  compensated 
Lambda.  It  allows  roll  attitude  control  along  with  providing  turn  coordination, 
mode  requested  by  the  Wright  Laboratory.  The  autopilot  modes  are  designed  to 
meet  military  specifications  and  the  simulations  prove  the  design.  For  example,  the 
wings  return  to  level  in  about  2  seconds  after  a  20°  disturbance  with  little  or  no 
overshoot  and  meets  MIL-STD-18244A.  The  QFT  controller  effectively  minimized 
the  parameter  variations  in  the  aircraft.  Therefore,  the  design  is  simple,  and  only 
one  set  of  controllers  is  necessary.  The  turn  coordination  modes  work  well  but  are 
limited  by  the  performance  of  the  aircraft.  Control  authority  is  exceeded  in  many 
cases,  and  though  the  design  is  stable  throughout  the  envelope,  pilot  intervention  is 
necessary. 

6.2  Conclusions 

The  QFT  design  technique  is  powerful  and  perfect  for  the  design  of  a  rate  con¬ 
trolled  autopilot  for  URVs  requiring  robust  autopilot  systems.  Lambda  simulations 
show  a  single  QFT  controller  can  perform  well  throughout  the  entire  flight  envelope 
of  Lambda  without  additional  sensors.  The  technique  is  straightforward.  Synthesis 
and  simulations  can  be  performed  using  a  computer  aided  design  package  such  as 
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MATRIX^  using  minimal  computer  resources  compared  with  the  hours  of  computer 
time  devoured  by  fellow  students  using  other  control  techniques.  The  method  is  com¬ 
pletely  transparent,  in  that  the  designer  is  able  to  predict  the  ultimate  performance 
early  in  the  design.  More  importantly,  the  designer  is  able  to  make  engineering 
tradeoffs  during  loop  shaping.  For  example,  performance  at  certain  frequencies  can 
be  increased  by  reducing  the  performance  at  frequencies  of  lower  interest  to  the 
designer. 

6.S  Recommendations 

Working  on  a  controller  for  Lambda  is  a  unique  opportunity.  Simply  knowing 
that  the  controller  design  will  be  used  is  highly  moti^'ating.  The  Wright  Laboratory 
should  continue  to  challenge  AFIT  students  with  thesis  proposals  involving  Lambda. 
The  students  benefit  by  facing  a  real  life  challenge  and  the  laboratury  gets  needed 
research.  The  effect  is  synergistic.  This  design  should  be  implemented  on  Lambda 
and  the  performance  should  be  compared  with  predictions.  The  Lambda  model 
should  be  improved  as  more  information  becomes  available  and  should  be  extended 
to  include  the  effects  of  operating  the  ten  flight  control  surfaces  independently.  The 
weight  of  Lambda  should  be  reduced  or  ,at  least,  held  to  its  current  value  since  the 
heavier  weights  are  the  main  source  of  pl^lnt  variation  and  contain  open  loop  unstable 
poles.  An  automatic  throttle  control  should  also  be  added.  The  Wright  Laboratory 
should  acquire  the  ability  to  run  MATRIX^  in  order  to  use  the  computer  software 
generated  during  this  thesis  to  upgrade  the  controller  if  it  becomes  necessary.  Ad¬ 
ditionally,  other  QFT  controllers  should  be  synthesized  for  Lambda.  This  design 
is  aimed  at  stability  throughout  the  envelope  assuming  no  failures.  Future  efforts 
should  consider  the  effects  of  component  failure  or  battle  damage  on  the  controlla¬ 
bility  of  the  aircraft.  Because  of  limited  computer  resources  aboard  Lambda,  this 
controller  is  limited  to  fifth  order  compensation.  Superior  performance  was  achieved 
in  simulation  with  higher  order  compensation  due  the  higher  available  loop  trans- 
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mission  gain.  Future  efforts  might  add  the  computer  resources,  without  additional 
weight,  to  Lambda  to  allow  higher  order  compensation  and  a  faster  sampling  rate. 
Finally,  Lambda’s  performance  could  be  improved  by  installing  faster  actuators  since 
the  present  actuators  severely  limit  the  loop  transmission  gain. 
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Appendix  A.  Model  Generating  MATLAB  macro  File 

The  following  MAT  LAB  “m”  file  was  created  by  Lt  Gerald  Swift  to  generate 
a  state  space  representation  of  Lambda  [24]; 


y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

% 

y. 

y. 

y. 

y. 

y. 

% 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 

y. 


This  macro  calculates  the  longitudinad  and  lateral- 
directional  state  space  equations  of  motion  for  the 
Lambda  URV  given  Xcg,  Ul,  ql,  W  and  theta. 

Input:  X= [Xcg, Ul,ql,W, theta] 

where,  Xcg  =  inches 

Ul  =  feet /second 
ql  =  pounds/sq  foot 
H  =  pounds 
theta  =  degrees 


Outputs: 

Clong= [CMu , CHa , CMad , CMq , CLu , CLa , CLad , CLq , CDa ,  CDu , 

CLde,CDde,CMde,CLdf ,CDdf ,CMdf] ' 

Clat=[Clb,Clp,Clr, Clda , Cldr , Cnb , Cnp , Cnr , Cnda , Cndx , 

Cyb , Cyp , Cyr , Cyda , Cydr] ’ 

Dlong  =  corresponding  vector  of  dimensional 

derivatives  for  Clong 

Dlat  =  corresponding  vector  of  dimensional  derivatives  for  Clat 

Along  ==  A  matrix  for  longitudined  equations  of  motion 

Alat  =  A  matrix  for  lateral -direct ionzd  equations  of  motion 

Blong  =  B  matrix  for  longitudined  equations  of  motion 

Blat  =  B  matrix  for  lat eral -direct ioned  equations  of  motion 


CMu=0; 
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CMa=0 . 326643*X(1 , 1) -16 . 5189 ; 

CMad=0 . 238235*X ( 1 , 1 ) "7 . 8977 ; 

CMq=0 . 628420+X (1 , 1) -45 . 0582 ; 

CLu=-0.00017l9»X(l,l)+0. 01518; 

CLa=5 . 820 ; 

CLad=-0 . 03838*X(1 , l)+3 . 83508; 

CLq=-0 . 65202+X ( 1 , 1) +37 . 8884 ; 

CDa=0 ; 

CDu=0; 

CLde=0.2908; 

CDde=0 ; 

CMde=0 .013464*X(1 , 1)-1 .48039; 

CLdf =1.419; 

CDdf =0.08489; 

CMdf =0 . 070932*X(1 , 1) -3 . 63838 ; 

Clongl=  [CMu , CMa ,  CMad , CMq ,  CLu , CLa, CLad ,  CLq , CDa ,  CDu]  ; 

Clong2= [CLde .CDde ,CMde , CLdf ,CDdi , CMdf]  ; 

Clong= [Clongl  Clong2] ' ; 

X 

Clb=-0. 01451 ; 

Clp=-.5538; 

Clr=0. 08763; 

Clda«0.2608; 

Cldr=0 . 000213*X(1 , 1) -0 . 00783 ; 

Cnb=-0 . 00038*X ( 1 , 1 ) +0 . 07834 ; 

Cnp=-0. 03601; 

Cnr=0 . 003087*X(1 , 1) -0 . 31072 ; 

Cnda=-0. 01368; 

Cndr»0. 00 1633*X( 1,1) -0.15208;  %  updated  4/22/91 

Cyb=-0.4372; 

Cyp=-0.001600; 

Cyr»-0 . 00424+X ( 1 , 1) +0 . 46047 ; 

Cyda=0 ; 

Cydr=0 . 2865 ; 

Clat * [Clb , Clp , Clr , Clda , Cldr , Cub , Cup , Cnr , Cnda , Cndr , 

X  Cyb , Cyp , Cyr , Cyda , Cydr] ' 

X 

X  Calculation  of  Dimensional  Derivatives 

X 

X  Dlong=  [Mu ,  Ma ,  Mad ,  Mq ,  Zu ,  Za ,  Zad ,  Zq ,  Xa ,  Xu ,  Zde , 

X  Xde,Mde,Zdf,Xdf,Mdf]' 

X 
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y,  Dlat=[Lb,Lp,Lr,Lda,Ldr,Nb,Np,Nr,Nda,N<ir,Yb,Yp,Yr,Yda,Ydr] ' 

y. 

xcg=X(l,l); 

u=X(l,2); 

q=X(l,3); 

m=X(l,4)/32.174; 

s=21.1; 

c=1.51; 

b=14.07; 

Ixx=^2.502; 

Iyy=26.666; 

Izz=40 . 939 ; 

Ixz=0.5922; 

CLl=X(l,4)/(q*s) ; 

CD1=0.027; 

CM1«0.0; 

Al=Ixz/Ixx; 

Bl=Ixz/Izz; 

X 

Mu®q*s*c*CCMu+2*CMl)/(Iyytu) ; 

Ma»q*s*c*CMa/Iyy ; 

Mad«q*s»c*c*CMad/(2*Iyy*u) ; 

Mq=q*s*c*c*CMq/(2*Iyy*u) ; 

Zu=“q*s*(CLu+2*CLl)/(ni*u) ; 

Za=“q*s* (CLa+CDl) /m ; 

Zad=-q*s*c*CLad/(2*a*u) ; 

Zq=-q*s*c*CLq/(2*m»u) ; 

Xu=-q*s*(CDu+2*CDl)/(in*u) ; 

Xa*-q*s* (CDa-CLl ) /m ; 

Zde=-q*  s  *CLde/m ; 

Xde=-q*s*CDde/m ; 

Mde=q* s*c*CMde/Iyy ; 

Zdf =-q*  s*CLdf /m ; 

Xdf =-q*s*CDdf /m ; 

Mdf =q*s»c*CMdf /lyy ; 

Dlong=  [Mu  ,Ma ,  Mad ,  Mq ,  Zu ,  Za ,  Zad ,  Zq ,  Xa ,  Xu ,  Zde ,  Xde ,  Mde ,  Zdf ,  Xdf ,  Mdf  ]  » 

y. 

Lb=q*s*b*Clb/Ixx ; 

Lp=q*s*b*b*Clp/(2*Ixx*u) ; 

Lr=q*s*b*b*Clr/(2*Ixx*u) ; 

Lda=q*s*b*Clda/Ixx ; 

Ldr»q*s*b*Cldr/Ixx ; 
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Nb»q*s*b*Cnb/Iz2 ; 

Np=q*s*b*b*Cnp/ (2*Izz*u) ; 

Nr=q*s*b*b»Cnr/(2*Izz*u) ; 

Nda=q*s*b*Cnda/Izz ; 

Ndr=q*s*b*Cndr/Izz ; 

Yb=q*s*Cyb/in; 

Yp=q*s*b*Cyp/(2*m*u) ; 

Yr=q*s*b*Cyr/(2*m*u) ; 

Yda=q*s*Cyda/m ; 

Ydr®q*s*Cydr/in ; 

Dlat=[Lb,Lp,Lr,Lda,Ldr,Nb,Np,Nr,Nda,Kdr,Yb,Yp,Yr,Yda,Ydr]  ’ ; 

% 

%  Calculation  of  the  Longitudinal  A  and  B  matrices 

% 

thet a*X( 1 , 5) /57 . 29578 ; 

Along*[Xu,Xa,0,“32. 174*cos(theta) ; 

Zu/(u-Zad) ,Za/(u-Zad) , (u+Zq)/ (u-Zad) ,-32 . 174*sin (theta) ; 
Mu+Mad*Zu/  (u-Zad)  ,Ma+Mad*Za/  (u-Zad)  ,Mq+Mad*(u+Zq)/(u-Zad) , . 

-Mad*32. 174*sin (theta) ; 

0,0, 1,0]; 

Blong=CXde,Xdf; 

Zde/ (u-Zad) ,Zdf/ (u-Zad) ; 

Mde+Had*Zde/ (u-Zad)  ,Mdf+Mad*Zdf/ (u-Zad) ; 

0,0] ; 

% 

%  Calculation  of  the  Lateral  Directional  A  and  B  matrices 

% 

BA«1-B1*A1; 

Alat = [Yb/u , Yp/u , 32 . 174*cos (theta) / u , Yr/u- 1,0; 

(Lb+Al*Nb)/BA, (Lp+Al*Np)/BA,0, (Lr+Al*Nr)/BA,0; 

0,1, 0,0,0; 

(Nb+Bl+Lb) /BA , (Np+Bl^Lp) /BA , 0 , (Mr+Bl*Lr) /BA , 0 ; 

0,0,0,1,01; 

Blat= [Yda/u , Ydr/u ; 

(Lda+Al*Nda)/BA, (Ldr+Al*Ndr)/BA; 

0,0; 

(Nda+Bl*Lda)/BA, (Ndr+Bl*Ldr)/BA; 

0,0]; 
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Appendix  B.  liransfer  Functions 


This  appendix  includes  all  transfer  functions  used  in  the  rate  controller  design. 


B.l  Longitudinal  Transfer  Functions 

The  flight  conditions  produce  nineteen  separate  plants.  The  s-domain  and  w'- 
domain  transfer  functions  used  in  the  longitudinal  SISO  design  are  included  below: 
Nominal  s-Domain  Transfer  E\inction  (v-^W  > ): 

Poit  pill 

K{s){s  +  0.0548)(s  +  3.7340) 

(«  +  0.0178  ±  i0.2298)(«  +  3.2050  ±  i8.0129)(«  +  9  ±  i6.2450)(«  +  50) 

{K  =  1.5774  X  10*) 

Off-Nominal  s-Domain  Transfer  F\inctions 

P2ii  — 

_ A'(g)(g  +  0.0519)(5  +  3.1423) _ 

(5  +  0.0150  ±  i0.2323)(«  +  2.8788  ±  jSM2Z){s  +  9  ±  i6.2450)(«  +  50) 

(A' =  1.5766  X  10*) 


P3ii  ~ 


K(s)(s  +  0.0583)(«  +  4.2689) 

(5  +  0.0206  ±  i0.2277)(»  +  3.6639  ±  i8.5618)(«  +  9  ±  j6.2450)(5  +  50) 


(K  =  1.8045  X  10*) 


P*Xl  = 


K{s}{s  +  0.0540)(s  +  3.5946) 

(5  +  0.0173  ±  i0.2304)(s  +  3.2913  ±  j8.5978){«  +  9  ±  j6.2450)(»  +  50) 


{K  =  1.8035  X  10*) 


Ps„  = 


Kis){s  +  0.0542)(s  +  3.87.54) 

{s  +  0.0173  ±  j0.2201)(»  +  3.0638  ±  i6.2369)(s  +  9  ±  j6.2450){«  +  .50) 


{K  =  1.5310  X  10*) 
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Pen  = 


I<{s){s  +  0.0511)(s  +  3.2619) 

(s  +  0.0145  ±  i0.2238)(s  +  2.7374  ±  i6.2664)(s  +  9  ±  i6.2450)(s  +  50) 
{K  =  1.5301  X  10®) 


P7n  = 


K{s)(s  +  0.0577){s  +  4.4299) 

(s  +  0.0202  ±  i0.2169)(s  +  3.5023  ±  i6.6723)(s  +  9  ±  j6.2450)(s  +  50) 


(a:  =  1.7515  X  10®) 


P8n  = 


K{s){s  +  0.0534)(s  +  3.7308) 

(s  +  0.0169  ±  ;0.2210)(s  +  3.1295  ±  i6.7084)(s  +  9  ±  ;6.2450)(s  +  50) 


{K  =  1.7504  X  10®) 


P9u  = 


_ /t:(s)(a+0.1141)(s-f  2.1412) _ 

(s  +  0.0033  ±  i0.4028)(s  +  1.9054  ±  j4.3264)(s  +  9  ±  i6.2450)(fi  +  50) 


{K  =  5.5923  X  10^) 


PlOii  — 


K{s)is  +  0.0729)(s  +  2.7725) 

(s  +  0.0102  ±  i0.3199)(s  +  2.3980  ±  i5.4539)(s  +  9  ±  i6.2450)(5  +  50) 


(K  =  8.9051  X  10“) 


Piiii  = 


K{s){s  +  0.0596)(s  +  3.4480) 

(s  +  0.0152  i  j0.2647)(s  +  2.9539  ±  j6.6i)30)(s  +  9  ±  j6.2450)(s  +  50) 


(K  =  1.3261  X  10®) 


Pnn  = 


ii:(s)(s  +  0.2872)(s  +  1.3270) 

(s  +  0.0178  ±  j0.5275)(fi  +  1.2835  ±  i2.5894)(5  +  9  ±  ;6.2450)(s  +  50) 


{K  =  2.5836  X  10“) 


Pl3„  = 


/t:(5)(5  +  0.3825)(s  +  0.9785) 

(s  +  0.0166  ±  i0.5372)(«  +  1.1470  ±  ;2.5972)(s  +  9  ±  i6.2450)(s  +  50) 


(K  =  2.5822  X  10“) 


Pl4„  = 


_ /C(g)(5  +  0.2427)(5  +  1.6011) _ 

(s  +  0.0147  ±  j0.5207)(s  +  1.4620  ±  j2.7686)(s  +  9  ±  j6.2450)(s  +  50) 


{K  =  2.9564  X  10“) 
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Nominal  w -Domain  Transfer  Function 


Pan  —  Pill  — 

/-i:(w')(w'  +  0.0548)(w'  +  3.7323)(w'  -  100)(w'  +  116.1585)('v'  -  130.3470)(w'  +  718.6759) 
(w'  +  0.0178  ±  i0.2298)(w'  +  3.2246  ±i8.0218)(w'  +  9.0106  ±  j6.2026)(w'  +  46.2117) 


(/!:  =  1.3414  X  10-^) 


Off-Nominal  w'-Domain  Transfer  Functions 


P2ii  = 


A'(w')(w'  +  0.0519)(w'  +  3.1413)(w'  -  100)(w'  +  116.1627)(w'  -  130..S37S)(w'  +  719.3433) 
(w'  +  0.0150  db  j0.2323)(w'  +  2.8967  ±i8.0529)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 

{K  =  1.3397  X  lO-'*) 


P3ii  = 


A(w')(w'  +  0.0583)(w'  +  4.2663)(w'  -  100)(w'  +  116.1373)(w'  -  130.4125)(w'  +  714.8263) 
(w'  +  0.0206  ±  j0.2277)(w'  +  3.6892  ±  j8.5712)(w'  +  9.0106  ±  j6.2026)(w'  +  46.2117) 

(K  -  1.5432  X  10-*) 


Pill  = 


A:(w')(w'  +  0.0540)(w'  +  3.5931)(w'  -  100)(w'  +  116.1422)(w'  -  130.4019)(w'  +  715.5924) 
(w'  +  0.0173  ±  i0.2304)(w'  +  3.3145  ±  j8.6096)(w'  +  9.0106  ±  ;6.2026)(w'  +  46.2117) 

(K  -  1.5408  X  10"^) 

Poll  - 


K(w')('v'  4. 0.0542)(w'  +  3.8735)(w'  -  100)(w'  +  1 16.1758)(w'  -  130.2617)(w'  +  722.8078) 
(w'  H-  0.0173  ±  ;0.2201)(w'  +  3.0748  ±  j6.2391)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 

{K  =  1.2931  X  10-4) 


P6„  = 


/C(w0(w'  +  0.0511)(w'  +  3.2607)(w'  -  100)(w'  +  116.l785)(w'  -  130.2555)(w'  +  723.2503) 
(w'  +  0.0145  ±  i0.2238)(w'  +  2.7475  ±  i6.2699)(w'  +  9.0106  db  i6.2026)(w'  +  46.2117) 

{K  =  1.2917  X  10-4) 
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P7ix  = 


K{-w>){w'  +  0.0577)(w'  +  4.4270)(w'  -  100)(w'  +  116.1572)(w'  -  130.3152)(w'  +  719.5097) 
(w'  +  0.0202  ±  j0.2169)(w'  +  3.5165  ±  i6.6740)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 


(K  =  1.4863  X  10-^) 


PSu  = 


K{vt'){w'  +  0.0534)(w'  +  3.7291)(w'  -  100)(w'  +  116.1603)(w'  -  130.3080)(w'  +  720.0226) 
(w'  +  0.0169  ±  j0.2210)(w'  +  3.1426  ±  i6.7119)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 


(K  =  1.4844  X  lO--*) 


P9„  = 


+  0.1141)(w'  +  2.1408)(w'  -  100)(w'  +  116.2107)(w'  -  130.1436)(w'  +  729.6175) 
(w'  +  0.0033  ±  i0.4028)(w'  +  1.9088  ±  i4.3275)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 

(K  =  4.6762  X  10-®) 


PlOu  = 


+  0.0729)(w'  +  2.7718)(w'  -  100)(w'  +  116.1888)(w'  -  130.2142)(w'  +  725.4450) 
(w'  +  0.0102  ±  )0,3199)(w'  +  2.4047  ±  i5.4562)(w'  +  9.0106  ±  ;6.2026)(w'  +  46.2117) 

(K  =  7.4917  X  10-®) 


Pu„  = 


^!:(w')(w'  +  0.0596)(w'  +  3.4467)(w'  -  100)(w'  +  116.1650)(w'  -  130.2956)(w'  +  720.8195) 
(w'  +  0.0152  ±  j0.2647)(w'  +  2.9662  ±  j6.6570)(w'  +  9.0106  ±  j6.2026)(w'  +  46.2117) 

(K  =  1.1233  X  10-“*) 


Pl2u  = 


K{vi'){v!'  +  0.2872)(w'  +  1.3269)(w'  -  100)(w'  +  116.2475)(w'  -  130.0370)(w'  +  736.3425) 
(w'  +  0.0178  ±  j0.5275)(w'  +  1.2842  ±  i2.5895)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 

(K  =  2.1401  X  10-®) 
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Pl3u  = 


Kivf'){W  +  0.3825)(w'  +  0.9784)(w'  -  100)(w'  +  116.2486)(w'  -  130.0345)(w'  +  736.5274) 
(w'  +  0.0166  ±  i0.5372)(w'  +  1.1478  ±  i2.5974)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 


(/i:  =  2.1385  X  10-®) 


Pl4ji  = 


+  0.2427)(w'  +  1.6010)(w'  -  100)(w'  +  116.2393)(w'  -  130.0578)(w'  +  734.9239) 
(w'  +  0.0147  ±  j0.5207)(w'  +  1.4630  ±  i2.7687)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 

{K  =  2.4535  X  10-5) 


Pl5„  = 


/i:(w')(w'  +  0.3032)(w'+  1.2516)(w'-100)(w'  +  116.2406)(w'-  130.0548)(w'+ 735.1389) 
(w'  +  0.0141  ±i0.5311)(w'  +  1.3077  ±i2.7794)(w'  +  9.0106  ±i6.2026)(w'  +  46.2117) 

{K  =  2.4513  X  10-5) 


Pi6„  = 


Kiv/'){v/'  +  0.3037)(w'  +  1.2527)(w'  -  100)(w'  +  116.2380)(w'  -  130.0666)(w'  +  734.5254) 
(w'  +  0.0114  ±  i0.5544)(w'  +  1.3368  ±  i3.3072)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 

(/r  =  2.2107  X  10-5) 


Pl7ii  = 


+  0.4180)(w'  +  0.8941)(w'  -  100)(w'  +  1 16.2397)(w'  -  1 30.0628)(w'  +  734.8052) 
(w'  +  0.0104  ±  j0.5608)(w'  +  1.2004  ±  i3.3168)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 

{K  =  2.2089  X  10-5) 


Pl8,i  = 


K(y;’){v/'  +  0.2544)(w'  +  1.5234)(w'  -  100)(w'  +  116.2283)(w'  -  130.0016)(w'  +  732.8558) 
(w'  +  0.0093  ±  i0.5495)(w'  +  1 .5250  ±  i3.5334)(w'  +  9.0106  ±  i6.2026)(w'  +  46.21 17) 

{K  =  2.5352  X  10-5) 
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K(w')(w'  +  0.3217)(w'  +  1.1774)(w'  -  100)(w'  +  116.2303)(w'  -  130.0871)(w'  +  733.1781) 
(w'  +  0.0087  .-t  i0.5565)(w'  +  1.3697  ±  j3.5457)(w'  +  9.0106  ±  j6.2026)(w'  +  46.2117) 

(K  =  2.5328  X  10-®) 
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B.2  Lateral-Directional  Transfer  Functions 


The  laberal-directional  dynamics  are  represented  by  a  two-by-two  transfer 
function  matrix: 

p  _  P22  P23 

P32  P33 

The  chosen  flight  conditions  produce  nineteen  separate  transfer  function  matrices. 
The  s-domain  and  w'-domain  transfer  functions  are  included  below: 

Nominal  s-Domain  Transfer  Function  po22  =  {saa}j]n(s)) 


(B.l) 


P02U  —  PI33 


/:(«)(«  + 0.9516  ±j3.8448) 

(s  +  0.0141)(s  +  0.9734  ±  i3.9181)(s  +  6.7492)(s  +  9  ±  j6.2450)(s  +  50) 


{K  =  3.5051  X  10®) 


Off-Nominal  s-Domain  Transfer  Functions  =  (g 

P2a3  — 

K(s)(s  +  0.9271  ±  j3.8429) 

(s  +  0.0142)(s  +  0.9490  ±  i3.9158>;5  +  6.7487)(c  +  9  ±  i6.2450)(s  +  50) 
{K  =  3.5051  X  10®) 


P333  — 


K(s)(s  +  1.0882  ±i4.1010) 

(s  +  0.0141)(s  +  1.1179  ±  i4.1789)(«  +  7.7060)(5  +  9  ±  j6.2450)(s  +  50) 


{K  =  4.0081  X  10®) 


Piaa  — 


X(s)(s  +  1.06Cl.-fcj4.0987) 

(5  +  0,0141  )(5  +  1.0900  ±j4.1761)(s  +  7.7056)(s  +  9  ±  i6.2450)(s  +  50) 


(K  =  4.0081  X  10®) 


PS33 


A:(s)(5  +  0,9230  ±  i3.8259) 

(s  +  0.0144)(s  +  0.9448  ±  j3.8997)(s  +  6.7492)(5  +  9  db  ,76.2450)(.s  +  50) 


(A'  =  3.5051  X  10®) 


PS33  — 


/<■(«)(«  + 0.8984  ±j3.8241) 

(s  +  0.0144)(s  -I-  0.9205  ±  i3.8975)(i  +  6.7488)(3  +  9  ±  ;6.2450)(5  +  50) 


(K  =.  3.5051  X  10®) 
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P733 


A:(s)(s  +  1.0554  ±i4.0816) 

(s  +  0.0144)(s+T0352  db  j4.1600)(s  •'  7.7061)(s  +  9  d:  i6.2450)(s  +  50) 
{K  =  4.0081  X  10®) 


P833 


K{s)(- -  ; '>973±i4.0793) 

(5  +  0.01 44)(s  +  1.0574  dr  ;4  'V  •.  ~i+  7.7056)(s  +  9  d:  j6.2450)(s  +  50) 


{K  =  4.0081  X  10®) 


P933  - 


Kis){s  ■  0.5961  d:j2.2934) 

(s  H-  0.0251)(s  +  0.6107  ±  ji./-  7  +  4  3763)(s  +  9  ±  i6.2450)(s  +  50) 


(A'  =  1.2610  X  10®) 


P1O33  — 


’i'(s)(s  -1-0  7521  ±j2.8941) 

(s  -1-  0.0201)(s  -J-  0.7723  ±  i2.96lO)(£  +  5.5062)(s  -f  9  ±  j6.2450)(s  -t-  50) 


(A  =  2.00b0  X  10®) 


Pllaa  = 


A(s)(s-f  0.92/3  ±j3.5301) 

(s  d-  0.0167)(5  -t-  0.9540  ±  i3.6026)(s  +  6.7772)(s  +  9  ±  y6.2450)(5  -t-  50) 


(A  =  2.9901  X  10®) 


Pi  2,3  = 


A(s)(5d- 0.3791  ±jl.5717) 

(s  -r  0.0333)(5  +  0.3819  d:  il.6351)(5  +  2.8124)(5  -f  9  d:  i6.2450)(5  +  50) 


(A  =  5.9151  X  10“) 


P13j3  = 


A(s)(s-l- 0.3691  ±;1.5709) 

(5  -t-  0.0334)(s  -I-  0.3720  d:  jl.6341)(5  +  2.8120)(s  -f  9  ±  ;6.2450)(s  d-  50) 


(A  -  5.9151  X  10“) 


P1433  = 


_ A(5)(sd- 9.4336  1.6769) _ 

(s  d-  0.0334)(s  d-  0.4413  ±  il.7419)(5  -f  3.2028)(8  -f  9  d:  j6.2450)(s  d-  50) 


(A  =  6.7651  X  10“) 


P1033 

A(«)(s  d-  0.4221  ±  j  1.67.59) 

(s  d-  0.0335)(.s  d-  0".4300  i  il.7407)(5  d-  3.2024)(s  d-  9  ±  i6.2450)(6  -f  50) 
(A  =  6.7651  X  10“) 
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Pl6aj  — 

/!:(s)(s  + 0.3909  ±il.5795) 

(s  +  0.0326)(s  +  0.3933  ±  il.6424)(s  +  2.8125)(s  +  9  ±  j6.2450)(s  +  50) 

{K  ~  5.9151  X  i0“) 

Pl~22  — 

ii:(s)(s  + 0.3808  ±  il.5787) 

(s  +  0.0328)(s 

+  0.3835  ±  il.6414)(s  +  2.8121)(s  +  9  i  i6.2450)^s  +  50) 

{K  =  5.9151  X  10“) 

Pl8aj  — 

/i:(5)(s  + 0.4471  ±jl.6849) 

(s  +  0  0327)-s 

0.4544  ±  il.7495)(s  +  3.2028)(s  +  9  ±  ■S.2450)(s  +  50) 

{K  =  6.7651  X  10“) 

Pi  922  = 

/(■(s)(s  + 0.4356  ±jl.6839) 

(s  +  0.0329)(s  +  0.4431  ±  il.7482)(s  +  3.2025)(s  +  9  ±  i6.2450)(s  +  50) 


(/<•  =  6.7651  X  10“) 
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Nominal  s-Domain  Transfer  Function 

Voi3  —  Pl23  “■ 

_ /i:(5)(5  +  3.9245)(s  -  83.4616) _ 

(5  +  0.0141)(s  +  0.9734  ±  i3.9181)(s  +  6.1996)(5  +  6.7492)(s  +  50) 

{K  =  -72.7664) 


Off-Nominal  s-Domain  Transfer  Functions  (■:  A 

\Orudder(v  / 

P233  — 

K{s)is  +  3.9202)(s  -  82.9677) 

(s  +  0.0142)(s  +  0.9490  ±  i3.9158)(s  +  6.1996)(s  +  6.7487)(s  +  50) 

{K  =  -72.7664) 


P333  = 


K(s){s  +  3.3710)(s  -  94.9224) 

(s  +  0.0141)(s  1.1179  ±  j4.1789)(s  +  6.1996)(s  +  7.7060)(s  +  50) 


{K  =  -83.2091) 


P433  = 


K{s)(s  +  3.9631)(s  -  94.3546) 

(s  -f-  0.0141)(s  +  1.0900  ±  i4.1761)(5  +  6.1996)(5  +  7.7056)(5  +  50) 


{K  =  -83.2091) 


PSjs  — 


/i:(s)(s  +  3.7525)(5  -  52.1859) 

(5  +  0.0144)(5  +  0.9448  ±  j3.8997)(5  +  6.y996)(s  +  6.7492)(s  +  50) 


{K  =  -113.3549) 


P6a3  “ 


_ /i'(5)(.5  +  3.7479)(s  -  51.8844) _ 

(5  +  0.0144)(s  +  0.9205  db  j3.8975)(s  +  6.1996)(s  +  6.7488)(s  +  50) 


(K  =  -113.3549) 


P7a3  = 


/C(s)(s  +  3.8071)(5  -  59.1913) 

(s  +  0.0144)(s  +  1.0852  ±  j4.1600)(s  +  6.1996)(5  +  7.7061)(s  +  50) 


(K  =  -129.6224) 


PS33  ~ 


{!< 


K{s){s  +  3.7992)(s  -  58.8138) 

(5  +  0.0144)(.?  +  1.0574  db  i4.1573)(5  +  6.1996)(s  +  7.7056)(s  +  50) 
-129.6224) 
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K{s){s  +  2.1658)(s  -  41.1729) 

(s  +  0.0^)(s+ 0.6107  ±  i2.3573)(s  +  4.3763)(s  +  6.1996)(s  +  50) 

{K  =  -33.4806) 


PlOj3 


K{s){s  +  2.7334)(s  -  51.9471) 

(s  +  0.0201)(s  +  0.7723  ±  j2.9610)(s  +  5.5062)(s  +  6.1996)(s  +  50) 


{K  =  -53.3141) 


P1I23 


K{s)(s  +  3.3076)(s  -  63.9864) 

(s  +  0.0167)(s  +  0.9540  ±  j3M26){s  +  6.1996)(s  +  6.7772)(5  +  50) 


{K  =  -79.3870) 


P1233  = 


K(s){s  +  1.5416)(s  -  21.4378) 

(5  +  0.0333)(s  +  0.3819  ±  jl.6351)(s  +  2.8124)(s  +  6.1996)(s  +  50) 


{K  =  -19.1297) 


P13j3  ~ 


i<:(s)(s  +  1.5397)(s  -  21.3139) 

(s  +  0.0334)(s  +  0.3720  ±  il.6341)(s  +  2.8120)(5  +  6.1996)(s  +  50) 


{K  =  -19.1297) 


P1433  = 


K{s){s  +  1.5640)(5  -  24.3194) 

(5  +  0.0334)(s  +  0.4413  ±  il.7419)(5  +  3.2028)(5  +  6.1996)(«  +  50) 


{K  =  -21.8786) 


P15j3  = 


A'(s)(s  +  1.5608)(5  -  24.1766) 

(s  +  0.0335)(s  +  0.4300  ±  jl.7407)(5  +  3.2024)(5  +  6.1996)(s  +  50) 


{K  =  -21.8786) 


P1623  — 


_ /i:(g)(a  +  1.6122)(s  -  34.2857) _ 

(5  +  0.0326)(s  +  0.3933  ±  jl.6424)(5  +  2.8125)(5  +  6.1996)(s  +  50) 


(K  =  -12.2800) 


PlTa,  = 


K{s)[s  +  1.6105)(s  -  34.0828) 

(s  +  0.0328)(s  +  0.3835  ±  il.6414)(5  +  2.8121)(5  +  6.1996)(s  +  50) 


{I<  =  -12.2800) 


PlS23  — 

PlQaa  — 

{K 


/r(s)(s  +  1.6314)(s-39) 

(s  +  0.0327)(s  +  0.4544  ±  jl.7495)(s  +  3.2028)(s  +  6.1996)(s  +  50) 

-14.0446) 


K{s){s  +  1.6281)(s  -  38.7667) 

(s  +  0.0329)(s  +  0.4431  ±  il.7482)(s  +  3.2025)(s  +  6.1996)(s  +  50) 

-14.0446) 


B-13 


Nominal  s-Domain  Transfer  Function  = 


Poii  ~  PI33  ~ 

K{s  -  1.8300)(s  +  2.3301)(s  +  22.9473) 

(s  +  0.0141)(5  +  0.9734  db  i3.9181)(s  +  6.7492)(s  +  9  ±  i6.2450)(s  +  50) 

{K  =  -9.5335  X  10^) 


Off-Nominal  s-Domain  Transfer  Functions  ( 

\‘>aiUron(V  J 

P233  — 

K{s  -  1.8520)(s  +  2.3024)(s  +  22.9479) 

(s  +  0.0142)(s  +  0.9490  ±  i3.9158)(s  +  6.7487)(s  +  9  ±  j6.2450)(s  +  50) 

(A'  =  -9.5335  X  10^) 


PSjj 


A'(s  -  1.8190)(s  +  2.3396)(s  +  26.2918) 

(s  +  0.0141)(s  +  1.1179  ±  j4.1789)(s  +  7.7060)(s  +  9  ±  j6.2450)(s  +  50) 


{K  =  -1.0902  X  lO"*) 


Pill  = 


A'(s  -  1 .8440)(s  +  2.3079)(s  +  26.2923) 

(f  +  0.0141)(s  +  1.0900  ±  i4.1761)(s  +  7.7056)(s  +  9  ±  ;6.2450)(s  +  50) 


(K  =  -1.0902  <  lO'*) 


PSjj  — 


K{s  -  1.8185)(s  +  2.3163)(s  +  22.9496) 

(«  +  0.0144)(s  +  0.9448  ±  j3.8997)(5  +  6.7492)(5  +  9  ±  j6.2450)(s  +  50) 


(A  =  -9.5335  X  10^) 


PSjj  — 


K{s  -  1.8405)(s  +  2.2885)(s  +  22.9502) 

(s  +  0.0144)(s  +  0.9205  ±  i3.8975)(s  +  6.7488)(5  +  9  db  j6.2450)(s  +  50) 


(A  =  -9.5.335  X  10^) 


PTjj  = 


A(s  -  1.8074)(s  +  2.3260)(s  +  26.2938) 

(5  +  0.0144)(s  +  1.0852  db  i4.1600)(s  +  7.7061)(s  +  9  db  i6.2450)(s  +  50) 


(A  =  -1.0902  X  lO'’) 


P833 


A(.s  -  1.8324)(s  +  2.2942)(s  +  26.2943) 

{s  +  0.0144)(5  +  1 .0574  ±  i4.1573)(s  +  7.70.56)(s  +  9  ±  j6.2450)(.s  +  50) 


(A  =  -1.0902  X  lO^*) 


B-M 


P93J  = 


K{s  -  1.8454)(s  +  2.3241)(s  +  14.6538) 

(s  +  0.0251)(s  +  0.6107  ±  j2.3573)(s  +  4.3763)(5  +  9  ±  i6.2450)(s  +  50) 

{K  =  -3.4299  X  10®) 


PlO,j  = 


K{s  -  1.8452)(s  +  2.3065)(s  +  18.6307) 

(s  +  0.0201)(s  +  C.7723  ±  j2.9610)(s  +  5.5062)(s  +  9  ±  j6.2450)(.'!  +  50) 


{K  =  -5.4617  X  10®) 


Pllsa 


I<{s  -  1.8371)(s  +  2.3067)(s+  23.0703) 

(s  +  0.0167)(s  +  0.9540  ±  i3.6026)(s  +  6.7772)(s  +  9  ±  j6.2450)(s  +  50) 


{K  =  -8.1327  X  10®) 


P1232  = 


K{s  -  1 .8105){s  +  2.4350){s  +  9.0075) 

(s  +  0.0333)(s  +  0.3819  ±jl.6351)(s  +  2.8124)(s  +  9  ±  j6.2450)(s  +  50) 


(/C  = -1.6089  X  10®) 


P133J  “ 


I<(s  -  1.8195)(s  +  2.4225)(5  +  9.0089) 

(s  +  0.0334)(s  +  0.3720  ±  jl.6341)(s  +  2.8120)(s  +  9  ±  i6.2450)(«  +  50) 


{K  =  -1.6089  X  10®) 


PHsj  = 


K{s  -  1.8200)(s  +  2.3884)(s  +  10.4478) 

(s  +  0.0334)(s  +  0.4413  ±  il.7419)(s  +  3.2028)(«  +  9  ±  j6.2450)(s  +  50) 


(K  =  -1.8401  X  10®) 


P15,3  = 


K{s  -  1.8304)(s  +  2.3746)(5+  10.4489) 

(s  4-  0.0335)(5  +  0.4300  ±  jl  .7407)(s  +  3.2024)(«  +  9  ±  j6.2450)(.5  +  50) 


{!<  =  -1.8401  X  10®) 


P1633  — 

_ K(s  -  1.8210)(5  +  2.4524)(g  +  9.0007) _ 

(«  4-  0.0326)(s  4-  0.3933  ±  il.6424)(«  4-  2.8I25)(s  4-  9  ±  j6.2450)(*  +  50) 

{!<  =  -1.6089  X  10®) 

Pl"33  ~ 

A  -  1 .8301  )(s  4-  2.4399)(5  4-  9.0020) 

(*  4-  0.0328)(.<  4-  0.3835  ±  jl  .64 14){.5  4-  2.8 121  )(5  +  9  ±  j6.2150)(  •  +  50) 

{K  =  -1.6089  X  10®) 


B-15 


P18j3  — 


K{s  ~  1.8308)(s  +  2.4048)(s  +  10.4421) 

(s  +  0.0327)(s  +  0.4544  ±  jl.7495)(s  +  3.2028)(s  +  9  ±  i6.2450)(5  +  50) 

{K  =  -1.8401  X  103) 


P1933  = 


_ K{s  -  1.8411)(5  +  2.3910)(s+  10.4433) _ 

(s  +  0.0329)(s  +  0.4431  ±  il.7482)(s  +  3.2025)(s  +  9  ±  j6.2450)(s  +  50) 


{K  =  -1.8401  X  103) 


B-16 


Nominal  s-Domain  Transfer  Function  P033  =  (grudllrC*)) 

P033  —  P133 

/f(s  + 0.0658  ±j0.3201) 

(s  +  0.0141)(s  +  0.9734  ±  j3.9181)(s  +  6.1996)(s  +  50) 

{K  =  5.4760  X  lO^) 

Off-Nominal  s-Domain  Transfer  Functions 


P233  = 


/<:(«  + 0.0542  ±i0.3222) 

(s  +  0.0142)(s  +  0.9490  ±  j3.9158)(s  +  6.1996)(s  +  50) 


{K  =  5.4760  X  103) 


P333  — 


/r(s  + 0.0774  ±i0.3175) 

(s  +  0.0141)(s  +  1.1179  ±  j4.1789)(s  +  6.1996)(s  +  50) 


(K  =  6.2618  X  103) 


P433  — 


/!:(«  + 0.0641  ±i0.3205) 

(s  +  0.0141)(s  +  1.0900  ±  j4.1761)(s  +  6.1996)(«  +  50) 


(K  =  6.2618  X  103) 


P533  = 


K(s  + 0.0637  ±j0.3156) 

{$  +  0.0144)(s  +  0.9448  ±  j3.8997)(s  +  6.1996)(s  +  50) 


{K  =  5.2531  X  10^) 


P633  = 


/s:(5  + 0.0524  ±j0.3176) 

(s  +  0.0144)(s  +  0.9205  ±  i3.8975)(s  +  6.1996)(s  +  50) 


{K  =  5.2531  X  10^) 


P733  — 


iffs  +  0.0749  dbjO.3132) 

(s  +  0.0144)(s  +  1.0852  ±  j4.1600)(s  +  6.1996)(s  +  50) 


{K  =  6.0069  X  10^) 


P833  = 


/<■(«  + 0.0620  ±j0.3160) 

(s  +  0.0144)(s  +  1.0574  ±  i4.1573)(s  +  6.1996)(s  +  50) 


{K  =  6.0069  X  10®) 


B-17 


P9«  = 


/!:(s  + 0.0303  ±i0.3224) 

(s  +  0.0251)(s  +  0.6107  ±  i2.3573)(s  +  6.1996)(s  +  50) 


(K  =  1.9300  X  10^) 


PlO,j  = 


/(:(s  + 0.0440  dbiO.3212) 

(s  +  0.0201)(s  +  0.7723  ±  i2.9610)(5  +  6.1996)(s  +  50) 


{K  =  3.0733  X  10®) 


P1I33 


iv:(s  + 0.0583  ±i0.3191) 

(s  +  0.0167)(s  +  0.9540  ±  i3.6026)(s  +  6.1996)(s  +  50) 


{K  =  4.5763  X  103) 


P1233 


/^(s  + 0.0108  ±i0.3200) 

(s  +  0.0333)(s  +  0.3819±  il.6351)(s  +  6.1996)(s  +  50) 


(K  =  886.5021) 


P1333  = 


{K  =  886.5021) 


/<■(«  + 0.0062  ±i0.3201) 

{s  f  0.0334)(s  +  0.3720  ±  il.6341)(s  +  6.1996)(s  +  50) 


P1433  = 


if(s  4- 0.0172  ±i0.3202) 

(s  +  0.0334)(s  +  0.4413  ±  jl.7419)(5  +  6.1990)(s  +  50) 


(K  =  1.0139  X  103) 


Pi  533  = 


1^(8  + 0.0120  ±j0.3204) 

0.0335)(s  +  0.4300  ±  il.7407)(5  +  6.1996)(s  +  50) 


{K  -  1.0139  X  10®) 


P1633  = 


(K  =  924.1210) 


K'(5  +  0.0112±i0.3247) 

(s  +  0.0326)(s  +  0.3933  ±  il.6424)(5  +  6.1996)(s  +  50) 


P1733  = 


_ /ir(g  +  0.0Q65±  jO.3249) _ 

(5  +  0.0328)(s  +  0.3835  ±  il.6414)(5  +  6.1996)(s  -h  50) 


(K  =  924.1210) 


B-18 


P1833 


i<r(s  + 0.0179  ±j0.3249) 

(s  +  0.0327)(s  +  0.4544  ±  il.7495)(s  +  6.1996)(«  +  50) 


(K 

P1933  = 

(K 


=  1.0569  X  lO^) 

/i:(s  + 0.0124  ±i0.3252) 

(s  +  0.0329)(s  +  0.4431  ±  il.7482)(s  +  6.1996)(s  +  50) 

=  1.0569  X  10®) 


B-19 


Nominal  w'-Domain  Transfer  Function 


Po2i  —  PI32 


K{W){W  +  0.9530  ±  i3.8464)(w'  -  100)(w'  +  115.8649)(w'  -  130.9057)(w'  +  678.4573) 
(w'  +  0.0141)(w'  +  0.9748  ±  i3.9197)(w'  +  6.7389)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 


{K  =  3.1346  X  10-^) 


OfF-Nominal  w'-Domain  Transfer  Functions 


P2«  = 


/C(w')(w'  +  0.9284  ±  i3.8445)(w'  -  100)(w'  +  115.8649)(w'  -  130.9057)(w'  +  678.4573) 
(w'  +  0.0142)(w'  +  0.9505  ±  i3.9175)(w'  +  6.7385)(w'  +  9.0106  ±  j6.2026)(w'  +  46.2117) 


{K  =  3.1346  X  10-^) 


PSaa  — 


/<(w')(w'  +  1.0900  ±  j4.1029)(w'  -  100)(w'  +  115.8017)(w'  -  131.0513)(w'  +  669.7337) 
(w'  +  0.0141)(w'  +  1.1198  ±  i4.1808)(w'  +  7.6908)(w'  +  9.0106  ±  j6.2026)(w'  +  46.2117) 

{K  =  3.6274  X  lO-'*) 


P4m  = 


K{W){yi'  +  1.0618  ±  j4.1005)(w'  -  100)(w'  +  115.8017)(w'  -  131.0513)(w'  +  669.7337) 
(w'  +  0.0141)(w'  +  1.0919  ±  j4.1780)(w'  +  7.6904)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 

{K  =  3.6274  X  lO-'*) 


PS33  — 


A:(w')(w'  +  0.9243  ±  i3.8275)(w'  -  100)(w'  +  115.8649)(w'  -  130.9057)(w'  +  678.4571) 
(w'  +  0,0144)(w'  +  0.9462  ±  i3.9014)(w'  +  6.7390)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 


{K  =  3.1346  X  10"^) 


P633  — 


K(w')(w'  +  0.8997  ±  j3.8256)(w'  -  100)(w'  +  115.8649)(w'  -  130.9057)(w'  +  678.4571) 
(w'  +  0.0144)(w'  +  0.9219  ±  i3.8992)(w'  +  6.7385)(w'  +  9.0106  ±  i6.2026)(w'  +  46.21 17) 


{K  =  3.1346  X  lO-'*) 


B-20 


P7a3 


A'(w')(w'  +  1.0571  ±  ;4.0834)(w'  -  100)(w'  +  115.8017)(w'  -  13L0513)(w'  +  669.7335) 
(w'  +  0.0144)(w'  +  1.0870  db  i4.1620)(w'  +  7.6909)(w'  +  9.0106  db  i6.2026)(w'  +  46.2117) 


{K  =  3.6274  X  10-4) 


P833 


K{vf'){v/'  +  1.0290  db  i4.0812)(w^  -  100)(w^  +  115.8017)(w'  -  131.0513)(w'  +  669.7335) 
(w'  +  0.0144)(w'  +  1.0592  db  i4.1593)(w'  -b  7.6904)(w'  +  9.0106  ±  j6.2026)(w'  +  46.2117) 


(K  =  3.6274  X  10-4) 


P933  = 


K{v/'){v/'  +  0.5964  ±  j2.2938)(w'  -  100)(w'  +  116.0205)(w'  -  130.5472)(w'  +  700.9747) 
(w'  +  0.0251)(w'  +  0.61 10  ±  i2.3577)(w'  +  4.3735)(w'  +  9.0106  ±  i6.2026)(w'  +  46.21 17) 

(K  =  1.0940  X  10-4) 


P1033 


Kiw'){v/'  +  0.7527  db  i2.8948)(w'  -  100)(w'  +  115.9461)(w'  -  130.7184)(w'  +  690.0289) 
(w'  +  0.0201)(w'  +  0.7729  ±  j2.96jL7)(w'  +  5.5006)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 

{K  =  1.7678  X  10-4) 


P1I33  ~ 


A(w')(w'  +  0.9285  ±  j3.5313)(w'  -  100)(w'  +  115.8626)(w'  -  I30.9109)(w'  +  678.1364) 
(w'  +  0.0167)(w'  +  0.9552  ±  j3.6039)(w'  +  6.7669)(w'  +  9.0106  ±  ;6.2026)(w'  +  46.2117) 


{K  =  2,6752  X  10-4) 


Pi  233  •” 


A(w')(w'  +  0.3792  ±  il.o718)(w'  -  100)(w'  +  116.1239)(w'  -  130.3092)(w'  +  716.8074) 
(w'  +  0.0333)(w'  +  0.3820  ±  jl.6352)(w'  +  2.8117)(w'  +  9.0106  db  i6.2026)(w'  +  46.2117) 

{K  =  5.0247  X  10-5) 


B-21 


P1323 

K{v/'){yt'  +  0.3691  ±  il.5711)(w'  -  100)(w'  +  116.1239)(w'  -  130.3092)(w'  +  716.8074) 
(w'  +  0.0334)(w'  +  0.3721  ±  il.6342)(w'  +  2.8113)(w'  +  9.0106  ±  j6.2026)(w'  +  46.2117) 

{K  =  5.0247  X  10-5) 

P1433  “ 

/f (w')(V  +  0.4338  ±  il.6770)( w'  -  100)(w'  +  116.0980)(w'  -  130.3689)(w'  +  712.7699) 
(w'  +  0.0334)(w'  +  0.4414  ±  jl.7420)(w'  +  3.2017)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 

(K  =  5.7776  X  10-5) 

P1533 

K{vi'){vi'  +  0.4222  ±  il.6761)(w'  -  100)(w'  +  116.0380)(w'  -  130.3689)(w'  +  712.7699) 
(w'  +  0.0335)(w'  +  0.4301  ±  il.7408)(w'  +  3.2013)(w'  +  9.0106  ±  j6.2026)(w'  +  46.2117) 

(K  =  5.7776  X  10-5) 

Pi633 

/^(wO(w^  -i-  0.3910  ±  jl.5796)(w^  -  100)(w'  -f  116.1239)(w^  -  130.3092)(w'  +  716.8075) 
(w'  +  0.0326)(w'  +  0.3934±il.6425)(w'  +  2.8117)(w'+  9.0106  ±j6.2026)(w'  +  46.2117) 

(K  =  5.0247  X  10-5) 

P1733  = 

/f  (w')(w'  +  0.3809  ±  il.5788)(w'  -  100)(w'  +  116.1239)(w'  -  130.3092)(w'  +  716.8075) 
(w'  +  0.0328)(w'  +  0.3836  ±  jl.6415)(w'  +  2.8113)(w'  +  9.0106  ±  j6.2026)(w'  +  46.2117) 

{K  -  5.0247  X  10-5) 

PI822  — 

K{W){yf'  4-  0.4472 ±  il.6850)(w'  -  100)(w'  +  116.0980)(w'  -  130.3689)(w'  +  712.7700) 
(w'  +  0.0327)(w'  +  0.4545  ±  j  1 .7496)(w'  +  3.2017)(w'  +  9.0106  ±  ;6.2026)(w'  +  46.21 17) 

{K  =  5.7776  X  10-5) 


B-22 


Nominal  w'-Domain  Transfer  Function  pojs  —  (grudifr'K)') 

Pojs  —  pi  23  ~ 

K{vi'){v/'  +  3.9225)(w'  -  58.1828)(w'  -  10'/j(w'  +  133.9606)(w'  -  360.3213) 
(w'  +  0.0141)(w'  +  0.9748  ±  j3.9197)(w'  +  6.1917)(w'  +  6.7389)(w'  +  46.2117) 

{K  =  -1.7033  X  10-5) 

Off-Nominal  w'-Domain  Transfer  Functions 

P233  — 

K{w'){w'  +  3.9182)(w'  -  67.9158)(w'  -  100)(w'  +  134.0174)(w'  -  358.7702) 
(w'  +  0.0142)(w'  +  0.9505  ±i3.9175)(w'  +  6.1917)(w'  +  6.7385)(w'  +  46.2117) 

{K  =  -1.7065  X  10-5) 

P323 

/<:(w')(w'  +  3.9689)(w'  -  73.9404)(w'  -  100)(w'  +  132.6917)(w'  -  399.6571) 
(w'  +  0.0141)(w'  +  1.1198  ±  j4.1808)(w'  +  6.1917)(w'  +  7.6908)(w'  +  46.2117) 

(ii' =  -1.8586  X  10-5) 

P423  = 

AXw')(w'  +  3.9610)(w'  -  73.6741)(w'  -  100)(w'  +  132.7492)(w'  -  397.6120) 
(w'  +  0.0141)(w'  +  1.0919  ±  i4.1780)(w'  +  6.l917)(w'  +  7.6904)(w'  +  46.2117) 

(A  =  -1.8629  X  10-5) 

PS33  — 

A(w')(w'  +  3.7507)(w'  -  47.8329)(w'  -  100)(w'  +  137.6547)(w'  -  284.7205) 
(w'  +  0.0144)(w'  +  0.9462  ±  j3.9014)(w'  +  6.1917)(w'  +  6.7390)(w'  +  46.2117) 

(A  =  -2.9123  X  10-5) 

P623  — 

A(w')(w'  +  3.7462)(w'  -  47.6015)(w'  -  100)(w'  +  137.7062)(w'  -  284.0332) 
(w'  +  0.0144)(w'  +  0.9219  ±  j3.8992)(w'  +  6.l917)(w'  +  6.7385)(w'  +  46.21 17) 

(A  =  -2.9155  X  10-5) 

Pr23  = 

A(w')(w'  +  3.8053)(w'  -  53.0185)(w'  -  100)(w'  +  136.4773)(w'  -  301.8984) 
(w'  +  0.0144)(w'  +  1.0870  ±  i4.1620)(w'  +  6.1917)(w'  +  7.6909)(w'  +  46.2117) 

(A  =  -3.2407  X  10-5) 


B-24 


P1933  * 

/$:(w')(w'  +  0.4357  ±  jl.6840)(w'  -  100)(w'  +  116.0980)(w'  -  130.3689)(w'  +  712.7700) 
(w'  +  0.0329)(w'  +  0.4433  ±  il.7484)(w'  +  3.2014)(w'  +  9.0106  4:  j6.2026)(w'  +  46.2117) 

{K  =  5.7776  X  10-5) 
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P^2i  — 


(w')(w'  +  3.7974)(w'  -  52.7699)(w'  -  100)(w'  +  136.5313)(w'  -  301.0226) 
(w'  +  0.0144)(w'  +  1.0592  ±  i4.1593)(w'  +  6.1917)(w'  +  7.6904)(w'  +  46.2117) 

{K  =  -3.2449  X  10"®) 


K{vi')ivi'  +  2.1655)(w'  -  38.9545)(w'  -  100)(w'  +  139.6278)(w'  -  261.8448) 
(w'  +  0.0251)(w'  +  0.6110  ±i2.3577)(w'  +  4.3735)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  -8.9330  X  10"°®) 

PlOjj 

/£'(w')(w'  +  2.7328)(w'  -  47.6504)(w'  -  100)(w'  +  137.7989)(w'  -  283.1456) 
(w'  +  0.0201)(w'  +  0.7729  ±  i2.9617)(w'  +  5.5006)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  -1.3747  X  10"®) 

Pilsj  = 

K(v/')(w'  +  3.3064)(w'  -  56.3575)(w'  -  100)(w'  +  136.0025)(w'  -  311.2883) 
(w'  +  0.0167)(w'  +  0.G552  ±  j3.6039)(w'  +  6.1917)(w'  +  6.7669)(w'  +  46.2117) 

{K  =  -1.9645  X  10"®) 

P12j3  ~ 

K(vi'){vi'  +  1.5414)(w'  -  21.1115)(w'  -  100)(w'  +  143.6447\w'  -  231.0659) 
(y  +  0.0333)(w'  +  0.3820  ±  il.6352)(w'  +  2.8117)(w'  +  6.101?)(w'  +  46.2117) 

{K  =  -5.4021  X  10-®®) 

PlSjj 

Ky){v/'  +  1.53':6Hw'  -  20.9932)(w'  -  100)(w'  +  143.6756)(v.'  -  230.8901) 
(w'  +  0.0334)(w'  v0.3721  ± ‘jl.6342)(w'  +  2.8113)(w'  +  6.1- 17)(w'  +  46  2117) 

{K  =  -5.4041  X  10-®®) 

P14m  — 

A' (w')(w'  +  1.5639)(w'  -  23.8448)(w'  -  100)(w'  +  142.9244)(vv'  -  235.3338) 
(w'  +  0.0334", Hw'  +  0.4414  ±jl.7420)(w'  +  3.2017)(w'  +  6.1£17)(v/'  +  46.2117) 

(A  = -6.1233  X  10-®®) 

P15j3  — 

A(w0(w'4-  1.5607)(w^-23.7103)(w'-  100)(w'+  142.9584)(w*  -  235.1224) 
(w'  +  0.0335)(w'  +  0.4301  ±  il.7408)(w'  +  3.2013)(w'  +  6.1917)(w'  +  46.2117) 

(A  =  -6.1259  X  10-°®) 

Pi  633 

A(w')(w'  +  1.6121)(w'  -  32.9818)(w'  -  100)(w'  +  141.1287)(w'  -  248.6719) 
(w'  +  0.0326)(w'  +  0.3934  ±  il.6425)(w'  +  2.81 17)(w'  +  6.1917)(w'  +  46.21 17) 

(A  =  -3.3575  X  10-®®) 
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Pi  733  = 


K{yi')(v/'  +  1.6103)(w'  -  32.8013)(w'  -  100)(v''  +  141.1695)(w'  -  248.3407) 
(w'  +  0.0328)(w'  +  0.3836  ±  jl.6415)(w'  +  2.8ll3)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  -3.3595  X  10“°®) 

P1833  — 

/'!:(w')(w'  +  1.6312)(w'  -  37.1036)(w'  -  100)(w'  +  140.1946)(w'  -  256.7633) 
(w'  +  0.0327)(w'  +  0.4545  ±  il.7496)(w'  +  3.2017)(w'  +  6.1917)(w'  +  46.2117) 

(/<■  =  -3.7852  X  10-“®) 

Pi  933  = 

/f(w')(w'  +  1.6279)(w'  -  36.9029)(w'  -  100)(w'  +  140.2385)(w'  -  256.3556) 
(w'  +  0.0329)(w'  +  0.4433  ±  il.7484)(w'  +  3.2014)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  -3.7879  X  10-“®) 
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Nominal  w'-Domain  Transfer  Function  ^<>32  = 


\  ^aileron  (w^)/ 


P032  —  Pl3J  — 

K{w'  -  1.8298)(w^  +  2.3297)(w'  +  22.5528)(w^  -  100)(w'  +  117.4175)(w'  -  128.0109)(w'  +  956.6060) 
(w'  +  0.0141)(w'  +  0.9748  ±  i3.9197)(w'  +  6.7389)(w'  +  9.0106  ±  j6.2026)(w'  +  46.21 17) 

{K  =  -6.2160  X  10-°®) 


Off-Nominal  w'-Domain  Transfer  Functions 

P23J  = 

K{vf'  -  1.8518)(w'  +  2.3020)(w'  +  22.5534)(w'  -  100)(w'  +  117.4175)(w'  -  128.0109)(w'  +  956.6068) 
(w'  +  0.0142)(w'  +  0.9505  ±  i3.9175)(w'  +  6.7385)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 

{K  =  -6.2160  X  10-°«) 


P33J  = 

K{vi'  -  1.8188)(w^  +  2.3392)(w'  +  25.7023)(w^  -  100)(w'  +  117.6049)(w'  -  127.7843)(w*  +  996.5365) 
(w'  +  0.0141)(w'  +  1.1198  d:  i4.1808)(w'  +  7.6908)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 

{K  =  -6.8587  X  10-°«) 


P<33  — 

K{w'  -  1.8438)(w'  +  2.3075)(w'  +  25.7028)(w'  -  100)(w'  +  117.6050)(w'  -  127.7844)(w'  +  996.5376) 
(w'  +  0.0141)(w'  +  1.0919  ±  j4.1780)(w'  +  7.6904)(w'  +  9.0106  ±  j6.2026)(w'  +  46.21 17) 


(K  =  -6.8587  X  10-°°) 

P5,3  = 

K(v/’  -  1.8183)( w'  +  2.3158)(w'  +  22.5550)(w'  -  100)(w'  +  1 17.4215){w'  -  128.0031)(w'  +  957.6502) 
(w'  +  0.0144)(w'  +  0.9462  ±  i3.9014)(w'  +  6.7390)(w'  +  9.0106  ±  i6.2026)(w'  +  46.21 17) 

{K  =  -6.2094  X  10"°®) 


P633  = 

K{w'  -  1.8402)(w'  -f  2.2881)(w'  +  22.5556)(w'  -  100)(w'  +  1 17.4215)(w'  -  128.0031)(w'  +  957.6509) 
(w'  +  0.0144)(w'  +  0.9219  ±  j3.8992)(w'  +  6.7385)(w'  +  9.0106  ±  i6.2026)(w'  +  46.21 17) 

(K  =  -6.2094  X  10"°®) 
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P733  = 


K{w'  -  1.8072)(w'  +  2.3255)(w'  +  25.7042)(w'  -  100)(w'  +  117.6096)(w'  -  127.7755)(w'  +  997.8334) 
(w'  +  0.0144)(w'  +  1.0870  ±i4.1620)(w'  +  7.6909)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 

{K  =  -6.8499  X  10-"®) 


P8j3  — 

Kjvi'  -  1.8322)(w'  +  2.2938)(w^  +  25.7047)(w'  -  100)(w'  +  117.6096)(w'  -  127.7755)(w'  +  997.8345) 
(w' +  0.0144)(w' +  1.0592  ±  j4. 1593)(w' +  7.6904)(w' +  9.0106  ±  i6.2026)(w' +  46.21 17) 


{K  =  -6.8499  X  10-°«) 


P9j3  = 

^■(w'  -  1.8452)(w'  +  2.3236)(w'  +  14.5498)(w'  -  I00)(w'  +  116.9860)(w'  -  128.6073)(w'  +  870.0983) 
(w'  +  0.0251)(w'  +  0.6110  ±  j2.3577)(w'  +  4.3735)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 

{K  =  -2.4320  X  10-"®) 


PlOjj  = 

K(vi'  -  1.8450)(w'  +  2.3061)(w^  +  18.4181)(w^  -  100)(w*  +  117.1871)(w'  -  128.3112)(w'  +  909.7099) 
(w'  +  0.0201)(w'  +  0.7729  ±  j2.9617)(w'  +  5.5006)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 


{K  =  -3.7217  X  10-"®) 


PIIm  = 


K{vi'  -  1.8369)(w'  +  2.3063)(w'  +  22.6695)(w'-  i00)(w'+  117.4252)(w'-  127.9967)(w'  + 958.4968) 
(w'  +  0.0167)(w'  +  0.9552  ±  i3.6039)(w'  +  6.7669)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 


(K  =  -5.2925  X  10-"®) 


P12m  = 

K(w'  -  1.8103)(w'  +  2.4345)(w'  +  8.9833)(w'  -  100)(w'  +  1 16.7232)(w'  -  129.0453)(w'  +  820.4626) 
(w'  +  0.0333)(w'  +  0.3820  ±  jl  .6352)(w'  +  2.81 17)(w'  +  9.0106  ±  j6.2026)(w'  +  46.21 17) 

(/f  = -1.2033  X  10-°®) 


B-28 


P13s3  ~ 


K{vi'  -  1.8193)(w'  +  2.4220)(w'  +  8.9846)(w'  -  100)(w'  +  116.7232)(w'  -  129.0453)(w'  +  820.4627) 
(w'  +  0.0334)(w'  +  0.3721  ±  il.6342)(w'  +  2.8113)(w'  +  9.0106  ±  j6.2026)(w'  +  46.2117) 

{K  -  -1.2033  X  10-°®) 


PUsa  = 

K{w'  -  1.8198)(w'  +  2.3880)(w'  +  10.4099)(w'  -  100)(w'  +  116.7876)(w'  -  I28.9316)(w'  +  832.4270) 
(w'  +  0.0334)(w'  +  0.4414  ±  il.7420)(w'  +  3.2017)(w'  +  9.0106  ±  j6.2026)(w'  +  46.2117) 

(K  =  -1.3581  X  10“°°) 


PlSsa  = 

K(w'  -  1.8302)(w'  +  2.3742)(w'  +  10.4110)(w'  -  100)(w'  +  116.7876)(w'  -  128.9316)(w'  +  832.4271) 
(w'  +  0.0335)(w'  +  0.4301  ±  il.7408)(w'  +  3.2013)(w'  +  9.0106  ±  j6.2026)(w'  +  46.2117) 

(K  =  -1.3581  X  10-°«) 


Pi6„  = 

if  (w'  -  1.8208)(w'  +  2.4520)(w'  +  8.9764)(w'  -  100)(w'  +  1 16.7217)(w'  -  129.0487)(w'  +  820. 1487) 
(w'  +  0.0326)(w'  +  0.3934  ±  il.6425)(w'  +  2.8117)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 

(if  =  -1.2038  X  10-°®) 


Pl7ja  = 


A'(w'-  1.8299)(w'  +  2.4395)(w'+8.9778)(w'-  100)(w'+  116.7217)(w'-  129.0487)(w'+ 820.1488) 
(w'  +  0.0328)(w'  +  0.3836  ±  il.6415)(w'  +  2.81 13)(w'  -f  9.0106  ±  j6.2026)(w'  +  46.21 17) 

(A' = -1.2038  X  10-°®) 


Pl«ja  = 

A'(w'  -  1 .8306)(w'  +  2.4044)(w'  +  10.4044)(w'  -  100)(w'  +  1 16.7858)(w'  -  128.9354)(w'  +  832.0572) 
(w'  +  0.0327)(w'  +  0.4545  ±  jl  .7496)(«'  +  3.2017)(w'  +  9.0106  ±  j6.2026){w'  +  46.21 17) 

(if  =  -1.3587  X  10-°®) 
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= 


K{vi'  -  1.8409)(w'  +  2.3906)(w'  +  10.4055)(w'  -  100)(w'  +  116.7858)(w'  -  128.9354)(w'  +  832.0573) 
(w'  +  0.0329)(w'  +  0.4433  ±  il.7484)(w'  +  3.2014)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2117) 

{K  =  -1.3587  X  10-°®) 
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Nominal  w'-Domain  Transfer  Function  P033  = 

P033  “  PI33 

K{v/'  +  0.0658  ±  i0.3201)(w'  -  100)(w'  +  149.0423)(w'  -  208.1660) 
(w'  +  0.0141)(w'  +  0.9748  ±i3.9197)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  0.0016) 

Off-Nominal  w'-Domain  Transfer  FYmctions 

K{vi'  +  0.0542  ±  i0.3222)(w'  -  100)(w'  +  149.0496)(w'  -  208.1436) 
(w'  +  0.0142)(w'  +  0.9505  ±  i3.9175)(w'  +  6.1917)(w'  +  46.2117)“ 

{K  =  0.0016) 

P3„  = 

K{v/'  +  0.0774  ±  j0.3175)(w'  -  100)(w'  +  148.9724)(w'  -  208.3999) 
(w'  +  0.0141)(w'  +  1.1198  ±  j4.1808)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  0.0019) 


K{yi'  +  0.0641  d:  i0.3205)(w'  -  100)(w'  +  148.9807)(w'  -  208.3742) 
(w'  +  0.0141)(w'  +  1.0919  ±i4.1780)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  0.0019) 


ii:(w' +  0.0637  ±j0.3156)(w' -  100)(w'  +  149.0583)(w'-  208.1158) 
(w'  +  0.0144)(w'  +  0.9462  ±  j3.9014)(w'  +  6.19l7)(w'  +  46.2117) 

(K  =  0.0016) 

P633  — 

K{v/'  +  0.0524  ±  j0.3176)(w'  -  100)(w'  +  149.0658)(w'  -  208.0927) 
(w'  +  0.0144)(w'  +  0.9219  ±  i3.8992)(w'  +  6.1917)(w'  +  46.2117) 

(K  =  0.0016) 

P'^33  — 

K{v/'  +  0.0749  ±  j0.3132)(w'  -  100)(w'  +  148.9906)(w'  -  208.3423) 
(w'  +  0.0144)(w'  +  1.0870  ±j4.1620)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  0.0018) 

PS33  — 

K{vi'  +  0.0620  ±  j0.3160)(w'  -  100)(w'  +  148.9992)(w'  -  208.3159) 
(w'  +  0.0144)(w'  +  1.0592  ±  j4.1593)(w'  +  6.1917)(w'  +  46.2117) 


{K  =  0.0018) 


PlSaa 


K{vt'  +  0.0179  ±i0.3249)(w'  -  100)(w'  +  149,2950)(w'  -  207.2924) 
(w'  +  0.0327)(w'  +  0.4545  ±  il.7496)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  3.1530  X  10"“) 

A’(w'  +  0.0124  ±  j0.3252)(w'  -  100)(w'  +  149.2985)(w'  -  207.2820) 
(w'  +  0.0329)(w'  +  0.4433  ±  il.7484)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  3.1531  X  10“^) 
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PSis  — 


/<(w'  +  0.0303  dbi0.3224)(w'  -  100)(w'  +  149.2166)(w'  -  207.5518) 
(w'  +  0.0251)(w'  +  0.6110  ±  i2.3577)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  5.7543  X  10-*) 


Pl0,3  = 

K{w'  +  0.0440  ±  i0.3212)(w'  -  100)(w'  +  149.1370)(w'  -  207.8213) 
(w'  +  0.0201)(w'  +  0.7729  ±i2.9617)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  9.1579  X  10-^) 


P1133  — 

K{vr'  +  0.0583  ±  i0.3191)(w'  -  100)(w'  +  149.0488)(w'  -  208.1273) 
(w'  +  0.0167)(w'  +  0.9552  ±i3.6039)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  0.0014) 

P12,3  = 

K{v/'  +  0.0108  ±i0.3200)(w'  -  100)(w'  +  149.3314)(w'  -  207.1789) 
(w'  +  0.0333)(w'  +  0.3820  ±jl.6352)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  2.6453  X  10“^) 

P1333  — 

K{vr'  +  0.0062  db  j0.3201)(w'  -  100)(w'  +  149.3345)(w'  -  207.1695) 
(w'  +  0.0334)(w'  +  0.3721  ±  il.6342)(w'  +  6.1917)(w'  +  46.2117) 

(K  =  2.6454  X  lO"**) 


+  0.0172  ±  i0.3202)(w'  -  100)(w'  +  149.3026)(w'  -  207.2692) 
(w'  +  0.0334)(w'  +  0.4414  ±il.7420)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  3.0248  X  lO-**) 


P1533  = 

/C(w'  +  0.0120  ±  i0.3204)(w'  -  100)(w'  +  149.3062)(w'  -  207.2584) 
(w'  +  0.0335)(w'  +  0.4301  ±il.7408)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  3.0249  X  10-“*) 


P1633 

K{v/'  +  0.0112  ±  i0.3247)(w'  -  100)(w'  +  149.3248)(w'  -  207.1992) 
(w'  +  0.0326)(w'  +  0.3934  db  il.6425)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  2.7574  X  10-'*) 


P1733  — 

Kivi'  +  0.0065  ±i0.3249)(vv'  -  100)(w'  +  149.3278)(w'-  207.1901) 
(w'  +  0.0328)(w'  +  0.3836  ±  jl.6415)(w'  +  6.1917)(w'  +  46.2117) 

{K  =  2.7575  X  lO"*) 
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Appendix  C.  MESO  Equivalent  Transfer  Functions 


This  appendix  includes  the  MISO  equivalent  transfer  functions,  and  33, 
used  in  the  MIMO  part  of  the  rate  controller  design. 


C.l  MISO  Equivalent  Transfer  Functions 

The  flight  conditions  produce  nineteen  separate  plants.  The  w'-domain  MISO 
equivalent  transfer  functions  used  in  the  roll  channel  of  the  lateral-directional  MIMO 
design  are  included  below: 

Nominal  w'-Domain  MISO  Equivalent  Transfer  Function 

?033e  —  Pi  33«  — 

K{vt'  +  0.0500)(w'  +  0.1523)(w'  -  100.)(w'  +  115.8656)(w'  -  130.9069)(w'  +  678.3976) 
(w'  +  0.0649  ±  ;0.3224)(w'  +  6.7663)(w'  +  9.0107  ±  ;6.2026)(w'  +  46.3698) 

(K  =  3.1345  X  10-“*) 

Off-Nominal  MISO  Equivalent  Transfer  Functions 

P211  = 

K(v/'  +  0.0500)(w'  +  0.1280)(w'  -  100.)(w'  +  115.8644)(w'  -  130.9069)(w'  +  678.3922) 
(w'  +  0.0529  ±  j0.3245)(w'  +  6.7272)(w'  +  9.0106  ±  j6.2027)(w'  +  46.3866) 

{K  =  3.1345  X  10-^) 


P3u  = 


K{vi'  +  0.0500)(w'  +  0.1737)(w'  -  100.)(w'  +  115.8012)(w'  -  131.0527)(w'  +  669.6585) 
(w'  +  0.0769  ±  i0.3198)(w'  +  7.7278)(w'  +  9.0107  ±  i6.2024)(w'  +  46.3038) 

{K  =  3.6273  X  10-“) 


P-lu  = 


K{w'  +  0.0500)(w'  +  0.1519)(w'  -  100.)(w'  +  115.8012)(w'  -  131.0527)(w'  +  669.6710) 
(w'  +  0.0632  ±  .j0.3227)(w'  +  7.7272)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2943) 


{K  =  3.6273  X  10-“) 
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PSii  = 

/C(w'  +  0.0500)(w'  +  0.1519)(w'  -  100.)(w'  +  115.8645)(w'  -  130.9077)(w'  +  678.3620) 
(w'  +  0.0626  ±  i0.3181)(w'  +  6.7270)(w'  +  9.0106  ±  i6.2027)(w'  +  46.2802) 

{K  =  3.1344  X  lO-'*) 

P6„  = 

K(vi'  +  0.0500)(w'  +  0.1239)(w'  -  100.)(w'  +  115.8646)(w'  -  130.9077)(w'  +  678.3434) 
(w'  +  0.0510  ±  i0.3200)(w'  +  6.7857)(w'  +  9.0107  ±  i6.2025)(w'  +  46.2752) 

{K  =  3.1344  X  10-^) 

P7„  = 

K{w'  +  0.0500)(w'  +  0.1684)(w'  -  100.)(w'  +  115.8008)(w'  -  131.0536)(w'  +  669.6284) 
(w'  +  0.0742  ±  i0.3156)(w'  +  7.7316)(w'  +  9  0107  ±  i6.2025)(w'  +  46.2864) 

{K  =  3.6272  X  10-4) 

P8u  = 

K{w'  +  0.0500)(w'  +  0.1519)(w'  -  100.)(w'  +  115.8008)(w'  -  131.0536)(w'  +  669.6118) 
(w'  +  0.0610  ±  j0.3183)(w'  +  7.7310)(w'  +  9.0107  ±  ;6.2024)(w'  +  46.3587) 

{K  =  3.6272  X  10-4) 

P9n  = 

Kiv/'  +  0.0500)(w'  +  0.0875)(w'  -  100.)(w'  +  116.0212)(w'  -  130.5482)(w'  +  700.9376) 
(w'  +  0.0237  ±  ;0.3244)(w'  +  4.3902)(w'  +  9.0106  ±  J6.2026)(w'  +  46.2751) 

{K=  1.0940  X  10-4) 

PlOii  = 

Kivi'  +  0.0500)(w'  +  0.1 104)(w'  -  100.)(w'  +  1 15.9461)(w'  -  130.7197)(w'  +  689.9580) 
(w'  +  0.0426  ±  i0.3233)(w'  +  5.5255)(w'  +  9.0107  ±  j6.2025)(w'  +  46.2613) 

{K  =  1.7678  X  10-4) 
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Pllll  = 

K{y/'  +  0.0500)(w'  +  0.1361)(w'  -  100.)(w'  +  115.8625)(w'  -  130.9126)(w'  +  678.0630) 
(w'  +  0.0570  ±  j0.3217)(w'  +  6.7545)(w'  +  9.0107  ±  j6.2027)(w'  +  46.2947) 

{K  =  2.6751  X  10“'*) 

Pl2,i  = 

K{vt'  +  0.0500)(w'  +  0.0605)(w'  -  100.)(w^  +  116.1249)(w'  -  130.3100)(w'  +  716.7692) 
(w'  +  0.0089  ±  j0.3216)(w'  +  2.8146)(w'  +  9.0106  ±  i6.2025)(w'  +  46.3636) 

{K  =  5.0247  X  10-5) 

Pl3„  = 

K{v/'  +  0.0500)(w'  +  0.0509)(w'  -  100.)(w'  +  116.1235)(w'  -  130.3100)(w'  +  716.7925) 
(w'  +  0.0042  ±  i0.3216)(w'  +  2.8143)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2587) 

(K  =  5.0247  X  10-5) 

Pl4,i  = 

/C(w'  +  0.0500)(w'  +  0.0691)(w'  -  100.)(w'  +  116.0976)(w'  -  130.3698)(w'  +  712.7262) 
(w'  +  0.0154  ±  i0.3219)(w'  +  3.2090)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2723) 

{K  =  5.7775  X  10-5) 

Pl5„  = 

K{w'  +  0.0500)(w'  +  0.0582)(w'  -  100.)(w'  +  1 16.0976)(w'  -  130.3698)(w'  +  712.7281) 
(w'  +  0.0101  ±  j0.3221)(w'  +  3.2087)(w'  +  9.0106  ±  j6.2026)(w'  +  46.2707) 

(K  =  5.7775  X  10-5) 

Pl6,,  = 

K{vf'  +  0.0500)(w'  +  0.0624)(w'  -  100.)(w'  +  116.1237)(w'  -  130.3097)(w'  +  716.7915) 
(w'  +  0.0093  ±  j0.3264)(w'  +  2.8142)(w'  +  9.0106  ±  j6.2025)(w'  +  46.2756) 

(K  =  5.0247  X  10-5) 
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Pl7n  = 

K{vi'  +  0.0500)(w'  +  0.0525)(w'  -  100.)(w'  +  116.1237)(w'  -  130.3097)(w'  +  716.8044) 
(w'  +  0.0045  ±  i0.3264)(w'  +  2.8139)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2768) 

{K  =  5.0247  X  10-®) 

PlSii  — 

K{w'  +  0.0500)(w'  +  0.0714)(w'  -  100.)(w'  +  116.0977)(w'  -  130.3694)(w'  +  712.7576) 
(w'  +  0.0160  ±  i0.3267)(w'  +  3.2083)(w'  +  9.0106  ±  j6.2026)(w'  +  46.3536) 

{K  =  5.7775  X  10-®) 

Pl9ii  = 

K{vr'  +  0.0500)(w'  +  0.0601)(w'  -  100.)(vv'  +  116.0994)(w'  -  130.3694)(w'  +  712.7559) 
(w'  +  0.0105  ±  j0.3269)(w'  +  3.2079)(w'  +  9.0106  ±  i6.2026)(w'  +  46.2643) 

(K  =5.7775  X  10-®) 
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C.2  qi33  MISO  Equivalent  Transfer  Functions 

The  flight  conditions  produce  nineteen  separate  plants.  The  w'-domain  MISO 
equivalent  transfer  functions  used  in  the  yaw  channel  of  the  lateral-directional  MIMO 
design  are  included  below: 

Nominal  w'-Domain  MISO  Equivalent  Transfer  Function 

9o33e  —  Plaje  ~ 

A:(w'  +  0.1519)(w'  -  100.)(w'  +  149.0505)(w'  -  208.1358) 

(w'  +  0.9748  ±  i3.9197)(w'  +  6.7389)(w'  +  46.2122) 

{K  =  0.0016) 


Off-Nominal  MISO  Equivalent  Transfer  Functions 


P2il 

{K  =  0.0016) 
P3n  = 

{K  =  0.0019) 
P4ii  = 

[K  =  0.0019) 

P5,i  = 

{K  =  0  0016) 
P6ii  = 

{K  -  0.0016) 


K{vi'  +  0.1279)(w'  -  100.)(w'  +  149.0577)(w'  -  208.1136) 
(w'  +  0.9505  ±  j3.9175)(w'  +  6.7385)(w'  +  46.2120) 


K{yi'  +  0.1737)(w'  -  100.)(w'  +  148.9816)(w'  -  208.3650) 
(w'  +  1.1198  ±  i4.1808)(w'  +  6.1917)(w'  +  46.2129) 


A'(w'  +  0.1462)(w'  -  100  )(w'  +  148.9898)(w'  -  208.3395) 
(w'  +  1.0919  ±  j4.1780)(w'  +  6.1917)(w'  +  46.2120) 


K(yi'  +  0.1471)(w'  -  100.)(w'  +  149.0668)(w'  -  208.0844) 
(w'  +  0.9462  ±  i3.9014)(w'  +  6.7390){w'  +  46.2127) 


/<’(w'  +  0.1239)(w'-  100.)(w'+  149.0743)(w'-  208.0615) 
(w'  +  0.9210  d-  j3.8992)(w'  +  6.7385)(w'  +  46.2121) 
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P7„  = 


{K  =  0.0018) 


{K  =  0.0018) 


K{vi'  +  0.1682)(w'  -  100.)(w'  +  149.0003)(w'  -  208.3061) 
(w'  +  1.0870  ±  i4.1620)(w'  +  6.19i7)(w'  +  46.2120) 


K{w'  +  0.1416)(w'  -  100.)(w'  +  149.0088)(w'  -  208.2798) 
(w'  +  1.0592  ±i4.1593)(w'  +  6.19r(')(w'  +  46.2127) 


K{w'  +  0.0875)(w'  -  I00.)(w'  +  149.2221)(w'  -  207.5327) 
(w'  +  0.6110  ±i2.3577)(w'  +  6. 1917)(w'  +  46.2128) 


{K  =  5.7523  X  10-“*) 


Pio„  = 


/i:(w'  +  0.1103)(w'-  100.)(w'+  149. 1439)(w'-  207.7969) 
(w'  +  0.7729  ±  i2.9617)(w'  +  6.1917)(w'  +  46.2126) 


{K  =  9.1550  X  lO-'*) 


Pllii  = 


K{vi'  +  0.1361)(w'  -  100.)(w'  +  149.0571)(w'  -  208.0966) 
(w'  +  0.9552  ±  i3.6039)(w'  +  6.7669)(w'  +  46.2121) 


(K  =  0.0014) 


Pl2u  = 


A'(w'  +  0.0604)(w'  -  100.)(w'  +  149.3351)(w'  -  207.1668) 
(w'  +  0.3820  ±  jl.6352)(w'  +  6.1917)(w'  +  46.2133) 


{K  =  2.6440  X  10-'’) 


Pl3,,  = 


K{vi'  +  0.0509)(w'  -  100.)(w'  +  149.3382)(w'  -  207.1575) 
(w'  +  0.3721  ±  jl.6342)(w'  +  6.1917)(w'  +  46.2120) 


(A  =  2.6441  X  10-'') 


K{w'  +  0.0691)(w'  -  100.)(w'  +  149.3069)(w'  -  207.2552) 
(w'  +  0.4414  ±il.7420)(w'  +  6.1917)(w'  +  46.2121) 


{K  =  3.0234  X  10-'') 


/\  (w'  +  0.0582)(w'  -  100.)(w'  +  149.3104)(w'  -  207,2446) 
(w'  +  0.4301  ±  il.7408)(w'  +  6.l917)(w'  +  46.2125) 


{K  =  3.0234  X  10-'') 


Pl6n  “ 


Kiv/'  +  0.0624)(w'  -  100.)(w'  +  149.3283)(w'  -  207.1876) 
(w'  +  0.3934  ±  jl.6425)(w'  +  6.1917)(w'  +  46.2123) 


{K  =  2.7566  X  10-4) 


Pirn  = 


K{w'  +  0.0525)(w'  -  100.)(w'  +  149.3313)(w'  -  207.1785) 
(w'  +  0.3836  ±  jl.6415)(w'  +  6.1917)(w'  +  46.2120) 


{K  =  2.7567  X  10-4) 


Pl8ii  = 


Kivi'  +  0.0714)(w'  -  100.)(w'  +  149.2991)(w'  -  207.2790) 
(w'  +  0.4545  ±  il.7496)(w'  +  6.1917)(w'  +  46.2124) 


{K  =  3.1521  X  10-4) 


Pl9n  — 


K{v/'  +  0.0601)(w'  -  100.)(w'  +  149.3025)(w'  -  207.2687) 
(w'  +  0.4433  ±  jl.7484)(w'  +  6.1917)(w'  +  46.2120) 


(/if  =  3.1522  X  10-4) 
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Appendix  D.  Template  Data 


This  appendix  includes  the  phcise  and  magnitude  information  used  to  generate 
the  plant  templates.  Templates  are  produced  for  P,„,  and  %3. 


D.l  Template  Data  for  Longitudinal  SISO  System  P,j, 


1  Template  Data  for  u  =1.0000x10”'’ trad/sec|  | 

1  Template  Data  for  v  =0.0050  [rad/sec|  | 

Flight  Condition 

B  3i  m  iTXffm  ffWTTB 

1 

451.0064 

455.0611 

2 

451.0766 

455.4127 

3 

450.9358 

454.7075 

-41.5528 

4 

451.0241 

455.1499 

-43.6898 

5 

451.0132 

-52.8514 

455.0949 

-38.7571 

6 

451.0884 

-54.8535 

455.4712 

-40.7549 

7 

450.9388 

-51.1476 

454.7224 

-37.0575 

8 

451.0320 

-53.3126 

1  8 

455.1891 

-39.2172 

9 

450.5065 

452.5526 

-49.9785 

10 

450.7775 

■KSMn 

10 

453.9150 

-47.5968 

11 

450.9371 

11 

454.7149 

WKBSSMMtii 

12 

90.2220 

12 

91.1196 

13 

90.1883 

HKSSSQHi 

13 

90.9497 

-49.1445 

14 

90.2507 

14 

91.2642 

15 

90.2143 

15 

91.0806 

16 

90.2163 

16 

91.0909 

17 

90.1836 

17 

90.9258 

-53.4042 

18 

90.2443 

18 

91.2318 

-53.0714 

19 

90.2088 

19 

91.0530 

II  Template  Data  for  v  =0.0100  [rad/sec]  | 

i  Template  Data  for  v  =0.0500  [rad/seci  I 

iaiiSalESI 

1  Flight  Condition 

1 

460.0183 

-37.1092 

490.3924 

-20.2859 

2 

460.7050 

-39.0696 

2 

492.5161 

-22.0410 

3 

459.3264 

-35.4391 

3 

488.2022 

-18.8280 

4 

460.1918 

-37.5609 

490.9300 

-20.6850 

5 

460.0827 

-32.6279 

5 

490.4427 

-15.7258 

6 

460.8173 

-34.6125 

6 

492.7055 

-17.4967 

7 

459.3539 

7 

488.1333 

-14.2558 

8 

460.2667 

-33.0847  1 

8 

491.0115 

-16.1285 

9 

455.0868 

9 

473.9133 

-29.1450 

10 

457.7794 

-41.5229  1 

10 

484.0482 

-25.7826 

11 

4.59.3461 

-38.0890  1 

11 

488.7716 

-21.6523 

12 

92.2348 

12 

101.0223 

-28.8213 

13 

91.8959 

HBEESESSHlj 

13 

99.3893 

-29.0436 

14 

92.5231 

TTttFFMI 

14 

102.3949 

-28.5955 

15 

92.1571 

15 

100.6492 

-28.8766 

16 

92.1777 

16 

100.7510 

-.33.1071 

17 

91.8483 

'mamamm 

99.1602 

■BOs&liiaHII 

18 

92.4586 

■EKMI 

18 

102.0931 

-32.8908 

19 

92.1021 

■Bo^an 

19 

100.3872 

-33.1599  1 
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1  Template  Data  for  v  =0.1000  {rad/sec]  | 

1  Template  Data  for  v  =0.2000  [rad/secj  || 

Flight  Condition 

BUiirnBrnB 

Magnitude  [dB] 

Flight  Condition 

Magnitude  [dB] 

1 

506.5487 

-9.1847 

1 

496.4792 

11.1798 

2 

509.0584 

-10.8101 

2 

503.4635 

9.5797 

3 

504.0086 

-7.8789 

3 

489.7485 

12.3330 

4 

507.1652 

-9.5495 

4 

497.9968 

10.8122 

5 

506.1464 

-4.4537 

5 

486.3716 

17.5230 

6 

508.8676 

-6.1094 

6 

496.7083 

15.7863 

7 

503.3913 

-3.1227 

7 

476.2673 

18.6225 

8 

506.8083 

8 

488.6013 

17.1076 

9 

491.6471 

-21.0033 

9 

510.9488 

-9.4540 

10 

502.9172 

-16.1772 

10 

516.7388 

-2.0770 

11 

506.3262 

-11.1374 

11 

512.0757 

5.7209 

12 

111.4676 

-22.1873 

12 

129.4728 

-13.8549 

13 

108.5054 

-22.5555 

13 

125.2952 

-14.6555 

14 

113.8739 

-21.8282 

14 

132.5298 

-13.1423 

15 

110.7954 

-22.2787 

15 

128.5211 

-14.0506  1 

16 

110.9964 

-26.5306 

16 

128.8913 

-18.3980 

17 

108.09(8 

-26.8786 

17 

124.7165 

-19.1536 

18 

113.3714 

-26.1891 

18 

131.9918 

-17.7196 

19 

110.3377 

-26.6164 

19 

127.9562 

-18.5803 

Template  Data  for  v  =0.4000  [rad/sec] 


Template  Data  for  v  =0.3000  [rad/sec] 


Phase  [deg] 


Magnitude  [dB] 


Phase  [deg] 


Magnitude  [dB] 


Flight  Condition 


Flight  Condition 


365.4678 


9.8066 


1 


359.8305 


5.9965 


364.9302 


8.7091 


360.3961 


4.7916 


366.2341 


10.5360 


359.7401 


6.9041 


365.6140 


9.5837 


360.1761 


5.7313 


362.7797 


12.7658 


358.2594 


9.4800 


362.5770 


11.9180 


358.9374 


8.4094 


363.3796 


13.3692 


358.1928 


10.2750 


363.1142 


12.5785 


358.7452 


9.2456 


519.6604 


1.9078 


488.3650 


30.9604 


10 


497.9735 


18.7854 


10 


358.8061 


10.4306 


11 


372.3707 


15.3369 


11 


359.2303 


8.3145 


12 


144.0963 


-6.4697 


12 


156.4607 


1.3931 


13 


140.3111 


-7.7293 


13 


153.4548 


-0.3208 


14 


146.4553 


-5.4260 


14 


157.6146 


2.7774 


15 


143.0991 


-6.7770 


15 


155.0899 


0.9651 


16 


143.4982 


-11.3478 


16 


154.9124 


-4.1379 


17 


18 


19 


139.6959 


-12.5282 


145.9830 


-10.3634 


142.5842 


-11.6284 


17 


18 


19 


152.1620 


156.2506 


153.9532 


-5.6717 

-2.9051 

-4.5049 


1  Template  Data  for  v  =50  [rad/sec]  | 

1  Template  Data  for  v  =100  [rad/sec]  (1 

Flight  Condition 

Esnsasgni 

Plight  Condition 

Phase  [deg] 

iBsssassii 

1 

45.8746 

-33.1618 

1 

2.4878 

-50.1788 

2 

45.7741 

-33.1633 

2 

2.4368 

-50.1821 

3 

46.4182 

-31.9618 

3 

2.7825 

-48.9989 

4 

46.2990 

-31.9630 

4 

2.7229 

5 

45.2562 

-33.5169 

5 

2.1597 

-50.4719 

6 

45.1906 

-33.5209 

6 

2.1262 

-50.4765 

7 

45.6964 

-32.3283 

7 

2.4047 

8 

45.6178 

-32.3327 

8 

2.3652 

^HEEESSSiil 

9 

44.3382 

9 

1.6524 

-59.2470 

10 

44.8797 

-38.2592 

10 

1.9580 

-55.1912 

11 

45.5107 

-34.7465 

11 

2.3061 

-51.7132 

12 

43.5168 

-49.0937 

12 

1.1767 

13 

43.4947 

-49.0984 

13 

1.1639 

-65.9751  1 

14 

43.6835 

-47.9197 

14 

1.2746 

15 

43.6576 

-47.9251 

15 

1.2597 

16 

43.7371 

-48.8183 

16 

1.3048 

-65.7053  1 

17 

43.7035 

-48.8220 

17 

1.2855 

18 

43.9363 

18 

1.4212 

MBSSEQHli 

19 

43.8973 

-47.6471  1 

19 

1.3988 
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D.2  Template  Data  for  Lateral-Directional  MIMO  System 


1  Template  Data  for  ti  =1.0000x10"'’ [rad/secj  | 

1  Template  Data  for  v  =0.0050  [rad/sec]  || 

Plight  Condition 

Magnitude  [dB] 

Plight  Condition 

tsmmM 

Magnitude  [dB] 

1 

-183.6963 

-16.2144 

1 

-197.8654 

-15.7056 

2 

-183.5840 

-17.6546 

wtmEmmm 

-197.3086 

-17.1470 

3 

-183.7780 

-15.2906 

3 

-198.2714 

-14.7812 

4 

-183.6729 

-16.7719 

4 

-197.7523 

-16.2648 

5 

-183.6000 

-16.1040 

5 

-197.4227 

-15.6130 

6 

-183.4840 

-17.5659 

6 

-196.8474 

-17.0760 

7 

-183.6827 

-15.8949 

-197.8344 

-15.4035 

8 

-183.5759 

-17.3111 

8 

-197.3056 

-16.8217 

9 

-181.5717 

-21.4814 

9 

-187.7862 

-21.2950 

10 

-182.3100 

-19.3960 

10 

-191.3732 

-19.1272 

11 

-183.0014 

-18.1925 

11 

-194.6621 

-17.8203 

12 

-180.6445 

-25.1780 

12 

-183.2049 

-25.0411 

13 

-180.4235 

-26.6121 

13 

-182.1025 

-26.4596 

14 

-180.7842 

-24.1569 

14 

-183.9043 

-24.0299 

15 

-180.5879 

-25.5838 

15 

-182.9229 

-25.4451 

16 

-180.7180 

-25.3050 

16 

-183.5700 

-25.1673 

17 

-180.5042 

-26.7476 

17 

-182.5031 

-26.5954 

18 

-180.8530 

-24.2779 

18 

-184.2460 

-24.1495 

19 

-180.6628 

-25.7200 

19 

-183.2948 

-25.5808 

1  Template  Data  for  v  =0.0100  [rad/sec)  | 

1  Template  Data  for  v  =0.0500  [rad/sec]  || 

Plight  Condition 

Plight  Condition 

laisggmi 

Magnitude  [dB]  || 

1 

-211.9206 

-14.4141 

1 

-237.4629 

2 

-210.8502 

-15.8550 

-233.2268 

3 

-212.7059 

-13.4901 

3 

-240.7390 

-3.3048  II 

4 

-211.7145 

-14.9763 

4 

-236.7608 

5 

-211.2159 

-14.3580 

5 

-236.4382 

-4.0912 

6 

-210.1098 

-15.8202 

6 

-232.1125 

-5.3386 

7 

-212.0139 

-14.1496 

7 

-239.7550 

-4.0280 

8 

-211.0010 

-15.5696 

8 

-235.6969 

-5.2834 

9 

-194.7545 

-20.7528 

9 

-212.3706 

-12.6818 

10 

-201.2093 

-18.3758 

10 

-223.6120 

-9.5362 

11 

-206.8336 

-16.8243 

11 

-232.0372 

-7.2664 

12 

-186.1339 

-24.6302 

12 

-194.1450 

-16.9202 

13 

-184.0060 

-26.0023 

13 

-188.8801 

-17.5870 

14 

-187.4991 

-23.6490 

14 

-198.1251 

-16.4477 

15 

-185.5910 

-25.0290 

15 

-192.8136 

-17.2093 

16 

-186.8311 

-24.7546 

16 

-195.6982 

-17.0889 

17 

-184.7673 

-26.1396 

17 

-190.4239 

-17.7912 

18 

-188.1530 

-23.7649 

18 

-199.6693 

-16.5884 

19 

-186.3008 

-25.1635 

19 

-194.3699 

-17.3930 

D-7 


Template  Data  for  v  =0.1000 


Flight  Condition  Phiise  [deg]  Magnitude 


-233.6027  3.7009 


-227,3772  2.8763 


-238.7709  4.1769 


-232.5140  3.2122 


-232.4350  3.7947 


-226.1861  2.9722 


-237.5786  3.5341 


-231.2864  2.6598 


-208.6173  -3.8827 


-219.2547  -1.17C2 


-228.1300  0.7124 


-193.0056  -7.4680 


-188.5884  -7.7129 


-196.7192  -7.3699 


-191.9003  -7.6745 


-194.1949  -7.7205 


-189,6905  -7.9902 


-197.9737  -7.5996 


-193.0725  -7.9396 


Template  Data  for  v  =0.2000  (rad/sec 


Flight  Condition  Phase  [deg]  Magnitude  [dB 


-233.8232  15.5879 


-225.9895  15.4074 


-240.7399  15.4551 


-232.2900  15.1950 


-233.1271  15.8983 


-225.3309  15.7182 


-239.8901  15.0255 


-231.4838  14.8597 


-208.5504  9,7262 


-218.7133  11.8289 


-228.2552  13.0486 


-191.7414  6,8247 


-190.8458  6.7806 


-198.4580  6.6956 


-194.0588  6.6449 


-195.3866  6.4049 


-191.4437  6.3468 


-199.1298  6.2904 


-194.6796  6.2158 


Template  Data  for  v  =0.3000  [rad/sec] 


Flight  Condition  Phase  [deg]  Magnitude 


-274.8380  26.5272 


-267.4827  27.7744 


-280.3593  25.1738 


-273.1934  26.3378 


-277.6025  26.9561 


-270.8736  28.2655 


-275.8648 


-249.8870 


-262.5350 


-271.6581 


-222.6590 


-209.3516 


-235.1797 


-223.4363 


-218.2281 


-206.5229 


-229.9543 


-218.8768 


26.1303 


25.8330 


25.6124 


24.8462 


27.5537 


28.3509 


25.8709 


27.0086 


25.8120 


26.3559 


24.5052 


25.3427 


Template  Data  for  o  =0.4000  [rad/sec 


1 

-338.0638 

2 

-340.9402 

3 

-335.3087 

-337.6750 


-340.1445 


-343.0338 


-339.5341 


-352.2171 


-345,6565 


-340.0357 


-365.6071 


-368.1115 


-360.4935 


-363.2823 


-365.1232 


-367.9133 


-359.6124 


-362.7001 


Magnitude 


26.2476 


26.9106 


25.3054 


25.9980 


26.0936 


26.6631 


24.4972 


25.1889 


22.2017 


23.5912 


23.9366 


19.3073 


19.3764 


19.1019 


19.2300 


19.7573 


19.8302 


19.5390 


19.6735 


D-8 


Flight  Condition 


Phase 


-396.6341 


-396.7455 


-393.2318 


-393.2808 


-397.1074 


-397.1964 


-392.8044 


-393.0016 


-396.0515 


r  V  =2 

[rad/sec]  || 

1  Template  Data  for  v  =3  [rad/sec] 

1  Flight  Condition  | 

1  Phase  [deg]  ! 

1  Magnitude  [dB 

r 

-414.6380 


-414.7145 


-410.6723 


-410.7610 


18.9464 


18.9596 


18.7413 


18.7097 


18.9318 


18.9198 


18.1346 


18.1588 


12.7686 


15.0486 


16.4095 


8.8804 


8.8970 


9.1044 


9.1235 


8.8697 


8.8781 


9.1028 


9.1056 


18.2702 


18.1824 


18.1517 


18.2402 


18.2278 


17.6569 


17.6745 


11.6307 


14.1542 


15.7560 


7.1873 


7.2016 


7.6473 


7.6635 


7.1741 


7.1801 


7.6435 


7.6434 


ISSIII 

iBSEEaii 


Template  Data  for  v  =4  (rad/sec 


Flight  Condition  Phase  [deg]  Magnitude 


17.4577 


17.4700 


17.5449 


17.5168 


17.4310 


17.4187 


17.0635 


17.0760 


10.3602 


13.1212 


14.9678 


5.4559 


5.4686 


6.0960 


6.1102 


5.4422 


5.4465 


6.0920 


6.0898 


Template  Data  for  v  =5  [rad/sec 


Flight  Condition  Phase  [deg]  Magnitude 


-444.1831 


■aBWir.mii 


maiaifmmt 


16.5443 


16.5567 


16.8093 


16.7847 


16.5020 


16.4901 


16.3244 


16.3333 


9.0247 


11.9826 


14.0410 


3.7875 


3.7992 


4.5497 


4.5625 


3.7735 


3.7769 


4.5457 


4.5421 


D.3  Template  Data  for  Lateral- Directional  MIMO  System 


1  Template  Data  for  v  =1.0000x10“^  [rad/sec)  | 

1  Template  Data  for  v  =0.0050  [rad/sec]  || 

Flight  Condition 

Magnitude  [dB] 

Flight  Condition 

Magnitude  [dB] 

1 

0.3593 

-15.3414 

1 

1.8114 

-15.3368 

2 

0.4303 

-16.8028 

mamamm 

2.1690 

-16.7963 

3 

0.3120 

-14.2072 

3 

1.5730 

-14.2036 

4 

0.3741 

-15.6638 

4 

1.8859 

-15.6589 

5 

0.3717 

-15.8706 

5 

1.8736 

-15.8657 

6 

0.4450 

-17.3327 

6 

2.2428 

-17.3258 

7 

0.3228 

-14.7357 

7 

1.6275 

-14.7320 

8 

0.3869 

-16.1932 

8 

1.9505 

-16.1879 

9 

0.6320 

-20.2703 

9 

3.1835 

-20.2565 

10 

0.4986 

-18.2516 

10 

2.5129 

-18.2429 

11 

0.4022 

-16.4373 

11 

2.0274 

-16.4316 

12 

0.9204 

-23.5971 

12 

4.6309 

-23.5681 

13 

1.0988 

-25.0588 

13 

5.5235 

-25.0180 

14 

-22.4610 

14 

4.0342 

-22.4388 

15 

0.9574 

-23.9182 

15 

4.8163 

-23.8869 

16 

0.8903 

-23.0680 

16 

4.4804 

-23.0407 

17 

1.0631 

-24.5290 

17 

5.3454 

-24.4906 

18 

0.7750 

-21.9325 

18 

3.9021 

-21.9117 

19 

0.9262 

-23.3889 

19 

4.6600 

-23.3595 

1  Template  Data  for  v  =0.0100  [rad/scc]  | 

1  Template  Data  for  u  =0.0500  [rad/sec]  | 

Flight  Condition 

lilESEISI 

Magnitude  [dB] 

1 

3.6125 

-15.3225 

1 

17.3525 

-14.8924 

2 

4.3238 

-16.7762 

2 

20.5000 

-16.1822 

3 

3.1379 

-14.1927 

3 

15.1896 

-13.8597 

4 

3.7609 

-15.6435 

4 

18.0193 

-15.1815 

5 

3.7353 

-15.8505 

5 

17.9091 

-15.3936 

6 

4.4704 

-17.3043 

6 

21.1342 

-16.6743 

7 

3.2464 

-14.7203 

7 

15.6888 

-14.3663 

8 

3.8894 

-16.1715 

8 

18.5929 

-15.6809 

9 

6.3348 

-20.2136 

9 

28.6393 

-19.0360 

10 

5.0066 

-18.2159 

10 

23.3809 

-17.4366 

11 

4.0422 

-16.4138 

11 

19.2615 

-15.8845 

12 

9.1615 

-23.4791 

12 

38.2667 

-21.3206 

13 

10.9208 

-24.8933 

13 

43.18.59 

-22.1108 

14 

8.01 12 

-22.3705 

14 

34..5473 

-20.6231 

15 

9.5439 

-23.7911 

15 

39  3578 

-21  5059 

16 

8.8867 

-22.9571 

16 

37  3568 

-20.9022 

17 

10.5748 

-24.3735 

17 

42.2573 

-21.7151 

18 

7.7515 

-21.8475 

18 

33  6687 

-20.1881 

19 

9.2384 

-23.2695 

19 

38.4413 

-21.0923 

I)-13 


Template  Data  for  v  =0, 


Flight  Condition  Phase 


31.5686 


36.2511 


28.1424 


32.5944 


32.4254 


37.1538 


28.9510 


33.4632 


46.5126 


40.1237 


34.4223 


56.0888 


60.3077 


52.5899 


57  0844 


55.2440 


59.5315 


51.7008 


56.2533 


1000  [rad/sec 


Magnitude  [dB 


-13.7740 


-14.7266 


-12.9609 


-13.9945 


-14.2175 


-15.1493 


-13.4187 


-14.4336 


-16.6274 


-15.6408 


-14.5571 


-17.8437 


-18.1627 


-17.5339 


-17.9255 


-17.5162 


-17.8519 


-17.1913 


-17.6018 


Template  Data  for  o  =0.; 


49.1602 


53.8154 


45.4068 


50.2399 


50.0614 


54.6649 


46.3338 


51.1317 


61.7122 


56.9525 


51.9446 


67.6055 


70.2268 


65.3690 


68.2905 


67.0323 


69.7219 


64.7378 


67.7320 


Magnitude 


-10.9681 


-11.4245 


-10.5375 


-11.0810 


-11.3195 


-11.7559 


-10.9058 


-11.4278 


-12.2894 


-11.9142 


-11.4312 


-12.7236 


-12.7969 


-12.6598 


-12.7585 


-12.4497 


-12.5288 


-12.3798 


-12.4858 


Template  Data  for  v  =0.3000  (rad/sec 


Plight  Condition  Phase  [deg]  Magnitude 


57.9010 


61.6964 


54.7076 


58.8194 


58.6642 


62.3829 


59.5665 


66.8205 


63.6957 


60.0334 


70.1911 


72.0605 


68.6591 


70.7566 


69.7441 


71.6637 


68.1667 


70.3200 


-9.1816 


-9.1946 


-9.2002 


-9.2231 


Template  Data  for  i;  =0.4000  [rad/se 


Phase 


61.9815 


65.1162 


59.3102 


62.7694 


62.6313 


65.6911 


60.0206 


63.4031 


68.2831 


66.2579 


63.5770 


69.9954 


71.5003 


68.8605 


70.5446 


69.6012 


71.1444 


68.4301 


70.1575 


Template  Data  for  v 


Flight  Condition  |  Phase 


1 


46.0750 


47.1086 


46.2437 


47.3530 


46.4361 


47.4599 


46.6258 


47.7219 


9.3879 


33.7625 


42.8217 


-54.2326 


-54.9231 


-52.5977 


-53.2480 


:2  [rad/sec 


Magnitude  [dB] 


8.8178 


8.8735 


8.4541 


8.S099 


8.6140 


8.6696 


8.2438 


8.2996 


12.3616 


10.3514 


9.0212 


8.3539 


8.4471 


9.7639 


9.9042 


8.6657 


8.7612 


10.0101 


10.1519 


Template  Data  for  v  =3  [rad/sec 


Plight  Condition  Phe^sc  [deg]  Magnitude 


15.5719  14.6874 


16.5100  14.8254 


19.7053  n.8156 


20.6973  13.9397 


15.6734  14.5875 


16.6215  14.7285 


19.9383  13.6904 


20.9376  13.8164 


-76.1418  10.6386 


-33.9544  15.4484 


4.1770  15.0503 


-102.5524  1.1037 


-102.8806  1.1169 


-98.8873  2.5823 


-99.2815  2.6019 


-101.9320  1.4763 


-102.2548  1.4900 


-98.1460  2.9490 


-98.5334  2.9693 


Template  Data  for  v  =4 


Flight  Condition  |  Phase 


1 


-30.9992 


-30.8354 


-106.5135 


-91.2245 


-117.5578 


-117.7539 


-115.4520 


-115.6841 


rad/sec 


Magnitude 


17.7938 


18.0143 


17.6060 


17.8372 


17.6795 


17.9043 


17.5434 


17  7819 


5.1702 


11.0165 


15.9040 


-3.1706 


-3.1659 


-1.8323 


-1.82.57 


-2.8003 


-2.7955 


-1.4634 


-1.4565 


Template  Data  for  v  =5 


Flight  Condition  j  Phase 


-93.714 


2  -94.3744 


3 


■84.6571 


121.4210 


113.6981 


100.6753 


128.5814 


128.7271 


127.1452 


127.3156 


128.2969 


128.4392 


126.8143 


126.9809 


rad/sec 


Magnitude 


14.1949 


14.2794 


15.6128 


15.7402 


13.8715 


13.9538 


15.3476 


15.4732 


1.4830 


6.6889 


11.6505 


-5.8872 


-5.8847 


1  Template  Data  for  v  =100  [rad/sec]  || 

Flight  Condition 

Magnitude  [dB] 

Flight  Condition 

Magnitude  [dB]  {{ 

1 

-240.1843 

-29.0456 

1 

-277.8020 

2 

-240.2155 

-29.0455 

2 

-277.8200 

-41.3494 

3 

-239.8716 

-27.8739 

3 

-277.6208 

-40.1827 

1  4 

-239.9073 

-27.8737 

4 

-277.6415 

-40.1826 

5 

-240.2493 

-29.4089 

5 

-277.8397 

-41.7125 

6 

-240.2814 

-29.4088 

6 

-277.8583 

-41.7124  1 

7 

-239.9460 

-28.2372 

7 

-277.6640 

8 

-239.9828 

-28.2370 

8 

-277.6853 

'f 

-240.9892 

-38.1403 

9 

-278.2683 

-50.4198  1 

t.O 

-240.6313 

-34.0878 

10 

-278.0603 

11 

-240.2273 

-30.6146 

11 

-277.8262 

12 

-241.4905 

12 

-278.5605 

-57.1812 

13 

-241.5035 

■KSKHi 

13 

-278.5681 

-57.1812 

14 

-241.3672 

■ES&QSHI 

14 

-278.4886 

-56.0147 

15 

-241.3822 

15 

-278.4973 

-56.0147 

16 

-241.4640 

16 

-278.5450 

-56.8183 

17 

-241.4766 

-44.5460  1 

17 

-278.5524 

-56.8183 

18 

-241.3369 

[■KKaHI 

18 

-278.4709 

-55.6518 

19 

-241.3514 

19 

-278.4793 

D-18 


Appendix  E.  A  Problem  With  MATRIX^  and  QFT 
Template  Generation 

A  problem  with  MATRIX^  arose  during  the  longitudinal  SISO  QFT  design. 
Eight  of  the  nineteen  plants  used  in  the  design  contain  unstable  poles,  while  the 
others  are  completely  stable.  Having  both  stable  and  unstable  plants  uncovered  a 
problem  with  the  way  MATRIXjf  handles  phase. 

An  apparent  problem  with  the  MATRIX^  Bode  plotting  routine  exists  if  used 
to  generate  data  for  plant  templates  necessary  for  a  QFT  design.  If  the  phase  of  the 
plant  at  the  lowest  plotted  frequency  exceeds  180°  or  is  less  than  -180°  ,  the  program 
removes  phase  information  from  the  Bode  data  so  that  the  initial  phase  data  will  be 
between  ±  180°  .  When  Bode  plots  of  stable  and  unstable  plants  are  used  together, 
as  in  forming  QFT  templates,  the  Bode  data  indicates  a  minimum  phase  difference 
for  low  input  frequencies,  and  360°  difference  at  higher  frequencies.  Note,  in  general, 
w'-domain  transfer  functions  are  equal  order  over  equal  order.  Thus,  analysis  of  the 
plants  at  very  low  frequencies,  that  is  by  taking  the  limit  as  the  frequency  approaches 
0,  predicts  a  wide  variation  in  phase.  Further,  by  taking  the  limit  as  the  frequency 
approaches  infinity,  minimum  phase  difference  should  exist  for  the  higher  frequencies 
(for  positive  gain,  an  angle  of  0°  must  result,  for  negative  gain,  an  angle  of  180°  must 
result). 

An  experiment  using  simple  plants  illustrates  the  problem.  Plants  Pi(5)  = 

and  p2(s)  =  are  entered  and  plotted.  MATRIX^  produces  the  expected 
results.  See  Figure  E.l.  The  Bode  data  for  pi(s)  produces  0°  phase  at  very  low 
frequencies  and  —90°  phase  at  very  high  frequencies.  The  unstable  plant  p2  produces 
a  -180°  phase  at  very  low  frequencies  and  -90°  phase  at  very  high  frequencies,  also 
as  expected.  When  the  number  of  poles  in  the  plants  are  increased,  the  results 
are  not  consistent  with  predictions.  p3(s)  =  and  p4(s)  =  are  entered 


E-1 


Figure  E.l.  Comparison  of  a  Single  Stable  Pole  and  a  Single  Unstable  Pole 

and  plotted.  MATRIX;^  produces  Bode  phase  data  for  the  stable  plant,  ps  exactly 
as  predicted  by  mathematical  analysis.  However,  MATRIX^  produces  Bode  phase 
data  for  the  unstable  plant  p4  with  360°  subtracted  from  the  predicted  values  (see 
Figure  E.2).  For  example,  limu,-,o  =  —360°,  and  lim,„_»oo  ^PaU^^)  =  —180°. 

Normally,  the  fact  that  MATRIXjr  starts  the  Bode  phase  data  between  180  and 
-180°  phase  is  not  a  problem  but  for  QFT  template  generation,  the  lost  phase  can 
greatly  change  the  results.  For  the  above  examples,  the  templates  formed  from  the 
Bode  data  of  the  single  pole  case,  pi  and  p2  can  be  used  for  QFT  synthesis.  However, 
in  the  double  pole  case,  the  MATRIX^  Bode  phase  data  will  produce  very  narrow 
templates  for  low  frequencies  and  wide  templates  for  the  higher  frequencies.  The 
phase  data  calculated  analytically,  without  MATRIXa^  ,  produces  narrow  templates 
at  the  high  frequencies. 

One  solution  is  to  calculate  the  initial  phase  of  the  plants  by  taking  the  limu,-,o- 
The  Bode  data  taken  from  MATRIX^  is  then  adjusted  by  the  proper  amount  before 
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Figure  E.2.  Comparison  of  Two  Stable  Poles  and  Two  Unstable  Poles 

the  data  is  used  to  generate  the  templates.  Another  solution,  which  is  applicable 
for  transfer  factions  of  equal  numerator  and  denominator  order  like  those  obtained 
from  bilinear  transformations,  is  to  obtain  the  frequency  response  data  by  plotting 
the  Bode  plots  “backward”  in  frequency.  For  example,  plotting  P4(s)  from  100 
(rad/sec)to  .1  (rad/sec)results  in  the  correct  phase  information  and  may  be  used  in 
forming  templates  without  adjusting  the  data  (see  Figure  E.3). 
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Figure  E.3.  Comparison  of  Two  Stable  Poles  and  Two  Unstable  Poles  After  Plot¬ 
ting  “Backward”  in  Frequency 
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Appendix  F.  Difference  Equations 

This  appendix  includes  the  difference  equations  needed  for  digital  implemen¬ 
tation  of  the  QFT  rate  controller  on  Lambda  assuming  a  sample  rate  of  50  Hz. 

Equations  (3.42)  and  (3.43)  are  z-plane  representations  of  the  QFT  prefilter 
and  compensator  for  the  Longitudinal  SISO  digital  QFT  design  and  Equations  (4.23), 
(4.24),  (4.21),  and  (4.22)  are  the  z-plane  representations  for  the  lateral-directional 
MIMO  digital  design.  The  z-plane  representations  are  converted  to  difference  equa¬ 
tions  for  implementation  on  the  Lambda  flight  control  computer  and  are  included  in 
Tables  F.l  and  F.2. 

The  z-plane  transfer  functions  are  of  the  form: 

-I-  •  •  •  +  CiZ  -f  Co  - 

U {z)  dnZ^  -t-  -1-  •  •  •  -f  diz  -f  do 

where:  n  >  m  and  dn  =  1.  The  difference  equation  is  formed  by  dividing  the 
equation  by  2",  cross  multiplying,  and  replacing  the  negative  powers  of  z  with  the 
appropiate  delay  to  obtain: 

y{kT)  =  -d„_xi/p-l)r] - d,y[{k^-l-n)T]-doy[{k-n)T] 

+Cmu[{k  +  m-  n)r]  +  Cjn-iu[{k  +  m-n  -  1)T]  +  ■•■ 

-f-Ciu[(A;  -I- 1  -  7i)T]  +  com[(^  -  Ti)T]  (F.2) 

fn{kT)  =  1.617064707083078/„[(/t-l)r] 

-0.648907940395999/u  [{k  -  2)T] 
-f-0.017911818738518^e,„rfP)T] 

+0.01 194 1212492345^ecmd[(^  ~  1)^] 

+0.001 9902020820585e,„,  [{k-  2)7’]  (F.3) 
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f22{kT)  =  0.960784313725490/22[(A-l)r] 

+0.084967320261435a,„,  [{k)T] 
-0.0457516339869285„,„, 

hsikT)  =  1.986073568578539/33[(ifc-l)r] 

-0.98609740160436/33  [(Jfc  -  2)] 
+0.014907557651096^,,^,  [{k)T] 
-0.014883724625275^,,^, [(it  -  l)r] 
-0.00^„„,[(fc  -  2)T] 


gn{kT)  =  0.768292682926834i?„[(it  -  l)r] 

+0.878048780487800^11  [(it  -  2)T] 
-0.170731707317077^11  [(it  -  Z)T] 
-0.365853658536584^^11  [(it  -  4)2’] 
-0.109756097560974511  [(it  -  5)7] 
+96.0079664634145ui  [{k)T] 
-371.4464908536588ui  [(it  - 1)7] 
+562.6357134146366U,  [(it  -  2)7] 
-414.0391280487847U1  [(it  -  3)7] 
+146.2375884146371U1  [(it  -  4)7] 
-19.3892103658544ui[(/t  -  5)7] 


922{kT)  =  2.473684210526316522[(it-l)7] 

-1.789473684210527522  [(it  -  2)7] 


+0.526315789473684522  -  3)T] 

-0.578947368421052522  [(A:  -  4)T] 
+0.368421052631579522  P  -  5)T] 
+021.7235319457573U2  [{k)T] 
-100.4706471025779u2[{fc  -  l)r] 
+186.6047527121375u2[(fe  -  2)T] 
-174.1142948925886u2[(ib  -  Z)T] 
+81.6824678490870u2[(ifc  -  4)T] 
-15.4257877749731u2[(&  -  5)r]  (F.7) 

gzzikT)  =  1.88323500491642l533[(fc-l)r] 

-0.719354310062275533  [(ife  -  2)T] 
-0.244428056374959533  [(A:  -  d)T] 
+0.08054736152081 3533  [(Ar  -  A)T] 

+0.00533[(&-5)r] 

+37.4838623948432U3  [( A:)r] 
-142.5540064277651u3[(A:  -  1)7] 
+214.1271236206708u3[(A:  -  2)T] 
-158.6615822681088u3[(A:  -  3)7] 
+57.9432763392695u3[(A:  -  4)7] 
-8.3349961761171u3p-5)7]  (F.8) 
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Table  F.l.  QFT  Prefilter  Difference  Equations 


Coefficients  of 

■BOH 

■m 

^33 

Equation  (F.2) 

y(kT)= 

Cm 

+0.013975 

+0.084967 

+0.014908 

Cm— 1 

+0.0093168 

-0.045752 

-0.014884 

Cm-2 

+0.0015578 

*^n-l 

-1.6439 

-0.96078 

-1.9861 

dn-2 

+0.66874 

+0.98610 

Table  F.2.  QFT  Compensator  Difference  Equations 


Coefficients  of 
Equation  (F.2) 

II 

o22 

y(kT)= 

g33 

y(kT)= 

Cm 

+96.749 

+21.724 

+37.484 

Cm— 1 

-374.31 

-100.47 

-142.55 

Cm— 2 

+566.96 

+186.60 

+214.13 

Cm— 3 

-417.22 

-174.11 

-158.66 

Cm— 4 

+147.36 

+81.682 

+57.943 

Cm— 5 

-19.538 

-15.426 

-8.3350 

^n— 1 

-0.78629 

-2.4737 

-1.8832 

dn-2 

-0.37805 

+1.7895 

+0.71935 

<ln-.3 

+0.17073 

-0.52632 

-i-0.24443 

dn— 4 

+0.36585 

+0.57895 

-0.080547 

dn-5 

+0.10976 

-0.36842 

0 
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Figure  F.l.  Block  Diagram  of  Rate  Controller  QFT  Design  for  Lambda 
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